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ABSTRACT 


This report describes the preliminary design specifications 
for an Advanced Space Transportation System consisting of a fully 
reusable flyback booster, an intermediate-orbit cargo vehicle, 
and a shuttle-type orbiter with an enlarged cargo bay. It 
provides a comprehensive overview of mission profile, 
aerodynamics, structural design, and cost analyses. These areas 
are related to the overall feasibility and usefulness of the 


proposed system. 
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INTRODUCTION 


The support of the United States’ future space station and 
further space ventures prompts the need for an advanced space 
transportation system which can transport cargo or passengers 
efficiently, safely, and economically. The existing Space 
Transportation System consists of three reusable shuttle orbiters 
fueled by an expendable external tank and boosted by two reusable 
solid rocket boosters. While this system is adequate for our 

present space needs, an increase in mission frequency, mission 
diversity, and payload quantity requires a more economical and 
flexible system with a greater payload capability. 

The boosters used in the present system descend into the 

ocean after staging. The expense of the retrieval and 

refurbishment of these boosters, along with the cost of replacing 

the external tank, adds substantially to the cost of each 

mission. When the number of missions is increased in the near 

future, lowering these expenses will be of utmost importance. A 

winged booster capable of flying back under its own power could 

be used for a number of different missions, thus cutting down on 

* 

the cost of each mission. 

With the construction and operation of a space station, a 
larger cargo capability will be needed than is available on the 
present shuttle. The creation of a new, unmanned cargo vehicle 
would substantially increase the quantity of cargo that could be 
transported to the station in one mission. The ability of the 
system to perform such "shipping" missions also increases its 


flexibility. Also* by designing a vehicle capable of 
transporting cargo to an intermediate orbit? "orbital warehouses" 
would be created that could store supplies for later use. 

With increased space activities? a larger manned orbiter 
with a larger return payload capability will be needed. An 
enlarged shuttle-type vehicle with a 15’ X 60’ cargo bay would 
improve the United States’ ability to perform passenger transport 
and repair missions. This vehicle could be boosted into orbit by 
the same flyback booster used to boost the cargo vehicle? thus 
doubling the role of the winged booster. 

The Advanced Space Transportation System described above 
would more readily fulfill the needs of America’s space program 
in the future than the existing system. The technology required 
to construct and operate this system is readily available today. 
Also? the savings incurred by this system? along with the 
increased compatibility with the needs of future space projects? 
would more than make up for the its cost. 




GROUNDRULES FOR FLYBACK BOOSTER STUDY 
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1) Two vehicles are to be developed: 

* Manned shuttle II 

* Unmanned carqo vehicle 

2) Both vehicle will be S-staqe vehicles utilizing the same 
fully reusable* winged * flyback booster. 

3) The payload capability for each vehicle will be as 
follows: 

* Shuttle 11-40 Klb. to space station orbit 

(39.5° inclination / 270 n.m. altitude) 
with 40 Klb. return to earth 

* Carqo vehicle - 125-150 Klb. to 2B.5 0 / 150 n.m. 

4) The payload envelope <bay) for each. vehicle is 

* Shuttle II - 15 ’ D X 60’L 

* Cargo vehicle - 25’D X 90 ! 'L and 33’D X 100'L 

5) Engine propellants: 

* l** stage - liquid oxygen / hydrocarbon 

* S"* 1 stage - liquid oxygen / liquid hydrogen 

6) Staging velocity - maximum of 7 * 000 fps 

7) The manned shuttle II will see a maximum of 3 g’s 

8) The carqo vehicle will be ready for its first flight in 
1998 therefore technology avid design freeze approximately 
1990 

9) The manned shuttle II will be ready for its first flight 
in 2005; therefore technology and design freeze 
approximately 1997 
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TRAJECTORY ANALYSIS 


The following analysis examines the complete trajectory 
of the flyback booster. In a general overview, the trajectory 
is broken into two phases. These phases include the 
trajectory to staging and the return trajectory. 

Trajectory to Staqinq 

The preliminary design of the Advanced Space 
Transportation System tASTS* is largely contingent on data 
related to the system’s trajectory to staging. Trajectory 
analysis allows the designer to generate various parameters 
essential for estimating the system’s thrust, mass and 
structural requirements. The most important of these 
parameters are those associated with the data gathered at 
staging. This data consists of the mass, velocity, altitude 
and time to staging. From this data the required propulsion 
parameters can be determined. It is then possible to calculate 
the mass and volume of the propellant burned. Ultimately, the 
dimensions of the structure can be computed from the size 
requirements of the fuel avid oxidizer tanks. 

As in any preliminary design project, it is necessary for 
the designer to follow certain groundrules and to make various 
assumptions. The groundrules and assumptions allow the 
designer to generate data that can then be studied for 
validity. For the ASTS, the groundrules are listed on page 3 


and the assumptions pertaining to the trajectory analysis are 
listed below: 

♦T/W— 1 .35 

*3 G’s acceleration limit 
♦Gravity turn into orbit 
♦Inertial and relative frames used 
♦Drag neglected during boost phase 
♦Parallel burn 

♦Crossfeed propellant from booster to 2nd stage 
♦Thrust of both stagey is constant 
♦Booster characteristics 
-L0x/JP4 propellant 
-Isp=320 sec 
-LGx/JP4 ratio=2.3 
♦End Stage characteristics 
-L0>:/L.H2 propellant 
-Isp=380 sec 
-L0x/LH2 ratio-=6.0 

Some of the assumptions require further examination. The 
thrust to weiqht ratio was chosen to help maintain a G level 
below 3.0 due to the human element involved. The parallel 
burn allows the system to achieve the required thrust level 
without putting the entire thrust demand on the booster. The 
propellant crossfeed system eliminates excess structural 
weight which the second stage would otherwise have to carry 
into orbit. Finally? the hydrocarbon fuel JP^ was chosen 



because it can be used in the liquid rocket engines as well as 
the air breathing engines. 

The G turn is used to aid the second stage in achieving 
the required relative angle of 90 degrees for a circular orbit 
at a particular altitude. Figure 1 shows how the ASTS is 
tracked in a relative and inertial reference frame during its 
flight. Notice how the range and altitude are measured. 

Coupled with the assumptions and groundrules , a variety 
of trajectory, propulsion and mass equations are incorporated 
into a BASIC program. Program 1, listed in Appendix B, is 
similar to the one used in the preliminary report, however the 
G turn is incorporated. A time increment of one second is 
used to determine all the parameters during the flight. The 
velocity restriction of 7000 feet per second at staging, 
together with the time increment, enables the designer to 
determine the time at which staging occurs. 

Three variables are entered at the start of the program; 

the total mass of the ASTS, the altitude for the G turn to 

begin, and the liftoff thrust of the second stage. From this 

data the total thrust and the thrust of the booster is 

dictated. Dividing the thrust of each stage by its respective 

specific impulse yields the mass flow of each stage. In order 

to determine the mass flow for the ASTS, an equivalent 

specific impulse is calculated by., 

Ispeq=T/(MD0T 1 +MD0T a ) (1) 

then the total mass flow is given by, 

L 






MOOT®*! / Ispeq 


( 2 ) 
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In order to calculate the change in velocity during each time 
interval, the new total mass at each interval must first be 
calculated by, 

DM=M-MDOT i *DT-MDOT e *DT (3) 

Using this result, the change in velocity can be determined 
by, 

DV=Ispeq*G*log(M/DM)-G*DT (4) 

where G is a function of altitude and the inertial angle 0 . 
Once the change in velocity has been determined, the altitude 
can be found using, 

Y=Y+.5*<2v+DV)*DT (5) 

The range is given by the relationship, 

RANG=RE*e <6> 

where RE is the radius of the earth. 

By using the above relations in a stepwise time interval 
loop, the parameters at staging <V=7000 fps) can be 
determined. From the time to staging parameter (TB) it is 
possible to find the total mass of propellant burned to 
staging by the equation, 

Mp=MDQT*TB (7) 

It follows that the total mass at staging is given by, 

Mb=M-Mp (8) 

It is also possible to examine the stages individually. 
The respective propellant masses are given as follows: 


Mp»-MDOTi*TB 

7 


< 9 ) 


Mp a =MDOT e *TB 


( 10 ) 


where the subscripts 1 and S correspond to the stage number. 
The masses of the fuels and oxidizer are then calculated with 
the use of the equation? 

Mp=M0x+Mf (11) 

together with the oxidizer to fuel ratios of 2.3 for the 
booster and 6.0 for the second stage. Putting these relations 
together yields: 


Mf t »Mpt/3.3 

( 12) 

M0x 1 =Mpi— Mf i 

( 13) 

Mf*=Mp«/7.0 

( 14) 

MOxe = MpE”Mfe 

( 15) 


It follows that the volumes of the propellants can be 
determined by dividing the various propellant masses by their 
respective densities. 

Using Program 1 and the above relations the results 
listed in Table 1 were obtained for the booster with shuttle 
II configuration. A similar running of Program 1 for the 
booster with the cargo vehicle yielded the results in Table 2. 
It is evident that for both configurations? the structural 
weight of the booster remains a constant. Due to the 
extremely small density values encountered by the ASTS at 
staging? some problems for the return trajectory of the 
booster were created. This problem could probably be solved 
by incorporating a perigee injection trajectory. This is a 
problem to be examined in further detail in the future. 


FVLBACK BOOSTER WITH SHUTTLE II 


GENERAL INFORMATION 

DENSITY OF OXIDIZER <O e ) = 71. E lbm/ft 
DENSITY OF SHUTTLE II FUEL <H*> = 4.4 lbm/ft 
DENSITY OF BOOSTER FUEL<JP4> = 47 . E99 lbm/ft 

ISP OF SHUTTLE II - 380 sec 
ISP OF BOOSTER = 320 sec' 

EQUIVILENT ISP = 325. 1 sec 

OXIDIZER TO FUEL RATIO FOR SHUTTLE II = 6.0 
OXIDIZER TO FUEL RATIO FOR BOOSTER = E.3 

TOTAL MASS ON LAUNCH PAD = 7,400,000 lbf 

TOTAL THRUST OF BOOSTER AND SHUTTLE II = 9,990,000 lbf 

THRUST OF SHUTTLE II = 1,000,000 lbf 

THRUST OF BOOSTER = 8,990,000 lbf 

ALTITUDE OF G-TURN = 2,6£5 ft 


STAGING DATA 

TIME TO STAGING = 167 sec 
VELOCITY AT STAGING ■ 7,076 ft/sec 
ALTITUDE AT STAGING * 344,676 ft 

MASS OF SHUTTLE II AT STAGING = 1,E9E,850 lbm 
MASS OF BOOSTER AT STAGING 

(STRUCTURE + RETURN FUEL) * 1,112,810 lbm 

MASS FLOW OF SHUTTLE II = 2,623 lbm/sec 
MASS FLOW OF BOOSTER = 28,094 lbm/sec 

MASS OF BOOSTER PROPELLANT BURNED = 4,691,698 lbm 
MASS OF SHUTTLE II PROPELLANT BURNED = 439,544 lbm 
MASS OF BOOSTER OXIDIZER BURNED = 3,269,971 lbm 
MASS OF BOOSTER FUEL BURNED = 1,421,727 lbm 
MASS OF SHUTTLE II OXIDIZER BURNED = 376,752 lbm 
MASS OF SHUTTLE II FUEL BURNED = 62,792 lbm 


TABLE 1. Shuttle II Configuration 
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FVLBACK BOOSTER WITH CftRGQ VEHICLE 


GENERAL INFORMATION 

DENSITY OF OXIDIZER (O a > = 71. E lbm/ft 
DENSITY OF CARGO VEHICLE FUEL < H a > = 4.4 lbm/ft 
DENSITY OF BOOSTER FUEL<JP4) = 47.E99 lbm/ft 

ISP OF CARGO VEHICLE = 380 sec 
ISP OF BOOSTER = 3E0 sec 
EQUIVILENT ISP = 3E7.7 sec 

OXIDIZER TO FUEL RATIO FOR CARGO VEHICLE =6.0 

OXIDIZER TO FUEL RATIO FOR BOOSTER =3.3 

TOTAL MASS ON LAUNCH PAD = 7,571,000 lbf 

TOTAL THRUST OF BOOSTER AND CARGO VEHICLE = 10,EE0,850 lbf 
THRUST OF CARGO VEHICLE = 1,531,000 lbf 
THRUST OF BOOSTER = 8,689,850 lbf 

ALTITUDE OF G-TURN = 1,800 ft 


STAGING DATA 

P 

TIME TO STAGING = 164 sec 
VELOCITY AT STAGING = 6,987 ft/sec 
ALTITUDE AT STAGING = E81,390 ft 

MASS OF CARGO VEHICLE AT STAGING - 1,491,380 lbm 
MASS OF BOOSTER AT STAGING 

(STRUCTURE + RETURN FUEL) = 1,113,810 lbm 

MASS FLOW OF CARGO VEHICLE = 4,039 lbm/sec 
MASS FLOW OF BOOSTER = 37,156 lbm/sec 

MASS OF BOOSTER PROPELLANT BURNED = 4,453,584 lbm 
MASS OF CARGO VEHICLE PROPELLANT BURNED = 660,756 lbm 
MASS OF BOOSTER OXIDIZER BURNED = 3,104,013 lbm 
MASS OF BOOSTER FUEL BURNED = 1,349,571 lbm 
MASS OF CARGO VEHICLE OXIDIZER BURNED = 566,363 lbm 
MASS OF CARGO VEHICLE FUEL BURNED = 94,394 lbm 


TABLE 3. Cargo vehicle Configuration 



Notice that the staging parameters for both configurations are 
quite similar. This makes the next phase of the trajectory 
analysis much easier. The return trajectory of the booster 
will be examined for the shuttle configuration only. However 5 
the results for both configurations are listed for comparision 
in Tables 4 and 5 on page 24 and 25. Again the results of 
this analysis would be similar for both configurations. 


Return Trajectory 

The return trajectory of the flyback booster was 
simulated using the fortran program DUAL? listed in Appendix B 
as Program # 2. The program was acquired by Auburn University 
from Marshall Space Flight Center. This very powerful program 
is outlined in the Users Manual? which accompanies this 
report. 

As in the trajectory analysis to staging? certain 
groundrules and assumptions must be examined during this phase 
of the flight. The assumptions are as follows: 

*95 percentile winds at 280 deg azimuth 
*MFSC flyback booster aerodynamics 
*Wing reference area=15500 sq ft 
♦Cruise engine start at 20000 ft 
♦High cruise alti tude= 10000 ft 
♦Low cruise altitude 5 * 1000 ft 


*500 ft /min descent to 1000 ft 

(U 


♦End trajectory within 1 N.MI. of KSC 
♦Maximum angle of attack=70 deg 
♦Cruiseback angle of attack=& deg 

♦Sea level thrust per air breathing engine=62500 lb 
♦Engine cantor --6 deg 
♦Launch/Landing at KSC 

The MSFC aerodynamics are used because the configuration of 
the booster is very similar to MSFC’s booster. The wing 
reference area is determined such that the wing loading does 
not exceed 70 pounds per square foot. The maximum angle of 
attack is related to the angle at which the most drag is 
experienced by the booster during reentry. Finally? The 
cruiseback angle of attack is determined by the angle which 
creates the maximum lift to drag ratio. 

The staging data which is of importance in the program 
DUAL is as follows: 

♦Booster weight (including fuel )=1 ? 1 12? BIO lbm 

♦GCLAT=28 . 46 deg L0NG=79.81 deg 

♦Relative veloc i ty=7076 ft/sec 

♦Alt i tude=344676 ft 

♦Relative azimuth=93.91 deg 

♦Relative flight path angle=53.6 deg 

♦TB=167 sec 

All of this data is input into the program DUAL through an 
input file called ZFLY.DAT which is listed in Appendix B. 
Once this data has been input? the program can run. In this 



manner the entire trajectory from staging to the return to KSC 
can be examined. 

The return trajectory is broken into eight phases by the 
program DUAL. These phases are outlined in Figure 2. The 
different phases after separation include coast to apogee? 
coast to reentry? altitude descent rate control? 45 degree 
bank turn? engine start? idle-engine descent? high cruise? 
letdown to low cruise? low cruise and arrival. During 
reentry? the altitude control phase begins when the dynamic 
pressure is sufficient to obtain a constant rate of descent of 
—100 feet per second. This is done by modulating the bank 
angle. The constant rate of descent down to low cruise is 
obtained by throttling back the cruise engines. Also included 
in this analysis is the capability of flying back with one 
engine out. 

The program DUAL had to be run numerous times before 
acceptable results were obtained. The problems encountered 
included excessive G limits? inability to achieve high cruise? 
insufficient thrust levels? and high Mach number when engines 
are started. All of these problems were solved by adjusting 
the angles of attack? imposing a dynamic pressure limit? 
increasing the thrust level? and lowering the altitude when 
the engines are started. The results of the program are 
listed in Appendix B. Table 3 shows the values of various 
parameters during the eight phases of the return trajectory. 



FIGURE 2. FLIGHT PROFILE 



SHUTTLE II MISSION 



PHASE 

t 

sec 

h 

ft 

V 

fps 

M 

a 

deq 

a 

deq 


Separation 

167 

344,676 

7076 

7.0 

0 . 0 

0 - 0 

1 

Apogee 

367 

905,151 

4019 

1.47 

70 . 0 

0 - 0 

a 

Alt control 

6E9 

34, 10E 

E683 

S . 69 

0.0 

0 . 0 

3 

begin turn 

637 

£7,771 

1799 

1.75 

0 . 0 

£8 . 7 

4 

engines on 

667 

SO , 049 

797 

0 . 75 

6.0 

45.0 

5 

high cruise 

839 

10,113 

47E 

0 . 43 

6.0 

0 . 0 

6 

letdown 

5889 

5 , 5E6 

4E0 

0.38 

6.0 

0 . 0 

7 

low cruise 

6483 

1 , 060 

391 

0 . 35 

6.0 

0 . 0 

8 

arrival 

6738 

933 

389 

0 . 34 

6.0 • 

1 .£ 


Table 3 


Return trajectory data for booster 



It is important to examine the maximum values of some of the 
parameters throughout the flight. The parameters of interest 
are as follows: 


Max 

QDOT 

89.488 

btu/ft 2 /s 

TIME 

613 

sec 

Max 

DYNP 

3674 

lb/ft 2 

TIME 

681 

sec 

Max 

ACC 

8.131 

g’s 

TIME 

603 

sec 

Max 

T 

805934 

tr 

TIME 

841*5 

f sec 

TOTAL FLIGHT TIME 

1 hr 

53 min 

18 

sec 


The time specificed for each parameter above is from launch to 
when the event occurs. It is evident that the dynamic 
pressure and the G level are both guite high. This is a 
result of the size of the booster and the high altitude at 
which staging occurs. The booster accelerates very quickly as 
it reenters the atmosphere and experiences large forces as it 
tries to deccelerate before entering the bank turn. As the 
designer becomes more familiar with the program DUAL these 
maximum values may be reduced by imposing additional 
limitations. Regarding the thrust requirements, the program 
was used to size the engines and determine the number of 
engines required for flight. Four 6E500 pound engines are 
utilized. The engines reach a maximum throttle of 85 percent, 
and can therefore maintain an engine out capability if needed. 
The total fuel requirement for the four engines is 136866 
pounds of JP4. 

To achieve an overview of the entire flight, five graphs 
are included, which plot various parameters versus time. 


Figure 3 shows altitude versus time. Notice* that the booster 
reaches an altitude of 905151 feet at apogee*. Fiqure 4 shows 
the mach number versus time. The mach number increases from 
apogee but then is decrease during the controlled rate of 
descent. Figure 5 shows the acceleration profile throughout 
the flight in terms of G’s. Figure 6 is a plot of the dynamic 
pressure during the flight. Notice that it is a maximum 
during the controlled rate of descent and then decreases to 
nearly a constant. Finally? Figure 7 shows the thrust levels 
achieved versus time. 

In conclusion? it should be evident that the booster can 
be flown back to its launch site? assuming that the structure 
can be built to withstand the large forces it will encounter 
during the critical controlled rate of descent. The booster 
will not be manned as originally planned? due to the high 
accelerations experienced after separation. 
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SHUTTLE II MISSION 



PHASE 

t 

sec 

h 

ft 

V 

f DS 

M 

« 

deq 

|3 

deq 


Separation 

167 

344,676 

7076 

7 . 0 

0 . 0 

0.0 

1 

Apogee 

367 

905,151 

4019 

1 .47 

70 . 0 

0.0 

a 

Alt control 

689 

34,108 

8683 

8.69 

0 . 0 

0.0 

3 

begin turn 

637 

37,771 

1799 

1 .75 

0.0 

88.7 

4 

engines on 

667 

80 , 049 

797 

0 . 75 

6.0 

45 . 0 

5 

high cruise 

839 

10,113 

478 

0 . 43 

6.0 

0 . 0 

6 

letdown 

5889 

5,586 

480 

0 . 38 

6.0 

0 . 0 

7 

low cruise 

6483 

1 , 060 

391 

0 . 35 

6.0 

0 . 0 

a 

arrival 

6738 

933 

389 

0.34 

6.0 

1.8 


PHASE 

CARGO VEHICLE MISSION 
t h V M 

sec ft f ds 

<x 

deq 

(3 

deq 


Separation 

164 

881,390 

6987 

10.1 

0 . 0 

0.0 

1 

Apogee 

360 

883,557 

3973 

1.49 

70 . 0 

0.0 

a 

Alt control 

608 

33 , 848 

8766 

8.77 

0 . 0 

0.0 

3 

begin turn 

618 

85 , 703 

1660 

1 . 60 

0 . 0 

88.7 

4 

engines on 

646 

18,867 

797 

0.74 

6.0 

45 . 0 

5 

high cruise 

800 

10,055 

467 

0 . 48 

6.0 

0 . 0 

6 

letdown 

5680 

5 , 438 

481 

0.38 

6.0 

0 . 0 

7 

low cruise 

6804 

1 , 053 

393 

0 . 35 

6.0 

0 . 0 

8 

arrival 

6454 

938 

391 

0.34 

6.0 

1.8 


Table 4 : Return trajectory data for booster 

for both missions 


SHUTTLE II MISSION 


PARAMETER 
Max ACCL 
Max QDOT 
Max STMP 
Max DYNP 
Max T 

FLYBACK FUEL 
TOTAL FLIGHT 

PARAMETER 
Max ACCL 
Max QDOT 
Max STMP 
Max DYNP 
Max T 

FLYBACK FUEL 
TOTAL FLIGHT 


VALUE UNITS FLIGHT TIME OCCURRED 


a. 131 

g’s 

603 

sec 

29.428 

btu/f t 2 /s 

613 

sec 

2310 

°F 

613 

sec 

3674 

lb/ft 2 

621 

sec 

205,934 

lb-r 

2419 

sec 

136, 

866 lb m 



IME 1 hr 

53 min 

18 sec 


CARGO VEHICLE MISSION 


VALUE 

UNITS 

FLIGHT TIME OCCURRED 

7 . 278 

g’s 

594 

sec 

29.113 

btu/f t 2 /s 

592 

sec 

2303 

°F 

592 

sec 

3764 

lb/ft 2 

602 

sec 

205 , 978 

lb-r 

2400 

sec 


131,266 lb m 

TIME 1 hr 47 min 34 sec 


Table 5 


Return trajectory for booster 
for both missions 



STRUCTURES 


Due to the concept of the booster burning in parallel 
with the other two vehicles and having to carry the fuel and 
oxidizer for the same vehicles plus its own propellants during 
the ascent portion of the boost j the booster needed to be 
built around the propellant tanks. Since the booster and other 
vehicle would be using the same oxidizer this tank could be 
shared. The preliminary design dictated the use of a three 
tank booster containing L.H a , JP*. and L0 a . From the given 
data; the following tank volumes were founds 

1) LH*» = 2 1,454 Ft® 

2) JP* = 32,952 Ft® 

3) L0 a * 51,551 Ft® 

The three tanks were designed to be cylindrical tanks 
with semi-prolate spheroid end caps. Originally spherical end 
caps were being used but by incorporating the semi— prolate 
spheroids the total tank height could be reduced thus 
decreasing the overall height of the booster. The volume of 
each tank is given by eguation 1. 

V=^/3irrb® + irr«h (1) 

Where r is the radius of the inside of the tank, b is the 
length of the semi-minor axis and h is the height of the 
cylindrical part of the tank. Knowing the volumes needed in 
each tank, eguation one can be rearranged to give the overall 
tank height once a radius and semi-minor axis has been chosen. 
The total tank height H can be found using equation 2. 



( 2 ) 


H= V - /»/3ir_ rb*_ + 3b 
irr* 


Equation 2 was used to find the heights of all three 
propellant tanks. 

A radius of 18 feet was used so that the booster would be 
approximately the same diameter as the Shuttle and cargo 
vehicles. The height of the semi-minor axis was determined to 
be six feet. The tank heights were then calculated as 
follows: The height of the LH*. tank is 31 feet, the height of 
the LOs tank is 60 feet and the height of the JP*. is 42 feet. 
These tank heights were calculated using the maximum volumes 
of fuel used for any mission plus all flyback fuel needed. 

The remaining heights of the booster sections were then 
determined. The three fuel tanks are placed with the LHa tank 
in the forward position, the LO® tank in the middle position 
and the JP*, tank being in the rear position near the boosters 
main engines. The nosecone section of the booster will house 
the avionics needed to run all of the systems on .the booster. 
This section will be 33 feet long. Second, the structural 
area behind the front tank will be used as attaching points 
for the canards, the air breathing engines and the mounts for 
the accompanying vehicle. Twenty feet of structure was needed 
to incorporate all of these attaching structures. Ten feet of 
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structure was placed between the L0« and JP*. tanks for 
plumbing and insulation reasons. Forty eight feet of was need 
behind the JP*. tank for the engine thrust structure, and aft 
attaching points for the vehicles. A diameter of 45 feet was 
need around the main enqines therefor a firing was placed 
around the engines. The overall height of the booster came 
out to be 244 feet with a diameter of 36 feet and a firing of 
45 feet in diameter placed around the engine’s and aft thrust 
structure. Figure 14 shows a three view drawing of the 
booster . 

The booster system will be determined by trends avid load 
analysis used on the X — 1 5 » the Saturn and the Space Shuttle 
programs. The ground rules for the analysis stated that no 
thermal protection will be used. fhe concept of a "Heat Sink" 
will be based upon the use of new high-temperature metals. 
Initial temperatures set for a warm windy day are as follows: 

—50 Deg F Liquid Hydrogen Tank Wall 
0 Deg F LOX Tank Wall 
120 Deg F Other Surfaces 

The maximum allowable temperatures are as follows,: 

350 Deg F Aluminum Lithium 

800 Deg F Titanium (Load Structures) 

lOOO Deg F Titanium (Non-Load Structures) 

1200 Deq F Inconel X 

The basic concept of the heat sink booster is to allow 
the booster to act as a sink to store heat as it is generated 
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during ascent and reentry 


The location of these materials on 


the booster are shown in fiqure 8. 

The basic structure is a semi-monocoique design* in which 
the primary loads are carried in the external skin of the 
fuselage. The fuselage skin also forms the outer shell of the 
propellant tanks. Thus* it must withstand the stresses from 
the propellant weight as well as from internal tank 
pr essur i z a t i o n . 

The primary structure design will be aluminum skin and 
concentric stringer shells which are straight columns tied 
together by skin and frame. The tanks should be separated by 
a honeycomb construction with aluminum face sheets and thermal 
barriers. The nose will be constructed of an aluminum 
honeycomb structure capable of withstanding the high 
aerodynamic pressures of reentry. The shell has to be 
designed to withstand external pressures introduced by 
shockwaves. The propellant tank walls also have to be 
designed to resist the following loads due to supporting mass* 
drag loads* bendinq compressive loads* engine thrust and 
internal tank pressure loads. Assuming maximum shear loads 
and maximum bending moments to be at liftoff prerelease and at 
maximum q(alpha) where the dynamic pressure is greatest. The 
required wall thickness and stiffness can be established form 
these maximum loads. 


■?© 



For further studies* the winqs? fin* and canards will be 
analyzed. However the fusel aqe* which makes up the majority 
of the booster surface and heat sink weight was analyzed in 
this study. 


AERODYNAMICS AND DESIGN OF THE FLYBACK BOOSTER 


Further changes in design required the recalculation of the 
aerodynamic and control characteristics. The main changes in the 
design were an increase in overall size* the addition of one 
rocket engine* and the addition of a faring over the rocket 
engines. The aerodynamic avid control characteristics calculated 
at subsonic speeds include the location of the aerodynamic 
center* the static margin* CL* CD and CM versus alpha. At 
supersonic speeds CL and CD versus alpha at different Mach 
numbers were determined. 


DESIGN 


It is necessary for the booster vehicle to carry three 
propellant tanks* avionics equipment* mating structures* internal 
supports and other necessary components. Upon sizing the 
components and assembling them* the final booster length is 244 
feet. The diameter of the booster is 36 feet. In addition* a 


faring with a diameter of 45 
eng i nes . 

The final configuration is 
important design dimensions are: 
Overall Length 
Wing Area 
Canard Wing Area 
Tail Area (wet) 


feet was placed around the rocket 

shown in figure 9. Some of the 

= 244 feet 
= 7,161 feet2 
= 625 f eet2 
- 1926 feet2 


Cross Sectional Body Area 


1017.9 feet2 



a 


a 
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Diameter of Body 
Wing Span (main) 

Wing Span (canard) 

Air Foil 
Thickness ratio 
Sweep Angle (LE) 
Aerodynamic Center 
Aspect Ratio 
Roott Chord 
Tip Chord 
Taper Ratio 

Mean Aerodynamic Chord 


= 36 feet 
= 175 feet 
» 106 feet 
= NASA 00 IS 64 
--= 12 % 

= 45 

« 45.055% of the chord 
= 5.04 

= 85 feet 
= 8 feet 
= .094 
= 47.5 feet 


SUBSONIC AERODYNAMICS 

CL, CD, and CM versus alpha were calculated for this 
configuration. The aerodynamic center and static margin were 
also calculated. The method for finding the aerodynamic center 
is outlined in Appendix A. The aerodynamic center was calculated 
to be 45.055 percent of the wing ' s mean aerodynamic chord. The 
static margin on the pad, before staging, after staging, and at 
landing are listed in Table 6. 


Tab 1 e 6 

Static Margin in Percent Mean Aerodynamic 


Fueled For 

At Lift-off 

At Staging 

Cargo Vehicle 

1 . 025 

0 . 7324 

Shuttle II 

1 . 00.6 

0 . 7324 


Chord 

At Landing 
0.8714 
0.8758 


33 



FIG 





These values were calculated using the c,.g. values found in the 
weights section. 

The method of calculating pitching moment coefficient with 
angle of attack was determined using the Air Force's DATCOM 
manuals. The major equations used in this method are included in 
Appendix A. The values of CM for the six different 
configurations are plotted versus alpha in Figure 10. These 
plots have neqative slopes and cross the x— axes at positive 
angles of attack? therefore the configurations are stable and 
trimable. 

In order to find the subsonic lift and drag coefficients of 
the booster at various angles of attacks it was necessary to 
write a computer program (shown in Appendix B) . The equations 
used in the program are referenced from Jan Roskam's Methods.. For 

Es ti mating Drag Polars_of Subsonic_Airelanes . < See Append i x A 

for the major equations used in the subsonic aerodynamics 
program. ) 

The calculations of CL and CD were made starting at a Mach 
number of 0.7 and an altitude of 35*000 feet. This is the point 
at which the jet engines are engaged . These values were 

calculated at angles of attack ranging from -5 degrees to +19 
degrees in 2 degree increments. 

The following assumptions were necessary in writing the 
programs 

(1) Mach 0.7 was maintained from 35*000 feet (the 
approximate height were the returning booster 



reaches mach 0.7) down to 1 5 , 000 feet < the 


altitude at which the booster will cruise to its 
destination ) . 

(£) The weight change at this point in the 
trajectory is negligible. 

< 3 ) The density is based on the standard atmosphere. 

(A) The -wave drag coefficient is negligible at Mach 
0.7. 

(Additional assumptions are listed in Appendix A) 

The results obtained for CL versus alpha and CD versus CL 
are listed in Table 7. These values are graphed in Figures 11. 

The maximum lift to drag ratio was found to be 5. £11 at. an 
angle of attack of 6.0 degrees. Therefore;, the suggested angle 
of attack for the cruise back of the booster is 6.0 degrees. 

SUPERSONIC AERODYNAMICS 

In the supersonic range, CL versus alpha and CD versus alpha 
were calculated using the Air Forces' DAJC0M. The method used is 
outlined in Appendix A. CL versus alpha was found for a wing 
body— tail configuration. A proqram was written* to calculate CL. 

versus Mach number. (For the proqram listing and calculations of 
constants for program, see Appendix B.) CL as a function of 
Mach number calculated using DAJCpM is plotted in Figure 1£. 



Moment Coefficient 


FIGURE 10. 

CM vs. Alpha 


Fim led for Cargo -3 Position* 
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T ab 1 e 7 

Subsonic CD and CL vs. Alpha 


Alpha <deq) 

CD 

CL 

-5.0 

0.874 

-1 .688 

-3.0 

0 . 873 

-1 .674 

-1 .0 

0 . 878 

~1 .686 

1 .0 

0.873 

1 . 686 

3.0 

0 . 873 

1 . 687 

5.0 

0 . 875 

1 . 689 

7.0 

0 . 877 

1 .691 

9.0 

0 . 88 1 

1 . 695 

11.0 

0 . 886 

1 . 700 

13.0 

0 . 893 

1 . 707 

15.0 

0.901 

1 .715 

17.0 

0.912 

1 . 725 

19.0 

0 . 924 

1.737 
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In the calculation of the supersonic zero lift draq 
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coefficient, four different calculations of Mach number ranqinq 
from 1.5 to 3.0 in increments of 0.5 were performed while leaving 
all equations as functions of alpha. Calculations of CD were 
performed in a proqram found in Appendix B where alpha was varied 
from 10 to 60 . The data found using DATCOM is found in Table 
8. Values of CL versus Alpha are plotted in Figure 13. 

From the supersonic booster group’s proqram, the maximum CL 
was found to be at an angle of attack of 50 degrees. Therefore, 
the suqqested angle of attack for reentry is 50 degrees. 

Table. 8 



Supersonic 

CD vs. Alpha 



Mach # = 1.5 

Mach # 

- 2 

CD 

Alpha (deg) 

CD 

Alpha (deg) 

. 384 

0.0 

0 . 404 

0 . 0 

0.714 

10.0 

0 . 90 1 

1 0 . 0 

0 . 493 

20 . 0 

0 . 556 

20 . 0 

0 . 990 

30 . 0 

1 .219 

30 . 0 

0 . 708 

40 . 0 

0 . 853 

40 . 0 

1.268 

50 . 0 

0.551 

50 . 0 

0.737 

60.0 

0.882 

60 . 0 


CONCLUSIONS 

According to this analysis, the chosen booster configuration 
is feasible with today’s technology. The values that were 
calculated for the 1 moment coefficient and static margin indicate 
that the design is stable and trimable. The lift and drag 
coefficient plots have fairly standard slopes and the suggested 
angles for flight are feasible. More detailed study of this 
design would lead to further improvements in its flight 


character istics 
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FIGURE 14. 3* Views 



LIQUID ROCKET ANALYSIS 


The principal objectives in the design of the rocket 
engines for the winged flyback booster are reliability and 
simplicity. With the forementioned in mind » along with the 
thrust requirement of approximately nine million pounds and 
the cost of developing an enqine from the conceptional stage? 
the group decided to modify Rocketdyne’s F-l engine for use as 
the booster’s launch to staging engines. The basic components 
of the F~1 enqine are a tubular-wall thrust chamber? a 
direct-drive turbopump? a gas generator? and their controls. 
Figure #15? on page 45 gives a schematic representation of the 
F-l enqine. As depicted in Fiqure #15? ,the turbo pump is 
mounted directly on the thrust chamber. All other components 

t 

are either mounted on these two assemblies or are in the 
plumbing system between them. We have determined that this 
layout is perfect for our application. One of the principal 
advantages of packaging the components in this manner is that 
the high-pressure propellant ducting does not need to be 
flexed as the enqine is qimbled. Figure #15 on page 45 shows 
a simplified schematic of the F-l engine. 

An ideal rocket analysis approach was chosen to decide on 
specific modifications to be made on the ordinary engine. The 
choice of ideal rocket analysis was appropriate since it has 
become accepted practice to use ideal rocket parameters which 
are then modified by suitable correction factors. 
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FIGURE 15. F-1 Engine Schematic Drawing 



A list of assumptions for ideal rocket analysis alonq with a 
derivation of the necessary equations can be found in Appendix 
A. Furthermore* we have designed the engines -to have an 
optimum thrust at both 30*000 and 60*000 feet. This is to be 
accompanied by the use of extending nozzles. Optimum thrust 
occurs when the exit pressure of the rocket nozzle equals the 
atmospheric pressure. 

fuel 

The fuel selected for the booster is JP—A * and the 
oxidizer chosen is liquid oxygen. The reason for choosing 
this combination is two-fold. First* the ground rules 
obtained from NASA stipulated the use of liquid oxygen and 
hydrocarbon as the propellant for the first stage. Second* 
the choice of JP-4 allows both the booster engines and the air 
breathing engines to utilize the same fuel* thus saving weight 
in the form of a fuel tank and structural supports. In 
addition to the reasons stated above* JP-4* a noncryoqenic 
fuel* has a good storability which means it can be stored in 
ordinary tanks over long periods of time and at various 
temperatures without decomposition or change of state. Based 
on information obtained from Jim Sanders* M8FC (NASA) 

Engineer* the propellant has a specific impulse (Isp) of 330 
seconds at sea level. As shown in Figure #16* the vacuum 
specific impulse is a function of the percent weight of fuel 



in the propellant 



FIGURE # 16. Vacuum Isp of J P-4 and Oxygen Rocket 

From this fiqure the best oxidizer to fuel ratio is E.3 and 
the corresponding vacuum Isp is approximately 340 seconds. 

The sea level Isp is approximately 3E0 seconds;, as stated by 
Jim Sanders. The ratio of specific heat is represented by 
gamma* * and is equal to 1.24. The propellant products have 
an average molecular weight of 22.0 lbm/ lbm-mo le * and the 
propellant itself has an adiabatic flame < theoretical 
combustion) temperature of 6230° R. 

INJECTOR ANALYSIS 

The injector sprays fuel and oxygen into the thrust 
chamber in a pattern to calculate the most efficient 
combustion. The injector face as shown in Fiqure #17 on the 
combustion side* contains the injector orifice pattern* 
determined by alternating fuel and oxygen rings. 


ORIGINAL PAGE IS 
OF POOR QUALITY; 

The injector configuration can affect the local wall 
temperature in the combustion chamber and throuqout the 
rocket. The injector pattern must be controlled. Proper 
control will result in a layer of relatively cool gas near the 
wall. Poor injection will result in hot spots or possible 
burnout of the wall material. 

riATt ) 
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FIGURE # 17. Flat-Faced Injector 

The mass flow rate of propellant and the pressure drop 
across the injector are related by the following equations: 
m * Cd*/*„ *Ai (£*(<SP/H ) xym <1> 

where /* denotes the propellant density, Cd denotes the 
discharge coefficient, Ai denotes the injector nozzle area, 
and SP denotes the pressure drop across the injector nozzle. 

The discharge coefficient, Cd , can vary in value 
depending on the type injector nozzle being used. Figure #18 
on paqe 49 illustrates various types of injector nozzles and 
their respective discharge coefficients. 
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FIGURE 18. Injector Discharge Coefficients 


In general 9 injector nozzles characterized by short tubes 
with round inlets have discharge coefficients in the range Of 
0.97 to 0.99 and injector nozzle characterized by sharp 
cornered tubes have a discharge coefficient in the range of 
0.6 to 0 . 8 . 

Typical pressure drops for most injectors are of the 
order of several hundred psi . A large pressure drop promotes 
good atomization and tends to control combustion instabilities 
that are associated with pressure oscillation i*n the 
combustion chamber and in the propellant supply system. A 
large pressure drop also induces high propellant inlet 
velocities. However > it should be noted? it is desirable to 
have a pressure drop that is consistent with good atomization 
and combustion stability. 

The injector analysis performed on the modified 
Rocketdyne F-l engine required three critical assumptions. 

The first assumption was that the total fuel area and total 
oxygen area were equal to one another. The second assumption 
was that the injector orifices wore sharp edged orifices with 
a discharge coefficient of 0.61. The third assumption was to 
assume that the injector utilized a doublet impinging stream 
injection pattern with the angle of the oxygen stream equal to 
20 . Based on the basic ideal rocket analysis results, a mass 
flow rate of A685.7 lbm/sec was found along with a fuel to 
oxidizer ratio <r> of a. 3. The mass flow of oxygen <m„) and 
the mass flow rate of fuel (i»f ) was found by using the 
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followinq equations. 

inr,*<r/r+l) m <S) 

m-r»(l/r+l> m (3) 

The result of equation (S) was mo— 3267.045 lbm/sec and the 
result of equation <3) was m *.* 1420.455 lbm/sec. The same 
number of fuel orifices and oxidizer orifices were used as in 
the original F— 1 engine. They are 3700 and 2600 respectively. 
The area of one fuel orifice was calculated to be 0.1698 
inches squared and the area of one oxygen orifice was 
calculated to be 0.2147 inches squared. Referring to Figure 
noting that the areas are above 0.1 inch squared a discharge 
coefficient of 0.61 was chosen. Using equation (1) and 
rearranqiriq terms* the pressure drop across the injector was 
found to be 82.2139 psia. In order to find the anqle of the 
fuel flow* the followinq equation was used. 

rn-r *Vf*sin(i>) = m CT *Vo*sin< t„ ) (4) 

In order to determine the fuel flow anqle > an arbitrary anqle 
of 20° was selected for the oxygen flow and an impingement 
anqle of 0® (axial flow) was selected at the point in which 
the fuel and oxyqen intersect. The velocity of the fuel and 
oxidizer was determined by the followinq equations: 

Vo * Cp* ( 2*qc*<S'P/ ) 4 ym <5) 

Vf = Cd* < 2*qc*<5'P/ J*t ) 1 <6) 

The velocity of oxidizer was found to be 10.516 ft /sec and the 
velocity of fuel was found to be 12.903 ft/sec. The fuel 
anqle was found to be 69.595 using equation (4). 
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GAS GENERATOR 

The power necessary to drive the turbine is obtained from 
the qas qenerator » which qenerates qases by a chemical 
reaction of propellants similar to those iri the thrust 
chamber. The qenerated qases have to be cooler than thrust 
chamber reaction gases* since excessive temperatures will 
cause failure in either the turbine nozzles or the turbine 
wheel. The gas generator burns a fuel-rich mixture of the 
same propellants used in the thrust chamber. It burns 
approximately two percent of the total propellants used in the 
engine. The qas generator used for the modified F— 1 enqine 
will be the same one used in the Rocketdyne F-l. It is 
partially spherical in shape and is approximately ten inches 
in diameter. The basic desiqn of the F-l qas generator 
incorporates a doubled wall combustion chamber* throuqh which 
fuel flows to reqenerat i ve.ly cool the body. 


TANK PRESSURIZATION AND PUMP ANALYSIS 

The main objectives of tank pressurization is to keep the 
components so small and as lightweight as possible. On 
account of this* the pressures in the tank must be much lower 
than pressures required at the injectors. To perform the task 
of developing these large pressures* turbo-centrifugal pumps 


are used to transport the propellants from their tanks to the 


injector. These pumps , one for the fuel and one for the 
oxidizer , provide a dynamic pump head <TDH) which is 
proportional to the pressure difference between the injectors 
and propellant tanks. 

The propellants are supplied to their respective pumps 
throuqh a sinqle inlet * in line with the main shaft* and is 
discharged radially through dual outlets. The dual outlet 
design provides a balance of the centrifugal loads in the pump 
which, in turn, minimizes the pump diameter, see Figure #19 . 

The pumps are powered by a turbine. The turbine is 
driven by hot fuel gases that are heated in a small gas 
generator, discussed in the for ementi oned section. The 
turbine gives approximately 55,000 brake horsepower to the 
oxidizer turbopump and about 25,000 brake horsepower to the 
fuel pump. 

The turbine exhaust gases, which are gaseous fuel, are 

collected by a manifold and ducted through a nozzle skirt 

into the engine exhaust plume. Using a circumferential 

manifold the turbine exhaust qases introduce a cooler boundary 

* 

layer to protect the wall of the nozzle. This idea allows for 
a comparatively liqht weiqht system with no need for turbine 
exhaust ducts and other attachments. 




where TDH is the dynamic pumphead, <£P is the pressure change 
from tank to injector, and Q is the volumetric flow rate. For 
the oxidizer liquid oxygen, the density is 71.2 sluqs/ft $ 
therefore, the specific gravity equals 1.141. Assuming a pump 
efficiency of 80%, <5'P = 2664.4 psia. 

P«t«r>*« SKldliar = P»nJ*B«ior - 1065.5 psia 

For the fuel JP-4, the density is 47.3*sluqsi/f t 5 therefore, 
the specific gravity equals 0.758. Assuming the same 
efficiency, SP ■ 2873.3 psia. 

Pl—r,* = Pvn.t-C*nr “ SP * 855.4 psia 

THRUST CHAMBER AND HEA T T RANSFER 

In the drive to produce large, high-pressure engines, a 
major problem was a satisfactory means to cool the thrust 
chamber. Therefore, this section will discuss the 
configuration of the thrust chamber, the strength of the 
thrust chamber material, and how the heat transfer of our 
rocket was determined. 

The thrust-chamber assembly consists of a tubulai — wall, 
reqeneratively cooled chamber with an uncooled extension, a 
double— inlet oxidizer dome, four integral fuel valves, and a 
flat-faced injector. The cooled portion of the thrust chamber 
has a 40 inch combustion chamber diameter , an 8 foot nozzle 
exit diameter, and is 11 feet in lenqth. The chamber is 
designed for an uncooled extension which will be deployed at 


45,000 feet 


Combustion chamber temperatures of rocket propellants 
typically are higher than the melting points of most common 
metals and alloys. However, at high temperatures the strength 
of most of these materials decline rapidly. This can be seen 
in the figure below. 



FIGURE # 20 Variation of Tensile Strength with Temperature 

Since the wall thickness of the thrust chamber depends 
strongly on the stresses which it can support, it is desirable 
to use highly stressed materials. The wall should be cooled 
to a temperature considerably below its melting point, and 
below the propellant stagnation temperature. 

After researching several different materials, the rocket 
propulsion group decided to use stainless steel as the 
material for the combustion chamber. Two major reasons for 
this choice are economically feasible and it is readily 



available 


Since a liquid rocket is being used* it tuts already been 
stated that it will be reqeneratively cooled. This means that 
the fuel or oxidizer will flow in tubular passaqes directly 
outside the chamber wall. These tubes reduce wall thickness 
and thermal resistance and* more importantly* increase the 
coolant velocity in the throat section. 

For our liquid rocket* we selected liquid hydroqen as the 
coolant. After careful research of the thrust chamber wall 
material* stainless steel, it was decided to use a wall 
thickness of 0.125 inches. 

In order to determine the heat transfer, it is first 
possible to find the total heat absorbed per second in the 
thrust chamber. The calculated value for this is: 
q - 8164.8 BTU/sec <7> 

Then the average thrust chamber area can be found by 
interqrating from the combusion chamber area to the exit area 
of the cooling section of the nozzle. The calculated average 
area is: 

A * 4536.0 in 2 <8) 

* 

Now, by dividing equation (7) by equation <EI) it is possible 
to find the average heat transfer of the combusion chamber. 

The calculated heat transfer is: 

q/A * 1.8 BTU/ in 2 -sec <9) 

From the calculated average value of the thrust chamber* it is 
easy to see that the heat transfer rate is high. 


COMBUST I ON I NSTAB I L I T. I ES 

Rocket engines are sometimes subject to combustion 
instabilities in the form of large pressure oscillations 
within the chamber. Such instabilities can cause engine 
failure either through excess pressure, increased wall heat 
transfer » or a combination of the two. Because combusion 
instabilities cannot often be observed on smal ler— scale 
models* the full scale model would have to be tested to detect 
any instabilities. Moreover* it is impossible to determine 
instabilities analitically and experimental test must be 
performed to insure that no instabilities exist* therefore* we 
cannot evaluate this area of the design at this stage of 
conceptual design. 

NO ZZLE AN ALYS I S 

The design of a nozzle requires taking into account 
variations in velocity and pressure on surfaces normal to the 
streamlines. Improper shaping of the nozzle can result in 
shock formation and substantial performance los3. It has been 
found empirically that simple conical divergent nozzles 
provide best performance when their half-angles are betwwen 12 
and 18 degrees. <10* 193) Furthermore* the bell shaped 
nozzle permits additional advantages in reducing size and 
weight when compared with the standard lb degree half-angle 
conical nozzle. Without any reduction in performance, the 


bell shape also permits a 20 percent reduction in length. 

< 10s 194 ) 

The F-l engine currently features a bell-shaped nozzle. 
The proposed modified F-l engine will feature a conical 
nozzle to satisfy the mission requirements. Through the use 
of the program in Appendix A, the area expansion ratios were 
found to be 55 and for altitudes of 30 >000 and 60 >000 feet 
respectively. As shown by Figure #21, on paqe , the thrust 
is to be optimized at 30,000 and 60,000 feet. The.ncrzzles 
will be optimized for an altitude of 30,000 feet at the 
launch. Once the vehicle has reached approximately an 
altitude of 45,000 feet, the nozzle extenssion will move down 
to an area expansion ratio which qives optimum conditions at 
60,000 feet. The nozzle will remain at this position until 
staging. Once the liquid rocket engines are shut down and 
staging is complete, the nozzles will be retracted to their 

initial area expansion ratio. In doing this the drag will be 
slightly reduced. 

THRUST VECTORING 

The thrust-vector control is achieved by gimballing the 
entire engine. The hiqh-pressure fuel is used as the 
hydraulic actuating medium. Although this is an 
unconventional approach, it was the method used on the Saburn 
V’s first stage. The use of the fuel has its drawbacks; 
however, it eliminates a separate hydraulic system. 





RESULTS 


Parameters at Optimum Altitude 


30,000 feet 
Pa-4.36 psi 
P«,«-3731.17 psi 
T %■— 6a30°R 
T_-1800«R 
T*«5563°R 
M..-4.53 
A*-8.01 ft2 
A«,*86 .59 ft* 
€-43 


60,000 feet 

Pa-1.04 psi 

Pi»«— 3731.17 psi 

T*„-6230°R 

T„=1572°R 

T*=5563°R 

M.,— 4.97 

A*— 2. 01 ft* 

A*-" 145. 37 ft 2 
€-78 


Nozzle Dimension 
Ha If -Angle 

Length of Cooled Nozzle 
Lenqth of Unextended Nozzle 
Length of Extended Nozzle 


i6.aa° 

11.00 ft 
15.30 ft 
20.63 ft 


CONCLUSION 

The basic parameters of the modified Rocketdyne F-l 
rocket engine have been determined. These parameters include 
the followings the temperature and pressure of the combusion 
chamber, the temperature, pressure and area of the throat, and 
the exit area of the rocket nozzle. The rocket features an 
expandable nozzle that is optimized at two different heiqhts. 
These heights are 30,000 feet and 60,000 feet. The analysis 
of the engine also includes and injector analysis, turbopump 
analysis, and a heat transfer analysis based on data obtained 
from trajectory analysis. The decision has been made to use 6 
F-l engines to satisfy the booster’s thrust requirement. JThe 
Rocketdyne F-l engine will be used; however, there will be a 
few modifications made to the nozzle, injectors, and turbopump 
based on our data. The F-l was chosen since it is economically 
feasible and readily available. A second reason the F— 1 was 
chosen was because the research and development cost for a new 
rocket engine are astronomical and the development time is in 
excess of five or more years.’ 

It is the recommendation of this group that Rocketdyne be 
awarded the contract to redesiqn the F-l engine to meet the 
requirements specified by both the cargo and shuttle missions. 
It is estimated that the redesigning and testing of the F— 1 
engine will take approximately two years. The proposed cost 
of one of these engines in 1986 dollars is approximately 6A.2 
million dollars. It is also recommended that the modified 



engines be tested at either Edward’s Air Force Base or 
Rocketdyne’s Santa Susana facility. Both facilities are 
capable of handling 1.5 million pound-force thrust levels. 


AIR BREATHING JET PROPULSION 


After the flyback booster reaches 7,000 feet per second, 
the rocket engines of the booster will shut down. At this - 
time the manned Shuttle II or the unmanned cargo vehicle will 
separate from the flyback booster. The flyback booster will 
then return to a landing destination powered by turbofan 
engines. While the flyback booster is in this boosting stage, 
it will reach velocities that would destroy the turbofan 
engines. Therefore, it is necessary to protect the turbofans 
from damage with the use of inlet/outlet enclosures. These 
protecting enclosures, during reentry, will open when the 
booster's velocity has decelerated to subsonic speeds. At 
this point, the turbofans will fire up arid power the booster 
to a predetermined landing destination. 

ANALYSIS!. The main information which will be found in this 
section of the report is: <1) the total thrust of, the weight 
of, and weight of fuel used by the engines, and (2) the design 
and placement of the specialized cowlings needed on the 
engines. To determine the this information. General Electric 
and NASA were consulted. See references 12 and 13. 

The first step in this analysis was to determine which 
type of jet engine should be used. Since one of the 
specifications for this vehicle was to fly at low velocities 


(Mach Numbers less than 1.0) during powered flight, a turbofan 
engine would be best suited for this particular case. This is 
due to the fact that turbofans have a high thrust to weight 
ratio* and at low specific fuel consumption at low velocities. 

The next step in the analysis was to determine which 
bypass ratio was best suited for this case. NASA suggested a 
medium bypass ratio of approximately two. General Electric 
has in operation two engines with a bypass ratio of two, with 
one augmented with an afterburner for hiqh thrust output. For 
comparison, further analysis was performed on three types of 
General Electric engines, with two different bypass ratios: 

(1) CF6-8032 <5 bypass), (2) F101 (2 bypass), (3) 
FlOl-augmented (2 bypass) . 

In determining the type of engine, the number of engines, 
and the amount of fuel use by the engines best suited for the 
flyback application, a program, obtained from NASA, was put 
into operation (program 2). This program considered a 
particular flyback vehicle from staging position and 
determined the flying characteristics of th€? return 
trajectory. To make the program fit this particular booster, 
engine performance and aerodynamic data were* needed to be 
inputted into the program. The engine performance consisted 
of thrust, and specific fuel consumption, at different Mach 
Numbers and at different altitudes. Also, the program had a 
built in engine-out arid ten minute qo-around the landing field 
capability. NASA’s program had the engines ignite at an 





altitude of 35,000 feet. The vehicle would then begin a high 
cruise at 15,000 feet and let down to a low cruise of 10,000 
feet before landing. The engines would run at a maximum time 
of 3,175 seconds (53.4 minutes). 

In determining the performance of each engine in the 
analysis, information again obtained form NASA and General 
Electric was used. General Electric furnished the dimensions, 
sea level static thrust, and specific fuel consumption for 
each of the engines. NASA furnished the equations and data 
necessary to change the sea level thrust and specific fuel 
consumption for different altitudes and Mach Numbers. The 
following equation was used to convert the sea level thrust to 
different altitudes: 

F Z f f ZS . *?.££T 

Aut SSL N ~T^ux 9 S5£- 

To find the thrust at different Mach Numbers, data from NASA's 
generic engines was used (reference 13). The following ratio 
equation was applied: 



The above relation is an approximation which does have some 
error, yet, according to NASA, the error should be small. The 
specific fuel consumption also used a ratio assumption to 
determine its change due to Mach Number and altitude. 



The next step in the analysis was to determine where the 
enqines would be installed on the booster. Because of 
aerodynamic stability, the enqines will be installed over the 
canards towards the front of the booster. Because, when the 
booster is powered flight, it is a flying empty fuel tank with 
500,000 pounds of dead rocket engines on the tail end. The 
weight of the jet will help bring the center of gravity 
towards the nose of the booster resulting in a more stable 
flight . 

One of the key points in the i vista 1 lat ion of these 
engines is the desiqn of the specialized cowlings. The 
cowlings must accomplish two different tasks. It should 
enclose bhe engines completely during the high Mach Number 
boosting stage, yet it must be made streamline to reduce 
aerodynamic drag. Figure 22 illustrates one possible design 
that this report recommends. The engines will resemble this 
figure until an altitude of 35,000 feet. At this point, the 
welded seams on the blow-off caps will split. Then the cone 
shaped enclosures on the inlets and exits of the cowling will 
divide up into four different slices. Immediately after, each 
piece of the cone will retract quickly inside the inner sleeve 
in the cowling (Figure 23). The enclosure cones and the 
cowling should be made of the same material as the booster’s 
outside skin, and sprayed with a carbon-carbon composite; 
since, the curvature of the cone is relatively blunt. This 
should protect the cones from the high temperatures during the 
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boosting staqe. 

The most important part of this section’s analysis is to 
decide on which engine to use on the booster. To make this 
decision? three important aspects should be viewed. The first 
aspect is the total weight of the engines and fuel. The 
second aspect is the aerodynamic drag. Since the aerodynamic- 
drag is directly proportional to the frontal area? the 
following equations were used: 

Drag*-~ Frontal Area or D <*+ S 

■ 2 - 

s - -rrd/4 

Drag » <0.50 C^> V 1 ) ird/4 

The last aspect is the number of engines. The more 
engines used? the more money? structures? electronics? pumps? 
feed lines? and mechanical parts will be needed to supply the 
system. So? the above three factors will have to be weighed 
against each other on deciding on the type of engine. 

In addition? the engines will be running on the same fuel 

(JP-4) as the main rocket engines. The JP-4 will be pumped 

* 

•from the booster’s main tank. When installing the engines* 
the nozzles should be angled six degrees downward from the 
longitudinal axis for maximum performance; since? the booster 
will cruise during flyback at an angle of attack of six 
degrees. 
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RJ^ULTSs After reading in the 
and specific fuel consumptions 
following data was outputed: 


Eng i ne 

Engine 

Weight 


Weight 

of Fuel 


( lbs) 

( lbs) 

CF6- 

38805 

108587 

80C8 

- 


F 1 0 1 
w/ AB 

31680 

581138 

F101 
w/o AB 

46338 

130449 


Table 9: Turbofan 


three different engines' thrust 
into NASA’s program , the 


Total 
Weight 
( lbs ) 

Frontal 
Area 
( ft-- a ) 

Number 
of Engines 

134798 

873 . 68 

4 

558818 

191.45 

8 

176781 

311.08 

13 


Engine Output Data 


From these results, the CF6-80CS is the best engine to use for 
the flyback booster. The CF6-80CS is better than the F101 
(without afterburner) in all categories, yet the augmented 
F101 seems to give the CF6-80C8 a challenge. The F101 with 
afterburner is lighter in engine weight by 585 lbs., and also 
has a smaller frontal area by 88 square feet. However, the 
CF6-80CS uses 80 V* less fuel, (418,545 lbs. less) which makes 
this the deciding factor in the determination of the engine 
type. After taking all the categories into account, the 
CF6-80CS is the better selection. 

Below is the data on the CF6 -80C8 furnished by General 
Electric (Figure 84). 
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Bypass Ratio 
Fan Tip Diameter 


5.09 


;>> 


> 


Length 

Weight 

Thrust 

Sp <sc i f i c Fue 1 Cnsp t . 

Thrust /Weight 

Fuel 

Max. Cowling Dia. 
Frontal Area 

( 10% change due to techn 


93 in. 

160.9 in. 

** 8,946 lbs. ** 
*•* 65,500 lbs. ** 
0 . 40 
7.763 
JP-4 
112 in. 

9 , 858 i n . •"'£ 
logical advancements ) 


3 

Some of the data above was changed from the 1986 data due 
to future technological advancements. These advancements were 
suggested through Ms. Mary Pryor at General Electric. See 

■> 

Appendix C for calculations. 

The booster will require four CF6-80CS turbofan engines. 
If all engines fire up, each engine will run at a maximum of 
85% power. If one engine fails, the remaining three will run 
at over 95% power for 30 minutes. After 30 minutes, the 
engines will throttle down to under 90% power for the rest of 
the flight. 

In table 10, a summary of the thrust performance is 
presented, this data was inputted along with the specific fuel 
consumption data, in table 11, into NASA's program. The 
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calculations for both tables is presented in Appendix C 


ALTITUDE S.F.C. 

<feet) < lbm/lbf -hr ) 


0 

1524 

3084 

4572 

6096 

1066a 


.3321 
. 3985 
.4151 
.4317 
.4317 
.4317 


Table 11: Specific Fuel Consumptions at Different Altitude 



CON CLU SION: With weighty cost* drag* electronics* and other 
data taking into account* the CF6-80CS is the obvious choice 
of the engines analyzed. 

In the designing of the enclosure cowlings* testing will 
have to be an integral part in the design. This is mainly due 
to the high complexity of the opening device. Also* the four 
piece cone may under go flutter as the blow-off cap releases 
its tension. The design of the enclosure cowling may also 
have some drawbacks in weight, yet they are completely 
reusable. Ori the other hand* if the cones were blown off 
completely, they might risk damaging the wing, fuselage, or 
vertical tail on their trajectory towards the ground. Because 
of the above reasons* enclosure cones which retract or just 
open would be prime answer to the problem of enclosing 
turbofan engines during supersonic flight. 



WEIGHTS 


When a major project like the flyback booster is 
attempted, a beginning point is needed. In most cases the 
starting point for a technical project begins with weights. 
Weights set the parameter each project engineer must work 
under and also are used when estimating the cost of the entire 
project . 

For the flyback booster, a weight specialist has two main 
objectives. First, a complete breakdown of component weights 
is needed. Weights are assigned to various components of the 
booster and that information is passed to the various 
engineers working on the project. Once a component is 
assigned a certain weight, the engineer must complete his 
project under that weight. If the engineer finds it 
impossible to stay under the assigned weight, the weight 
specialist must recalculate all numbers and assign new 
weights. The second objective is to calculate the cutoff 
weight. This is essential for the performance engineer. 

A computer program was obtained from MSFC called Ascent / 
Vehicle Interactive Design (AVID) and was used to accomplish 
both objectives. AVID is preliminary design software used in 
establishing performance capabilities over a wide range of 
aerodynamic designs. 

Before AVID can be used, a preliminary vehicle must be 
-constructed, and its dimensions must be fed into AVID. 
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Information on engine type* propellants used? material 

constants? arid modes of operation must also be fed to AVID. 

After this information has been entered? the different weights 

of the booster can be calculated. 

To calculate the different weights certain formulas are 

implemented. Some of the formulas are shown below. 

The fuselage weight is calculated by: 

W ( fuse ) »A < 5 ) *SB ( wet ) 1 - ( RS TR ) 

where A<5) is the fuselage weight constant 
SB (wet) is the body wetted area 

The thrust structure weight is calculated by: 

W< thrust ) «A ( 6 ) * < TH 1 *enq +TH2*enq > 1 * oaa *< RSTR ) 

where A(6) is the thrust weight constant 

TH1 is the mode 1 engine vacuum thrust 

eng is the number of engines 

THE is the mode 2 engine vacuum thrust 

The landing gear* weight is calculated by: 

W< lng)»A< 12>*W< lnd >«< As * » *< RSTR ) 

where A(12) and A(13) are landing gear 
constants 


The reaction control system weight is calculated by: 

W(rcs)=A( 16>*W(ent)*(RSS) 

where A(16) is the res weight constant 

W(ent) is the vehicle entry weight 
RSS is the subsystem weight reduction. 


@ 
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AVID can calculate the weights for thirty components 
of the flyback booster. See Table 12 for a list of the 
components and their weights. 


nr 


Center Of Gravity Calculations 


£ 
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As in all flight vehicles it is imperative that the 
vehicle's center of qravity (c.q.) location be within a 
certain range along the vehicle’s longitudinal axis. On the 
flyback booster this ranqe must be located such that the 
vehicle is controllable at all times during the trajectory. 

In order to keep the booster’s c.g. within this ranqe, the 
location of each item is carefully considered. 

The following method is used to find the c. -q. along the 
longitudinal axis. 

1. Determine a datum or reference from which the c.g. of each 
component will be measured. The tip of the booster’s 
nosecone is used for the datum. 

S. List each component and its weight. 

3. Measure each component’s c.q. location from the datum to 
obtain the lever arm. 

4. Multiply each component’s weight by its lever arm length 
to obbain the moment arm. 

5. Sum up all the weights to obtain total weights and 
moments to obtain total moment. 

6. Divide the total moment by the total weiqht to obtain 
the c.q. location as measured from the datum. 

These calculations are shown in the following table for 
the booster with all propellants removed. The location of the 
c.g. is shown in Figure 35 on page 81. 
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COMPGNEN r 

WEIGHT 

DISTANCE 

MOMENT 


(LBS) 

( FEET ) 

(FOOT-LBS/ 

No se S t r uc t ur e 




with Avionics 

8131 .0 

20.0 

162.2 

Canard 

4000 . 0 

74 . 0 

296.0 

Structure 1 

95700 . 0 

29.0 

2775.3 

Struc ture S 

820000 . 0 

74.0 

16280.0 

Structure 3 

156000.0 

149.0 

23244.0 

Structure A 

190000.0 

206 . 0 

39140.0 

Nose Wheel 

1513.0 

74.0 

112.0 

Wings 

84298 . 0 

220 . 0 

5344 . 0 

Tail 

3378.0 

229 . 0 

773 . 6 

Thrust Structure & 




Therm. Protect. 

18488.0 

1 44 . 0 

£663 w 3 

LH2 Tank 

12010.0 

48. S 

582.5 

LOE Tank 

28764 . 0 

1 14.0 

3279.0 

JP4 Tank 

17880.0 

175.0 

31 £9.0 

Ha in Gear 

7533 . 0 

204 . 0 

1536.7 

React. Cent. Sys. 

2000 . 0 

80.0 

40 . 0 

Air Breath. Eng. 

32064 . 0 

74.0 

2372 . 7 

Rocket Engines 

121200.0 

232 . 8 

28215.4 

Eng i ne G i mb ling 

2 1 00 . 0 

214.0 

449 . 4 

Vehicle Mounts 




For e 

7000 . 0 

o 

• 

518.0 

Aft 

7000 . 0 

£06 • 0 

1442.0 

Eng i ne P 1 umb i ng 

8344 . 0 

175.0 

1 460 . 2 

Hydr au lie Sys t . & 




Control Surfaces 

5547.0 

175.0 

970 . 7 

Rocket Eng ine 




Ferr inq 

3000 . 0 

224 . 0 

672.0 


TOTAL 

WEIGHT = 975 

944 lbs 

CENTER 

OF GRAVITY LOCATION = 138 

.8 feet 


TABLE IE 


List of Structural Components and Weight 




a 
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In addition the c.q. was calculated at lift-off? 
separation? and at landing. These calculations were necessary 
in order to see how the propellant weights would affect the 
c.g. at various times during the trajectory. Since the 
booster is used for two different missions? to boost the 
shuttle vehicle and to boost the cargo vehicle? it is 
necessary to find the c.q.’s at various times in the two 
missions. These c.g. calculations are shown on the following 
table. 


CARGO VEHICLE 


SHUTTLE VEHICLE 


AT LIFT-OFF. 


weight = 6227150 lbs 
c”.q. = 131.69 feet 


weight = 6224052 lbs 
c.g. — 132.61 f ee t 


BEFORE/ AFTER SEPARATION 


weight — 1112910 lbs 
c.g. * 145.61 feet 


weight ~ 1112810 lbs 
c.q. = 145.61 feet 


AT LANDING. 


weight = 991544 lbs 

c.g. — 139.01 feet 


weight = 975944 lb 

c.q. = 138.80 fee 


TABLE 13. Weights and C.G. At Lift-off? 
Separation? and Landing 
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COST ESTIMATION ANALYSIS 


This cost model contains normalised Cost Estimating 
Relationships <CERs> for the Flyback Booster Vehicle System. 
The cost model was developed for the Engineering Cost Group 
( PP03 ) of NASA's Marshall Space Flight Center (MSFC) by PRC 
system services, a Planning Research Corporation Company:, 
located in Huntsville* Alabama* and is applied to the Flyback 
Booster/Shuttle Orb iter < FBSG ) project of Auburn University's 
Aerospace Engineering* AEVi-9 FBSO group. 

This model provides the FBSO group with the CERs required 
to prepare the estimates of the Design* Development* Test and 
Evaluation < DDT & E) and Fliqht Hardware <FH) manufacturing 
costs of the Flyback Booster Vehicle System (FBVS). To this 
end CERs have been provided using data from historical 
operations of previous Launch Vehicle > Aircraft* Spacecraft* 
Boosters* Liquid Rocket Motors and Propellant Tank data. These 
data were used to augment launch vehicle data in several of 
the subsystem CERs. All data used in this model were extracted 
from MSEC'S REDSTAR (Resource Data And Retrieval > 4 data base 
system in Huntsville* Alabama. 

The recommended CERs contained under each individual 
subsystem area are the result of regression combinations of X 
and Y values for the entire data base. To formulate the 
recommended CERs the data recorded consisted of two 
parameters* cost and weight* which were graphically presented 



on log-log format. From these plots the CERs are determined as 
a result of curve fitting and determining the mathematical 
equation for each curve with supporting statistical data. It 
in the case due to the limited data available* slopes of 0.5 
for D & D and 0.7 for FH & A were assumed. 

With all the required CERs determined* the Launch Vehicle 
Cost Model <LVCM) was generated and is shown in Figure 86. The 
L.VCM is arranged in spreadsheet form using the program "LOTUS 
183". At this point it is necessary to analyze the basic 
structure of the LVCM* in order to gain insight on the 
fundamental operation. 

The LVCM structure allows the user to develop cost 
estimates at the subsystem and system— levels* separated 
between Non-recurring Design* Development* Test and Evaluation 
(DDT &< E) and Recurring Flight Hardware (FH) unit cost. It 
serves to categorize those hardware components* subassemblies* 
and assemblies into generic hardware groups or subsystems. Tire 
primary function of the cost model structure is to establish a 
functional commonality of both system and subsystem-level 
elements. 

Application of the model is based on the bui Iding— block 
approach. A sequence of calculations* utilizing the CERs* are 
made by the user for each subsystem and system-level element. 
The known performance characteristic (independent variable or 
CER parameter) is introduced as the input into the worksheet. 
All independent variables (CER input) are subsystem dry weight 
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STAGE ONE 
COST ESTIMATE 







DATE: 

22-Nar-87 






INFLATION : 

1.229 

FROM 82$ TO 86$ 




UT RESERVE: 

1.23 






CER'S 


COMPLEXITY 

COST 



1NDEP. VAR. 



FACTOR 

(86 $ MILLIONS) 

COST ELEMENT 

DDTfcE 

TFU 

DDT&E 

TFU 

DDT4E 

TFU 

STRUCTURES/TPS 
THERMAL CONTROL 

734354.00 

734354.00 

1.00 

1.00 

1112.42 

1104.14 

EJWIROItfNTAL CONTROL 3358.00 

3358.00 

1.00 

1.00 

17.32 

1.68 

BASE HEAT SHIELD 

18488.00 

18488.00 

1.00 

1.00 

44.15 

12.10 

UING/TAIL/LEAD1NG EDGE 

31670.00 

31670.00 

1.00 

1.00 

478.30 

86.92 

LANDING 6EAR 

9046.00 

9046.00 

1.00 

1.00 

63.68 

15.95 

AVIONICS 

8131.00 

8131.00 

1.00 

1.00 

1171.69 

119.98 

ELECTRICAL POWER 

10639.00 

10639.00 

1.00 

1.00 

439.15 

113.02 

PROPULSION (LESS ENGINES) 

13444.00 

13444.00 

1.00 

1.00 

262.40 

34.13 

SEPARATION PROVISIONS • 

2134.00 

2134.00 

1.00 

1.00 

41.51 

16.49 

SURFACE FLIGHT CONTROLS 

2129.00 

2129.00 

1.00 

1.00 

185.33 

50.78 

AUXILARY POWER UNIT 

6387.00 

6387.00 

1.00 

1.00 

252.51 

81.36 

HYDRAULICS 

5547.00 

5547.00 

1.00 

1.00 

112.64 

28.06 

SUBTOTAL 

815686.00 

815686.00 



3326.70 

1453.78 

STRUC. T00LIN6 

1104.14 


1.00 


1631.19 


SYS. TST. HROU. & ASSEMBLY 1453.78 


1.00 


1817.65 


SYS. TST. OPS. 

1817.65 


1.00 


520.45 


SUBTOTAL 





7295.99 

1453.78 

6SE 

7295.99 


1.00 


875.05 


SUBTOTAL 





8171.04 

1453.78 

SE&I 

8171.04 

1453.78 

1.00 

1.00 

663.66 

94.64 

SUBTOTAL 





8834.70 

1548.42 

PR06RAM MANAGEMENT 

8834.70 

1548.42 

1.00 

1. 00 

233.96 * 

41.12 

SUBTOTAL 





9068.66 

1589.53 

ENGINES (OHwTi»«ISP«Pc) 

21706454000.00 130238720000.00 

1.00 

1.00 

91.21 

0.21 CE 

SUBTOTAL 





9159.87 

1589.74 

FEE (147.) 





1282.38 

222.56 

PROG. SUPPORT (37. DEV., 27. PROD.) 




313.27 

36.25 

SUBTOTAL 





10755.52 

1848.55 

COST CONTINGENCY (157.) 





1613.33 

277.28 

TOTAL 





12368.84 

2125.84 


FIGURE 26. 

Launch Vehicle 

Cost 

Model 






with the exception of rocket engines which require a composite 
variable of engine dry weight X vacuum thrust X Isp X chamber 
pressure. The subsystem dry weights are the result of weight 
estimation performed on the FBSO* which are also documented in 
this report. The resulting output of each subsystem element is 
performed and a cost estimation of DO & E and Theoretical 
First Unit <TFU) cost is calculated. This is; presented in the 
far right hand column of the LVCM worksheet. The following 
documentations are a simplified break down of all the 
subsystems used in determination of the cost estimation of the 
FBVS. Each subsystem includes all the historical data utilized 
to determine the recommended CERs and includes a basic break 
down and description of what each subsystem consists of: 

(1) STRUCTURES /THERMAL PROTECTION SYSTEM: 

All CERs are the result of historical structures data 
from the S— IC* S-II* S-IVB and Centaur boosters* Shuttle 
Orbiter* External Tank and IUS. All visible tooling cost have 
been removed from the structures/TPS cost to give a more 
accurate Design and Development Cost <D & DC). 

This element includes the forward* mid* aft fuselage and 
crew module* frame structure* inter tank* thrust structure* 
tank support structure and combined LOX and fuel tanks. 

Thermal control related to these elements is included in the 



(2) THERMAL CONTROL 


The CERs are the result of historical data from S--IC* 

S-II and SO. This subsystem covers those thermal control items 
not in the other subsystem. Due to limited data available* all 
slopes have been assumed to be 0.5 for D?xD and 0.7 for FH ?x A. 

This subsystem is broken into two parts: 

Base Heat Shield: This element protects the thrust 
structure* propellant tanks and other elements from excessive 
heat and gases given off by the engines during their burn. 

Environmental.. Control : This element provides temperature* 
humidity* and hazardous qases control to the forward and aft 
mounted equipment containers. 

(3) WING/TAIL/LEAD I MG EDGE: 

These CLRs are from historical data from the C-5A* C-1A1* 
KC-.1.35* C-130A aircraft and Shuttle Qrbiter. 

this element provides the FRVS with aerodynamic .lift and 
stabilization functions. Major elements included are wings* 
tail and leading edqes. This system consists of the wing box 
structure* leading and trailing edge structure and leading and 
trailing edge control surfaces. 

<.'*) LANDING GEAR: 

The CERs are from historical data from the C-5A* C-1A1 » 
KC-135* C-130A aircraft and Shuttle Qrbiter. 

This subsystem provides the launch vehicle with the 



capability to land. Major elements included are; 

* Landing Gear 

* Supports 

* Landing Gear Hydraulics 
(5) AVIONICS: 

Historical data from numerous space programs were used to 
develop these CERs. 

This subsystem serves as the control element of the FBVS. 
Its functions include guidance* navigation and control; 
display and control of vehicle functions and status; sensing 
to monitor vehicle operations? and communications and data 
handling to assure proper operation of the vehicle subsystem 
during its mission. The Avionics Subsystem is composed of four 
main elements; 

* Guidance* Navigation* and Control 

* D i sp 1 ay and Con tr o 1 s 

* Instrumentation 

* Communications and Data Handling 
< 6 ) ELECTRICAL POWER: 

Historical data from S-IC, 8-1 1, S-IVB, Centaur— D, IUS 
and Shuttle Orb iter were used to develop these CERs. 

This element includes the power source; electrical power 
air-conditioning* distribution* and control equipment. The 
subsystem provides the power* from batteries or fuel cells 



through power distribution system? to operate the 
pressur ization? propellant management ? engine controls) 
staging and separation systems) controls for hydraulics 
systems? telemetry? display and controls? data handling? and 
communications systems. 

< 7) PROPULSION SUBSYSTEM (LESS ENGINES): 

These CERs are develop from historical data from 
Centaur-D? S~IC? S— II? S—IVB? External Tank and Shuttle 
Orb iter . 

This subsystem provides the launch vehicle rocket engine 
with fuel. Major elements included are? 

* Propellant Feed System 

* Recirculation System 

* Propellant Management System 

* Pressurization Sy s t em 

* Plumbing? Valves and Lines 

* Fuel Fill? Drain and Vent System 

(8) SEPARATION SUBSYSTEM: 

Historical data from the Shuttle Or biter -ET Separation 
System? S-II Separation System and the Shuttle SRB were used 
to develop these CERs. 

The subsystem provides a means of separating the shuttle 
orb iter or from the FBVS as their requirement for the mission 


<9> SURFACE FLIGHT CONTROLS! 


Historical data from the Shuttle Orb iter surface flight 
controls and actuators were used to develop these CERs. 

This subsystem provides the flight control surfaces and 
the mechanisms required for proper control of the FBVS during 
atmospheric flight. 

(10) AUXILARY POWER UNIT; 

Historical data from the Shuttle Orb iter auxiliary power 
unit were used to develop these CERs. 

This subsystem provides the power to drive the hydraulic 
pumps which provides the power for the TVC and other hydraulic 
power mechanisms. 

<11) HYDRAULIC SYSTEM! 

The required CERs are developed from historical data from 
the Shuttle Grbiter hydraulic system. Due to the limited data 
available assumed slopes of 0.5 for D & D and 0.7 for FH & A 
were performed . 

This subsystem consists of components required for 
generation? control distribution? monitoring? and use of 
hydraulics. The hydraulic power operates the aerosurf ace 
controls (aileron? eleron? rudder/speed brake? and body 
flaps)? the FBVS main engine control valves? lockup and unlock 
of landing gear; operate main wheel brakes and provide nose 



wheel steering. 


At this point it is necessary to recognize that all 
estimated cost of each subsystem are based upon certain 
assumptions and ground rules established prior to the initial 
study. This is mainly attributed to the amount of progress 
obtained so far in the design of the FBSO project and not all 
needed information is available at this time. For a list of 
all ground rules and assumptions made in determining the LVCM? 
see Appendix A . It is anticipated that ? as the Space 
Transportation System <STS) progresses? these data will become 
available to support a more accurate cost estimation. 

After all assumptions were established the LVCM needed to 
be updated to represent current dollars? since all the 
recommended CERs were determined using $88. This was 
accomplished by introducing an inflation factor of 1.889. Upon 
using this value all DDT & F and TFU cost estimation were 
automatically updated to $86. This factor was found using the 
NASA R&cD Inflation Indices provided by MSFC. 

Another item of importance of the LVCM is a factor 
titled? Complexity Factor <CF>. These CFs have a direct 
multiplicative effect upon the cost and determine the 
complexity of the structure to that of normal historical 
experience. It has been determined from previous cost models 
from NASA’s MSFC? the CFs were typically found to be 0.75 < 

1.0 <1.25. Since at this point in time little is known on the 



complexity of the structure* a CF of 1.0 was used. 

Once the basic procedures are defined for the cost model* 
the application of the CERs to obtain both DD & E and TFU cost 
is straight forward. First the analyst obtains both D?xD and 
FH&cA subtotals by exercising the subsystem-level CERs 
contained in the LVCM . Using these subtotals as a starting 
point* the analyst then estimates the total DDT & E and TFU 
cost by varying the performance characteristics (subsystem dry 
weights). Subsequently* this produces an accurate 
representation of the final cost estimation of the FE<S0 
project . 


ft- 


CONCLUSIONS 


The results of each of the subgroup’s studies have been 
tated. From the trajectory analysis;, it is evident that the 
>TS achieves a substantial altitude of nearly 52 nautical 
Ves before staging occurs. This high altitude tended to 
ate some difficulties in the return trajectory analysis, 
booster coasts to an apogee altitude of nearly 138 
cal miles. Due to low density values at such high 
\des. the booster experiences little aerodynamic drag and 
"ates rapidly during the ballistic reentry. As a result 
acceleration high G levels are experienced when the 
ed rate of descent is imposed. As mentioned earlier, 
'lent might be simplified by using a perigee injection 
' and should be examined in further detail in the 
?spite the minor problems, the booster’s return 
t was successfully modeled by the program DUAL. The 

b. ..as found to be capable of boosting either of the two 

second stages to a desired orbit and returning to the launch 
site for a horizontal landing. 

From the structures group it; is apparent that the shapes 
and dimensions of the booster’s propellant tanks can easily be 
determined. From the tank sizes the dimensions of the 
booster’s fuselage were dictated. The booster was calculated 
to be 193 feet in length with a diameter of 32 feet. The 
materials for the booster’s structure together with the 


allowable temperature extremes are also determined by the 
structures group. The booster is capable of completing the 
entire trajectory without exceeding the temperature extremes 
calculated . 

The performance analysis and configuration design of the 
booster is conducted by the aerodynamics group. The booster 
is very similar in design to the one created by NASA. Due 
this fact* the aerodynamic data used by NASA is also used in 
this analysis. The major difference in the design is the mass 
of the booster. The aerodynamics group gives a detailed 
report on the performance characteristics of the booster. It 
is apparent that the booster is capable of producing the 
required lift necessary for the return trajectory* despite its 
large mass of 800*110 pounds. An attempt is -made to calculate 
the required performance coefficients without the use of 
NASA’s values. This will be an area for future work. 

The propulsion groups successfully calculated the thrust 
requirements during both phases of flight* and found present 
engines to fulfill these requirements. Studies were conducted 
on the design and operating characteristics of the engines. 
Both the rocket engines and the air breathinq engines are 
capable of propel linq the booster as die: bated by the 
trajectory analysis. 

The final analysis deals with the weights and costs of 
the booster. Both areas of study are conducted with the use 
of NASA programs. The booster was broken into its various 


3 


components and each component examined for its value of 
weight. The weights were then totaled and found to match the 
880100 pounds dictated by the trajectory analysis. The cost 
analysis shows that the booster , despite its high cost, is 
finacially feasible. 

In an overview, the booster designers have generated the 
'■-* necessary parameters, equations and programs to provide the 

preliminary analysis of the flyback booster’s potential for 
design and development . This advanced booster design is found 
® to be very practical for the missions presently being designed 

in the space field. The booster is larger than expected, but 
the technology presently exists for it to be design to 
withstand the forces which result from its large mass. The 
results of this paper prove promising for future analysis of 
an Advanced Space Transportation System. 

@ 

© 

# 

! © 


96 



REFERENCES 


(1) Roy, A.E., Orbital Motion , John Wiley & Sons, 


(2) Hays, Mike, Trajectory Analyst, NASA Marshall Space 
Flight Center. 


(3) Limoge, Edmond J., Structures Analyst, NASA Marshall 

Space Flight Center. 

(4) Lockrow, A. Lynn, Model Design, Auburn University 

Auburn, Alabama. 

(5) Etkin, Bernard, Dynamics of Flight , John Wiley &= Sons, 

New York, 1982. 


(6) Lowery, Joe, Aerodynamics Analyst, NASA Marshall Space 
Flight Center. 


< 7 ) Struck, H.G., Butsko, J.E., "Booster Wing Geometry 
Trade Studies". 


(8) Olson, W.T., "Recombination and Condensation Processes 
in High Area Ratio Nozzles." Journal of American 
Rocket Society, Vol. 32 No. 5, May 1962, pps 672— 680. 


(9) Hill, F'.G. and Peterson, C.R., "Mechanics and 

Thermodynamics of Propulsion", Reading, Massachusetts 
Addison-Wesley Publishing Company, 1970. 


(10) Bilstein, R.E., "Stages to Saturn", Washington, D.C.: 
NASA, 1980. 


(11) Sutton, G.P., "Rocket Propulsion Elements", New York: 
John Wiley & Sons, 1956. 



( 12 ) Shevell, Richard S., Fundamentals of Flight, Printine- 
Hall, N.J., 1983. 


<13) Pryor, Mary, Ait — Breathing Engine Analyst, General 
Electric . 


< 14 ) Sanders, James L., Air-Breathing Engine Analyst, NASA, 
Marshall Space Flight Center. 


(15) Brothers, Bobby, Weights Analyst, NASA, Marshall Space 
Flight Center. 




(16) Hamaker, Joseph, Cost Analyst, NASA, Marshall Space 
Center . 


© 




© 

97 



APPENDIX A 

(Derivations and Assumptions) 



Aerodynamic Theory and Equations 
I.) Method for Finding Aerodynamic Center 

The aerodynamic center location of a wing-body configuration 
is given as follows: 

Xac ’ /Cre~ (Xac •' /Cre ) N+ ( Xac ’ /Cre ) W( B ) *CL W( B) + ( Xac ? / Cr e) B < W ) *CL. B < W 

CL N + CL W ( B ) + CL B(W> 

where 

CL N = (CN )B * < *d**2 ) / (4*Sw) 

CL W(B) * KW(B) * (CL >e * (Se/Sw) 

CL B< W) = KB ( W ) * (CL )e * (Se/Sw) 

Here KW(B) and KB(W) are interference factors obtained from 
Figure (27). The subscript W(B) refers to the exposed wing in 
the presence of the body and B(W) refers to the body in the 
presence of the wing. 

The static margin, (hn - h), is used as indication to the 
stability of an aircraft, where!! 

hn is defined as the non-dimensional distance from 
neutral point or aerodynamic center to some 
reference point. 

h is defined as the non-dimensional distance from 
center of gravity to the same reference point. 

If (hn — h) is positive the plane is statically stable. 

II>) Method for Finding Moment Coefficient 

The following equation is used to find the pitching-moment 
coefficient based on the product of the planform area and root 



dkiginae page is 

op POOR QUALITY 



body diameter d 
wing r-pan h 


FIGURE 27. Interference Factors 


Cm 


CmNO + < Xm/Cr 


Xcp/Cr ) * CN ’ 


< 6 ) 


where 

CmNO is the zero normal force pitching moment 

coefficient. For symmetrical configurations* 
as in our case* CmNO - 0.0 

Xm/Cr is the distance from the nose of the 

configuration to the desired moment reference 
center measured in root chords 

Xcp/Cr is the distance from the nose of the 
configuration to the center of pressure 

CN’ is the normal force variation with angle of 

attack 


III.) Method for Calculating Subsonic Character istics 

Due to the numerous equations necessary to write the 
program* only the most important equations will be mentioned 
here. The drag coefficient at a given angle of attack is given 
by : 

CD = (CDO)wb + (CDi)wb + < CD) mi sc (1) 


where 


<CD)misc takes into account drag contributors. 


The 



parasite drag coefficient for the wing body combination is found 
from: 

< CDo ) wb = < CDo ) w + (CDo)b * Sb / Sref <£> 
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The induced drag coefficient for the wing body combination is 
found from: 

< CD i ) wb = ( CDi ) w + CCD( ) 3b * Sb / Sref (3) 

After calculating CD, CL was found by rearranging the equation: 

CD = CDob + CLE / < pi * A * E > <4> 

to get: 

CL == sqrt C < CD — CDowb ) # p i * A * E (S> 

IV.) Method for Calculating Supersonic Characteristics 

The main equation for this calculation is given by* 

CL = < CL >’e CKN + KW(B> + KB < W ) 3 * Se ’ /S ’ 

+ (CL )”i CKW(B) + KB<W)3' , q ,!, S ,, 'Se'! 

q S’ S” 

+ (CL ) W •’ ' ( v ) 

where 

CL is the lift-curve slope of the wing-body tail 


/crz- 



combination 


(CL ) e 


q : 
q 


s 

Se 
< ’ ) 

< ' ’ ) 


is the lift-curve slopes of the exposed 
forward and aft surfaces. 

is the averaqe dynamic-pressure ratio actinq 

on the aft surface. 

is .the projected wing area. 

is the exposed wing area. 

denotes the canard. 

denotes the wing. 


The equation of supersonic zero lift drag coefficient used 


was: 


CD = (CDo)WB + (CDL)WB 


a 




where 

CD is the draq coefficient at supersonic speeds 

(CDo)WB is the zero lift drag of the winq body 

configuration 

(CDL)WB is the draq due to lift of the. wing body 

configuration 
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LIQUID ROCKET ANALYSIS 
ASSUMPTIONS & DERIVATION OF EQUATIONS 


The use of ideal rocket analysis takes into consideration the 
following assumptions for the thrust chamber: 

1. The working fluid (propellant products) is a 
perfect gas of constant composition. 

E. The chemical reaction is equivalent to a 
constant pressure heating process. 

3. The expansion process is steady * one- 
dimensional and isentropic. 

4. Chemical equilibrium is established within the 
combustion chamber, and does not shift in the 
nozzle. 


With these assumptions, equations for the analysis of the 
chamber can be developed. Through the application of the 
energy equation: 

Q + W = < h+V 2 /E+gz ) pV dA (1) 

in combination with the heating process, the? following 
equation is obtained. 

Tt a =Tti + Qr/Cp <2> 

This equation can be further developed assuming an isentropic 
nozzle expansion. 

Ue 2 » 2 Cp Tt a C l-(Pe/Pt) T “ l ' T ] (3) 

Assuming that the nozzle exit pressure is equal to the 
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atmospheric pressure at sea level? the exit velocity can now 
be related to the specific impulse by: 

Ue = Isp*g. (4) 

thus? an expression is obtained for the total chamber 
pressure. 

Pta = pe/(l - Ue 2 / a g c Cp Tte>- r ' T - 1 <5) 

From the calculation of the total chamber pressure, the exit 
Mach number can be evaluated by: 

Pta / Pe = Cl + C(t - 1)/ S3 Me 2 }' r '"' r ” 1 (6) 


solving for Me, the following relation is found. 

Me * £<£/ t -1) C(Pt a /Pe) T -‘ /T -1 3 > 1 (7) 

After calculating the exit Mach number, the nozzle areas at 
the throat and exit can be determined using the law of 
conservation of mass. 

m = p V A = constant <8) 

assuming: p = P / R T and V = a M 

m=(P/RT)aMA (9) 

Putting both static condition pressures and temperatures into 
stagnation quantities using these equations: 

P = Pta Cl + C(T - 1)/ S 3M 2 }-’-'"'- 1 (10) 

T « Tte (1 + C(T - 1)/ S 3M ) (11) 

the mass flow rate can be written now in terms of stagnation 
conditions. 

m = (Pt / R Tt) a M A Cl + C(t - 1)/ E 1M*}*'*-* (12) 


0 


Using the definition of the speed of sound 


© 
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a 


( t g 


R T) 1 '* 


( 13) 


a - t g e R Tt (1 + C(t - 1)/ £3M 2 3- 1 ''** (14) 

the mass flow rate can be written as: 

m = Ft A (t g c /R Tt) 1 '® £ 1 + C ( t-1 ) /23M 2 y*+*- '« < •*— 1 > (15) 
Realising the mass flow rate is: 

m = (T / Isp) ( g<= / g- ) (16) 

the area can be evaluated in terms of the mass flow rate. 

A ■ m /(Pt a M)*(R TWt g«= ) £1 +C(t -1)/ 23M 2 3- (17) 

After further algebraic manipulations? and using throat 
qualities of M 8 1 and A* = A* the throat area is equal to the 
following expression. 

A* = (m/Pta) CR Tta/g«= T ( 2/t + 1 ) 1 3 1 "® (18) 

Now using the ratio of mass flow rates at the throat and the 
exit an equation results which relates the exit area to the 
exit Mach number and A* 

Ae = (A~/Me) C(2 /t+1) £1 + C(t - 1 ) /23Me 2 > (19) 

From these areas the diameters of the throat and exit can be 
calculated . 

d = 2 (A/w) 1 '® (SO) 

Thrust chamber characteristics can be written in terms of the 
char acter i st ic velocity.* C** ? and the thrust coefficient? CT . 

C~ ■ (Pt a A* ) /(m) (SI ) 

and 

CT = T/(Pt A* ) (SS) 
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GROUND RULES AND ASSUMPTIONS 

3 

FOR COST ESTIMATION 

3 

* All cost and equations are shown in 1986 dollars. 

(in mill ions) . 

* The complexity of the structure has the same complexity of 

3 

normal historical experience. 

* Flight Hardware and Assembly costs are included in the 
Flight Hardware cost at the subsystem level. 

* Systems Test Hardware and Assembly costs, Structural Tooling 
and Systems Test operations are not addressed in the cost 
model at this time. 

Q ^ ..... . 

* Slopes for the FBSO subsystems were assumed by utilizing 

data from Aerospace Corporation SRM cost model and F'RC 
developed CERs from other cost yodels. 

O-. 

* The inflation factor is based on NASA Comptroller Indices, 
date January 1986. 

* An historical cost data utilized in this model are contained 

© 

within the REDSTAR Data Base System. 

* No institutional cost are included. 

* Contractor (Prime) fee and/ or profit is excluded. 

0 

© 


/© 7 


© 



APPENDIX B 


Prograrnsj Input St Output Files 



ALT -0 PlLOGrRAAl | 

DT -0 ♦ 1 __ ^ 

RANG=0 TVtAAEtf.'toft.Y “To 

ISPl ss 320 
ISP2=380 

INPUT 'ENTER INITIAL MASS' ?M1 

M'=M1 

INPUT 'ENTER ALTITUDE FOR 0 -TURN ' 5 A0T1 
i'.oT-AGTI 


ORIGINAL PAGS jg 
OF POOR QUAXH2J 
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so 

VO- 

©o 

GOT 


FRI 

?0 
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VC t DC .: ZD' 


' ~“Y+0 ♦ 5* ( VO+V ) #DT 
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I L- i. i I w *- ! v i i i 
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> v *•* ■* i . i ^ 
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RANG 


6 ' S 



20 

jO 

;o 

M-M+MBQTlfcBT+MBGT2#BT 
BPSI”<PI/130>*0. 1 
PSIQ-PI/120 

ORIGINAL PAGE IS 
OE POOR QUALITY 



.0 

PSI^PSIO 




:.0 

'0 

GAM-1 

ZO=SIN<PSI0)/<l+CGS(P:;i0> > 




U) 

PRINT "TIME V RANG 

ALT PS I GAM THRUST 

MASS 

6 ' 3 • 

s° 

F$=*lt# ***** ******* 2 

****** **.*• **-* ******* 

***!i ii 3 / 

!* » »*■ ** 

J -1. 

PRINT 

N=TH/M 




;• w‘ 

PRINT USING F$JT*V> RANG? ALT 
OVO/ < Z0~ ( N~1 > ) / ( l-fZO-2 ‘ 

r* t ;•> ’«• . . r\ 

7 P 3 1 1 1 3 0 / P I v GAM / Til? M •/ A G 



. ./ 

i W J. " 1 w * v l O 
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-0 
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'I- T rL: 
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Jim W 
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•' 
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j. i a- ■' \j a. 

: « A. 1 ..f i a. . i l .i. / .<■ .* t .i. 




0 

\ ♦. .*% -V« f . • • 1 — • • .**V 

1 i x :• tJ ; — I i'*i / J. ..*• i ‘ 




••$’20 

M M - M BOTl^BT- M 1 « G T 2 % 0 2 




•• %.} VJ 

IF MAO THEN 1390 




;• a A 

G » 3 0 % 0 / ( R 0 f A L.T ) ) ~2 





GOTO 720 

vi „•». if "+ ■■ . \ • t » -v ^ -n r> • / - a j y t\ w ■ > 

•r > U / w 1 \ \ J. A,. \t »,. ) / S J. 1 i.. J i- -• 




• ’ w 

g-vo/z:; : ;n ) / ( i : Z , o*'2 > 





v~C#Z''(N“l )#( 1+Z'*2) 




, 7 C 

di=c/G3kz”< N-i ) % a/ ( N-i > v/rz. 

/ j. > > -o/ojjczo" ( m - 1 : # ■ 5 / ( 24-1 > z 



. 30 

DX=0 ♦ 5* ( V0*SIN< PSIO ) + V %% I K < ! 

r<: 2 ;■ r 



* ~ 

BY =0 ♦ 5* ( VO*COS (PSIO) + V#C CS ( ! 

Poi:>*b; 



. Z w‘ 

X-X+BX 




,."l \/ 

RANG-RO*THETA 




Svo 

Y=Y+DY 




. •' j 0 

ALT « •' V fLO ) / COS ( THETA ) - R0 




| O 0 

T ~TtBT 




70 

PSIO-PSI 




■ O 

Z 0 z 




• r\ r> % 

V-V:V 




©00 

TH-=N»M 



| 

: i o 

THl«TSt~TH2 




.220 

MBCT1~TH1/] 3P: 




230 

M B 0 T 2 -- T H 2 / 1 S P 2 




240 

M * M ~ M B 0 T 1 * B 7 • M D 0 2' 2 * P T 

HO 



.230 

MC1»MC1+MB0T1*BT 



©60 

MC2-MC2+MBCT2'-KDT 




270 

MC-MC1+MC2 




2S0 

IF M<0 THEN 1390 





i‘\ / r*. « i i 3 g > v • US 





500 

510 

520 

530 

550 

560 

570 

580 

590 

591 

592 
500 

520 

521 

522 

523 

524 
K'5C 


PSI*PSI+DPSI 

THET A-ATN < X/ • ROT Y ) ) 

GAMMA-PSI-THETA 

GAM=GAMMA#1S0/P1 

PRINT USING F$» T » V » RANG » ALI ? PS I •> GAN » TH » M y AG 
IF V>7000 THEN 1390 
Z=SIN<PSI ) / < 1+COS < PSI ) > 

GOTO 1060 
PRINT 

PS I *PSI -DPS I 
PRINT VI 

PRINT 'TIME V RANG ALT PSI 

PRINT* 

PRINT 


ORIGINAL page is 

OP POOR QUALITY 


GAM 


THRUST MASS 


PRINT* 
PRINT 
PRINT* 

PRINT * Liiiie 
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t /•' r» ... - ,■% r • 


tsT, 


j 

>•: ■ .• • 
• 


: • ( : ■ >.j ' > ;••• p j* / * y . j- j ; j p t a 

■ '< 'A < ’• •> V *■' '«/ *♦ “*« / f. T 

* * i i ; i f : i *•/ ••• / r j. 

'*• .. r: •? >.• ■- c s t a ' * 

T ! i / » 

r- . I > I {■ . «•*,!:* *** * * It 'in ,*• • T . »;•» r? 1 - t O /y / * 

i ‘ 1 i < ; u 1 .k i •. v > »• / » »" 5 j*’*. ' 1 1 /“Hu I r t 55 .■ v* .< .* 

r : v "• / ‘ “ ^ i ; / ( 

, .. ... ...* ... . . . -r 

• k* • ■ ; ■—* .. V.* 4 . 

• 1 / * V I “V * I 

• i t S i i i « I I — 111 

•• - •/ i ,•> /• / • T 

I" Wl"t ’t.'l' I 'Jk v ,• * -I. 


AM Ms 


* 1.1 . M , 
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ALl « 5 Y PRO ) .'COS ( THETA) -R0 
. T v-d r 


%: z '5 ( N - j ; < i / ( m - j '• 5 7 . 0 ” 2 / ‘ 


LT * 
SO 

•31 

; :o 

L 10 

.. 5. 2 


•. MDC'I “ D ? 

I ,1<C THEN 2080 
0-3 OJK < RO/ < ROT ALT ? > ~2 

r* j”* M 
r*» r;« 7 1 -j 

PRINT * INITIAL MASS 
PR I'P " Ai ! FUR G TURN 


I " * M 1 

* m ' n n t •» 

.• ri •../ ! j. 


HI 


r*,* r-o'n \ *rruT i 'j r.:- r ♦ n * 

Wi’IHUu •/ w. I ! . 1 . *V Ll- ■■ V / I . I A. J. 


p R I N T * f- 1 : 5 A L V ELOCITY ♦ * $ V 

!'r:;:h t * final macs :*jm 

PR INI "FINAL PSI J * » PSI* 130/P I 


t ttim ir.r 


* 1**1 


* / « HZ' 


n y i 


PRINT'FINAL GAMM J-JUAMM 
NP'JT’TRY AGAIN Y^l N=2"*TRY 
F TRY <2 THEN 15 


ORIGINAL PAGE IS 
OF POOR QUALITY. 





PROGRAM DUAL 


26 AUG 1986 


IMPLICIT DOUBLE PRECISION <A-HfO-Z> 

INCLUDE 'DUAL.CMN' 

CHARACTERS IDUM 
CHARACTERS HNM ( 5 ) 

CHARACTER*12 FNMfGNM 

DATA HNM/' PLOT 1 'f'PL0T2 'f'PL0T3 'f'PL0T4 


'f 'PLOTS '/ 


WRITE< IOUT f200) 

READ <IINf210> FNM 
WRITE( IOUT f220> 

READ( IINf210) GNM 

OPEN < UNIT*ICR f FILE*FNM f STATUS® ' OLD ' ) 

OPEN ( UNIT*IPR f FILE*GNM f STATUS* ' UNKNOWN ' ) 

OPEN(UNIT*IMTfFILE*'SCRATCH' f STATUS* 'SCRATCH' f FORM* 'UNFORMATTED' ) 
WRITE< IOUT f230) 

READ< IINf210) DATE 


CALL EDATA 
NPL0T*0 

10 DO 20 1*1 f96 ! INITIALIZE INTEGRATION ARRAY 

20 A<I)“ZERO 
REWIND IMT 


READ TITLE CARD 1 

READ OPTIONS CARD 2 

POSITION AND VELOCITY CARD 3 

FLY ALTfENG ON ALTfGUIDE CHECKS CARD 4 

ALP fALPCf ALPHA PROFILE CARD 5 

ALIMfALPHADOT PROFILE CARD 6 

BET f BETC f BETT f BETM f BETL f BETX CARD 7 

BLIMfBETADOT PROFILE CARD 8 

THRUST PROFILE fENGMfENGF CARD 9 

WEIGHTS-EMPTYfFUELfON ORBITfRESERVE RATES-APUfVENT CARD 10 

> ACC LIMITfAERO PERTRBfFLY THRST PERTRBfSITE 8 CARD 11 

TIMEfDT fFLY DTfPRINT TIMEfDTPfFLY DTP CARD 12 


READ < I CR f 21 0 f END®50 ) TITLE 
WRITE ( IOUT f215) TITLE 

READ < ICR f 100 ) NABT f NATM f NFL Y f NLIM f NLND f NPRN f 
3 1NPLT f NPUT f NQAL f NTR Y fNTYPfNVNTfNWINfNMET 

READ < ICR f 101) FLAT f FLON f R f V f GAM f HDNG f 
1 ALTH f ALTE f RDTN f TCON f TSTG f VTST f 
2ALF f ALFD f BETA f BETAD f THRF f 
3WGTE f WGTF r WGTO f WGTR f FMD A f FMD V f 
3 4ACCL f DCDf DCL f DTH f XLL A < 8 ) f XLLO < 8 ) f 

5T f DT f DTF f TP f DTP f DTPF 

IF <NABT »EQ»2) FMDV*ZERO ! CHECK CROSS COUPLING 

IF (NLND»LE«0) NLND*1 
IF ( FMDA ♦ GT ♦ ZERO ) FMDA=-FMDA 
@ IF ( FMDV ♦ GT ♦ ZERO ) FMDV=-FMDV 
IF < FMDV ♦ EQ ♦ ZERO ) NVNT=1 


NPG*1 

CALL HEADER 

WRITE( IPRf217) FNMfGNM 
IF < NPLOT ♦ GE ♦ 5 > NPLT=1 


! RESET CASE PAGE NUMBER 
! PRINT TITLE AND CASE DATE 

//J 


000010 

000020 

000030 

000040 


000050 

000070 

000080 

000090 

000100 

000110 

000120 

000130 

000140 

000150 

000160 

000170 

000180 

000190 

000200 

000210 

000220 

000230 

000240 

000250 

000260 

000270 

000280 

000290 

000300 

000310 

000320 

000330 

000340 

000350 

000360 

000370 

000380 

000390 

000400 

000410 

000420 

000430 

000440 

000450 

000460 

000470 

000480 

000490 

000510 

000520 

000530 

000540 

000550 

000560 

000570 

000580 


i fs i*. t r>n 


000590 


grK a irri 

DO 25 J»1f2 

WRITE < JPR r 102 ) NABT f NATM f NFLY f NLIM f NLND f NPRN r ! OPTION TABLE 000600 

1NPLT f NPUT f NQAL f NTR Y f NTYP f NVNT f NWIN f NMET 0006 1 0 

25 JPR-IOUT 

WRITE( IOUT f 140) 

READ •< IIN f 210) IDUH 
WRITE< IOUT f 215) TITLE 

000620 


WRITE ( IPRf 103) FLATfFLONfRf VfGAMfHDNGf ! PRINT INPUT 000630 

1 ALTH f ALTE f RDTN f TCON f TST6 f VTST f 000640 

2 ALF f ALFD f BET A f BET AD f THRF 000650 

WRITE ( IPRf 104) WGTEfWGTFf WGTOf WGTRfFMDAfFMDVf 000660 

1 ACCL f DCD f DCL f DTH f XLLA ( 8 ) f XLLO ( 8 ) f 000670 

2T fDT fDTFfTPfDTPfDTPF 000680 

) 000690 

CALL HEADER 000700 

WRITE < IPRf 105) 000710 

000720 

TBEGIN-T ! SAVE INITIAL TIME 000730 

CALL TYPE 000740 

* CALL RDWIND 000750 

CALL INITIAL 000760 

CALL FUN 000770 

IF(NPRN»EQ«3) CALL HEADER 000780 

CALL DRIVER 000790 


IF ( NPLT * EQ * 1 » AND ♦ NPRN » EQ * 2 ) GO TO 40 
IF<NPLT *EQ»2) THEN 
NPL0T-NPL0T+1 

OPEN < UNIT-IFL f FILE-HNM < NPLOT ) f ST ATUS* ' UNKNOWN ' ) 

WRITE<IFLf ' ( A60f A12) ' ) TITLEfDATE 
IPTS-23 

WRITE < IFLf ' ( 16) ' > IPTS 
WRITE < IFLf ' ( 10A8) ' ) 'TIME-SEC' f 'LONG-D 'f'GCLAT-D 'f'ALT-KFT 'f 

1 'ALPHA-D 'f'BETA-D 'f'ACCT G"S' f 'DYNP-PSF' f 'Q-BTU/S 'f 

2 'QSUM-BTU' f 'MACH 'f'VREL-KT ' f ' THRUST-P' f ' WEIGHT-P' f 

3 'FUEL-P/S' f 'RANGE-NM' f 'FPA-D 'f'AZM-D 'f'ACCXG"S'f 

4 'WPOS G"S' f'ACCT I 'f'THR-% 'f'Q*ALPHA ' 

END IF 


ORIGINAL PAGE IS 
OF POOR QUALITY 


REWIND IMT 000830 

NBLKS-15 000850 

30 READ ( IMT fEND*40) Tf ACCf ACCXf ACCZfDACZf ALPf ALTfBETfCDfCLfCRNGf 000920 

1 CTMP f DRTG f DYNP f FLAT f FLON f GAM r GAMR f HDNG f NHDP f QALF f QDOT f QTOT f 000930 

$ 2RD0T f RHO t V r VMACH f VRL f WGT f WDOT f XDRGM f XLFTM f XTHRM f THRTL f ACC I 000940 

000950 

IF (NPRN»EQ«2) GO TO 34 

IF < NBLKS » GE » 10 ) THEN * 000860 

NBLKS »0 000870 

CALL HEADER 000880 

3 WRITE < IPRf 106) 000890 

END IF 000900 

NBLKS*NBLKS+1 000910 

WRITE < IPRf 107) Tf ALTf V f GAM f DYNP f QDOT f ALPf FLAT f DRTG f ACCZ 000960 

WRITE (IPRf 108) NHDPfRDOTf VRL f GAMR f VMACH f QTOT f BET f FLON f HDNG f QALF 000970 
% WRITE < IPRf 109) RHOfXTHRMf WGTfUDOTfCDfCL 000980 

001040 

PLOTTING 


34 IF (NPLT »EQ«2) THEN 
ALT-ALT/1000* 

VRL-VRL/1 *68781 ,.J 

GPOS^-ACCZ * 

WRITE ( IFL f ' ( 5E16 * 7 ) ' ) T f FLON f FLAT f ALT f ALP f BET f ACC f DYNP f QDOT f 

a#nrnT _ LLUxru_ i iqi _ vtucu _ Lir**r _ unnT. nnro _ camd _ unxift .Arrv.nDnc.ArrT . TUDTI . 


lull V rinwri 9 VIM. » A I I ll\n t WO I > WOO I t Ol\ I O r oniHS / » I*myw r nww/% r w» ww r riww* / 

*QALF 
END IF 

60 TO 30 
40 CALL HEADER 

DYNPM*D YNPM/47 ♦ 8802S ! LBS/SQFT 

QALFM*CNV#QALFM/47« 88025 • DEGS*LBS/SQFT 

FUEL-WGTE-WGT 

JPR=IPR 
WNEG=-ACZPM 
WP0S--ACZNM 
DO 45 K=1 f 2 

WRITE < JPR f 1 10 ) QD0 TH f TQD0T t CTPM f TCTPM f STPM f TSTPM f DYNP h f TDYNPM f 
1QALFM t TQALF f ACCTM f TACCT f ACCXM f TACCX f ACCYM f TACCY f 
2 ACCZR t TACCR f WNEG * TACZP f WP0S f TACZN 

IF<NFLY*NE»1) WRITE < JPRf 130) FUEL 


45 


JPR-IOUT 

CONTINUE 


CLOSE < IFL) 


ORIGINAL PAGE IS 
OF POOR QUALTE 33 


50 

100 

101 

102 


IF<NPLT »EQ«2) 

GO TO 10 
STOP 

FORMAT <2012) 

FORMAT <6E13«8) 

<//52X» 'OPTION TABLE'// 

NABT *f I4f4Xf 'MAIN ENGINES OFF / MAIN ENGINES ON'// 
=f I4f4Xf 'PRA63 ATMOS / STD USA2 ATMOS'// 


* 


FORMAT 
130Xf8H 
230Xf8H 
330Xf8H 
430X f 8H 
530X * 8H 
6 

730Xf8H 
830X f 8H 
930Xf8H 
*30Xf8H 
130Xf 8H 
230Xf8H 
430Xf8H 
530X r 8H 
630Xf8H 
103 FORMAT 
> 16H 

26H 
36H 
46H 
56H 
66H 
76H 
86H 
96H 
*6H 


NATM 
NFLY 
NLIM 
NLND 
/f46Xf 'NAS 
NPRN 
NPLT 
NPUT 
NQAL 
NTRY 
NTYP 
NVNT 
NWIN 
NMET 


=fI4f4Xf 'NO FLYBACK 
=fI4f4Xf'N0 G LIMIT 
«fI4f4Xf 'CAPE KEN / 
CECIL / CHERRY 
FULL PRINT 


/ FLYBACK / FLYBACK & G.S«'// 

/ EXP LIMIT / FIXED LIMIT'// 
HOUSTON / EDWARDS / HICKAM' 
POINT / AZORES / INPUT'// 

*f I4f4Xf 'FULL PRINT / STATE ONLY / TABLE ONLY'// 
=fI4f4Xf'N0 PLOT FILE / WRITE DATA TO PLOT FILE'// 
=fI4f4Xf 'INERTIAL / INERTIAL AND ENC / RELATIVE'// 
*fI4f4Xf'N0 LIMIT / Q ALPHA LIMIT'// 

=f I4f4Xf 'BEGIN AT THIS SEGMENT'// 

■ t 14 f4Xf 'BOOSTER / ORBITER'// 

=fI4f4Xf'N0 VENT / VENT'// 
m r I4f4Xf 'NO WIND / 95% INPUT / 95% 

*fI4f4Xf 'INPUT ENGLISH / METRIC') 


IN+HEAD' // 


</f5Xf 'INPUT VALUES'/ 


LATF12*4f6H 
VVVF12*3f AH 
ALTHF12 ♦ 2f 6H 
TC0NF12*3f6H 
ALPF12 ♦ 3 f 6H 


L0NF12»4f6H ALT/RF12.2» 
GAMF12* 3 f 6HEN/HDNF12 . 3/ 
ALTEF12*2f6H RDTNF12»2f 
TSTGF12.3f6H VTSTF12.3/ 
ALPCF12*3f4(6XfF12*3>/ 
ALIMF12»3f5(6XfF12»3)/ 

BETF12»3f6H BETCF12»3f6H 
BETMF12»3f6H BETLF12,3f6H 
BLIMF12*3f5(6XfF12«3)/ 
THFACF12»3f3<6XfF12*3) f6H ENGMF12.3f6H 


BETTF12*3f 

BETXF12.3/ 


ENGFF12.3) 


104 FORMAT < 


© 


WGTFF12»2f6H 
FMDAF12»2f6H 
DCDF12.2f6H 
LAT8F12»2f6H 
DTF12»2f6H 
DTPF12*2f6H 
105 F0RMAT<15Xf 'OUTPUT UNITS'// 

© 1 IOXf 'TIME' f16Xf 'SEC'/ 

2 IOXf 'ALTITUDE' f12Xf 'FEET'/ 

3 IOXf 'VELOCITY' f12Xf 'FEET/SEC'/ 


16H 

26H 

36H 

46H 

56H 

66H 


WGTEF12»2f6H 
WGTRF12* 2f6H 
ACCLF12 ♦ 2f 6H 
DTHF12*2f6H 
TF12»2f6H 
TPF12*2f6H 


WGT0F12.2f 
FMDVF12 ♦ 2/ 
DCLF12.2f 
L0N8F12 » 2/ 
DTFF12»2f 
DTPFF12.2) 


ft* 


AA.uet.rn a av ,/ nr.en — .. 


001060 

001070 

001080 

001090 

001100 

001110 

001120 

001130 

001140 

001150 

001160 

001170 

001180 

001190 

001200 

001210 

001220 

001230 

001240 

001250 

001260 

001270 

001280 

001290 

001300 

001320 

001330 

001340 

001350 

001360 

001370 

001390 

001420 

001430 

001440 

001450 

001460 

001470 

001480 

001490 

001500 

001510 

001520 

001530 

001540 

001550 

001560 

001570 

001580 

001590 

001600 

001610 

001620 

001630 

001640 

001650 

001660 
l\f\ 1 


H 1UAF HNULtd' F XHA F ' UCUS / 

5 IOXf 'RATES' f15Xf 'FEET/SEC'/ 

6 IOXf 'RANGES' t 14Xf 'NAUT HILES'/ 

9 7 IOXf' WEIGHT ' f 14X f ' LBS ' / 

8 IOXf 'THRUST' f14Xf'LBSF'/ 

9 IOXf 'HEAT IN'f13Xf'BTU/FT**2'/ 

* IOXf 'HEAT RATE' f IIXf 'BTU/FT**2/SEC'/ ORIGINAL* PAGE IS 

1 IOXf' PRESSURE ' f 12X f ' LBS/FT*2 ' / OF POOR QUALITY 

2 IOXf 'ACCELERATIONS' f7Xf'Gs'/ 

3 IOXf 'Q ALPHA ' f 13X f ' DEGS#LBS/FT*2 ' ) 

106 FORMAT ( 

* //f40Xf 'ENGLISH UNITS' 

1 /f31Xf 'TIME' f2Xf 'ALTITUDE' f3Xf'VNRT GAMNRT' f3Xf 'DYNP' f3Xf 

2 'QDOT' f2Xf 'ALPHA' f4Xf 'LAT' f3Xf 'DRTG' f3Xf ' ACCZ'/ 

3 41Xf'RD0T'f3Xf'VREL GAMREL' f3Xf 'MACH' f3Xf 'QTOT' f3Xf 'BETA' f4Xf 

? 4 ' LON ' f 3X f ' HDNG ' f 3X f ' Q ALF ' / 

5 42Xf'RH0'f8Xf 'THRUST' fHXf 'WGT' fIOXf 'UDOT' f5Xf'CD' f5Xf 'CL' ) 

107 FORMAT < /25XfF10»3fF10»0fF7,0fF7*2f4F7. 1 fF7»0fF7»2) 

108 FORMAT (33Xi A2 f3Xf 2F7»0f2F7.2fF7. 0f2F7. 1 f F7» 0 f F7 ♦ 0 ) 

109 FORMAT (35XfE10*3f2(5XfF9*0) f5XfF9«2f2F7 *3) 

110 F0RMAT(//15Xf 'MAXIMUM VALUES'// 

> *18Xf10HMAX QDOT fF7.3f16H BTU/FT2/S TIME fF8*1/ 

118Xf 10HMAX CTMP fF7*0f16H DEGS F TIME fF8»1/ 

218Xf 10HMAX STMP fF7»0f16H DEGS F TIME fF8*1/ 

318Xf 10HMAX DYNP fF7»0f16H LB/FT2 TIME fF8.1/ 

418X f 10HMAX QALF fF7»0f16H DJLB/FT2 TIME fF8»1/ 

518Xf 10HMAX ACCT fF7»3f16H GRAVS TIME fF8.1/ 

618X f 10HMAX ACCX fF7«3f16H GRAVS TIME fF8*1/ 

718Xf 10HMAX ACCY fF7»3f16H GRAVS TIME fF8.1/ 

818Xf 10HMIN ACZR fF7»3f16H GRAVS TIME FF8.1/ 

918Xf 10HUING -G fF7«3f16H GRAVS TIME fF8*1/ 

*18Xf 10HUING +G fF7*3f16H GRAVS TIME fF8»1> 

115 FORMAT <6Xf 'TIME' f7Xf ' VRL' f7Xf 'FPA' f2Xf 'ALTITUDE' f 
i 14 Xf 'ATTACK' f6Xf 'BANK' f6Xf 'DENSITY'/ 

25Xf ' (SEC) ' f2Xf ' (FT/SEC) ' f5Xf ' (DEG) ' f6Xf ' (FT) ' f 
35Xf ' (DEG) ' f5Xf ' (DEG) ' f6Xf ' (SLUGS) ' ) 

120 FORMAT(F10*3fF1043fF10»4fF10»2f2F10.3fE13.5) 

130 FORMAT (//15Xf 'FLYBACK FUEL' fF10.2//> 

140 FORMAT (//IOXf 'HIT RETURN TO CONTINUE') 

2200 FORMAT (//IOXf' INPUT FILE 'f*> 

205 FORMAT ( IX fA60) 

210 FORMAT (A) 

215 FORMAT ( / f 5X f A60/ ) 

217 F0RMAT(15Xf' INPUT FROM ' fA12f5Xf 'OUTPUT TO 'fA12) 

220 FORMAT ( IOXf 'OUTPUT FILE 'f*) 

3230 FORMAT ( IOXf 'DATE 'f*) 

END 

SUBROUTINE DRIVER 

IMPLICIT DOUBLE PRECISION (A-HfO-Z) 

INCLUDE 'DUAL»CMN' 

EXTERNAL FUN 

ILINE=' BEGIN ENTRY ' 

NHD=2 

NHDP=2HBE 


! VENT CONTROL 


//(, 


IF(NABT »EQ»2) GO TO 3 
® IF(NVNT »EG» 1 ) GO TO 6 
TT=-WGTR/FMDV 
TV*T+TT 

IF(TT*GT»ZERO) GO TO 4 
3 IF (FMDV»NE»ZERO) THEN 
NHD*2 

© NHDP=2HVN 

ILINE='END VENT 
FMDV=ZERO 


001680 

001690 

001700 

001710 

001720 

001730 

001740 

001750 

001760 

001770 

001780 

001790 

001800 

001810 

001820 

001830 

001840 

001850 

001860 

001870 

001880 

001890 

001900 

001910 

001920 

001930 

001940 

001950 

001960 

001970 

001980 


001990 

002000 

002010 

002020 

002030 

002040 

002050 

002060 

002070 

002080 

002090 

002100 

002110 

002120 

002130 

002140 

002150 

002160 

002170 

002180 

002190 

002200 



tnu ir 
NVNT-1 
GO TO 6 

? 4 IF(TV-T.GT.DT) GO TO A 

NVNT«3 

FMDV® ( TV-T ) *FMDV/DT 

6 IF<NABT*EQ*2> GO TO 30 
. 10 CALL GUIDE 

* CALL FUN 

' CALL TSTMAX 

IF<NEND«EQ»2) GO TO 40 
20 CALL PRINT 
TF=T+DT 

IF<TF*GT»TP> TF=TP 
> CALL RUN< A»T » TF fFUNf 7 ) 

IF< ALT *LT* 10* ) THEN 
ILINE* 7 CRASHED 
NHD®2 
NHDP*2HZZ 
GO TO 40 
END IF 

GO TO < A»4»3) fNVNT 

30 GO TO ( 31 f 32 » 34 ? 36 ) f NSTG ! THRUST CONTROL 

31 NSTG*2 

ILINE® 7 THRUST EVENT 7 
? NHD=2 

NHDP«2HTE 

WTST*WGTE+WGTR 

32 IF(TSTG*GT»T) GO TO 34 
ENGM®ENGM/TWO 

NSTG-3 
$ NHD*2 

NHDP*2HTE 

34 FISP*FISPL+BEX#DELP 

THRST® ( THRSL+AEX*DELP ) *THF AC*ENGM 

FMDM*-THRST/FISP 

WGT*FMASS*GO 

3 IF <WGT+<FMDA+FMDM)*DT.GT.WTST> GO TO 3? 

NSTG*4 

FMDM* ( WTST-WGT ) /DT-FMDA 
THRST«-FISP*FMDM 
GO TO 39 
36 NABT*1 
:§ FMDM«ZERO 
THRST-ZERO 

IF ( WGTR ♦ LE ♦ ZERO ) GO TO 38 

FMASS*FMASS-WGTR/GO ! DROP PERFORMANCE RESERVE 

ILINE* 7 WGT DROP 
O NHD*2 

NHDP=2HWD 
GO TO 39 

38 ILINE® 7 THRUST EVENT 7 
NHD=2 
NHDP=2HTE 

© 39 XTHRM=THRST /FMASS 
GO TO 10 
40 TP=T 

CALL PRINT 

RETURN 

END 

@ SUBROUTINE INITIAL 

IMPLICIT DOUBLE PRECISION <A-HfO-Z> 

INCLUDE 7 DUAL*CMN 7 //? 


! CONTROL LOOP 
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002220 

002230 

002240 

002250 

002260 

002270 

002280 

002290 

002300 

002310 

002320 

002330 

002340 

002350 

002360 

002370 

002380 

002390 

002400 

002410 

002420 

002430 

002440 

002450 

002460 

002470 

.002480 

002490 

002500 

002510 

002520 

002530 

002540 

002550 

002560 

002570 

002580 

002590 

002600 

002610 

002620 

002630 

002640 

002650 

002660 

002670 

002680 

002690 

002700 

002710 

002720 

002730 

002740 

002750 

002760 

002770 

002780 

002790 

002800 

002810 

002820 

002830 

002840 

002850 

002860 

/wYop/7n 


V v*-w / w 


TACCR“T 



TACCT»T 

TACCX-T 

TACCY-T 

TACZN-T 

TACZP-T 

TCTPM“T 



TDYNPM“T 

TGDOT“T 

TSTPM“T 

ORIGINAL PAGE IS 


ACCTM=ZERO 

OF. POOR QUALITY 

002910 

ACCXM=ZERO 

002920 

ACCYM=ZERO 

ACCZR“ACCL 


002930 

ACZPM“ZERO 


002940 

ACZNM* ONE 


002950 

CTMP=ZERO 


002960 

CTPM“ZERO 


002970 

DACZ“ACCL 

* 


DYNPM=ZERO 

QALFM“ZERO 

QDOTM“ZERO 


003010 

STPM“ZERO 


003060 


NBET=1 

NBLP=4 

NEND=1 

NSTG=1 

NTF*1 

NTR“1 

3 IF<NMET»EQ«1) THEN 
R=R*FTM 
V»V*FTM 
ALTH«ALTH*FTM 
ALTE*ALTE*FTh 
RDTN=RDTN*FTM 
END IF 

DO 2 I“l>6 
ALF < I ) “ALF < I ) /CNV 
ALFD ( I ) “ALFD < I ) /CN V 
BETA < I ) “BETA ( I ) /CNV 
3> 2 BET AD ( I > =BET AD ( I ) /CNV 

ALP“ALF < 1 ) 

SALP=SIN< ALP) 

C ALP=COS < ALP ) 

BET=BETA< 1 ) 

3 SBET=SIN<BET) 

CBET “COS < BET ) 

DALP=TWO/CNV 

XLAT=XLLA < NLND ) /CNV 
XLON“XLLO < NLND ) /CNV 
® SLLAT“SIN<XLAT> 

CLLAT “COS < XL AT ) 

RLN=RE/SQRT ( 1 ♦ +GF#SLLAT*SLLAT ) 
XL ( 3 ) =RLN*SLLAT 


003150 

003160 

003170 

003180 

003190 

003200 

003210 

! CONVERT ENGLISH TO METRIC 003220 

003230 

003240 

003250 

003260 

003270 

003280 

003290 

! CONVERT CONTROL ANGLES 003300 

003310 

003320 

003330 

003340 

003350 

! SET INITIAL CONTROL ANGLES 003360 

003370 

003380 

003390 

003400 

003410 

003420 

003430 

! LANDING SITE 003440 

003450 
003460 
003470 
003480 
003490 
003500 
003510 
003520 
003530 
003540 


© 


FLAT-FLAT/CNV 
SLAT=SIN<FLAT> 
CLAT “COS < FLAT ) 


INITIAL LATITUDE 


Ilf 


a i mi — p*i nil /mu i iru»nurn> 


i nkif?TTi mr ad meTrn 


nk.uix—r t_um/ uihvt I ^urit-uri : lui'iua iu*.'u nt'uuu 

SLON*SIN< ALON) 

CLON*COS ( ALON > 

GAM*GAM/CNV I FLIGHT PATH ANGLE 

,SGAM*SIN<GAM> 

* CGAM*COS ( GAH ) 

HDNG=HDNG/CNV ! HEADING (AZIMUTH) 

SHDN*SIN(HDNG> 

CHDN=COS < HDNG ) 

NPUT* 2f HEADING COMPUTED ASCENDING FROM INPUT INCLINATION 
-2» HEADING COMPUTED DESCENDING 
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! OBLATE EARTH RADIUS 

! ALTITUDE INPUT 
! KM 


! INITIAL POSITION 

X<2)=R*SL0N#CLAT 

X(3)=R*SLAT 

XDA(1)*-X< 2) IOMEGA ! ATMOSPHERE VELOCITY VECTOR 

XDA ( 2 ) * X<l)*OMEGA 
XDA<3)» ZERO 

CALL COMPASS 
DO 4 1*1 f 3 

4 XD<I)=CHDN#AA<I)+SHDN*BB<I> 

CALL VUNIT (XDrXD) 

DO 6 1*1 >3 

XD < I ) * < CGAMftXD < I ) +SG AM#CC ( I ) ) * V 
6 XDU( I >*ZERO ! NO WIND 

DIR*ZERO 
CDIR=ONE 
SDIR=ZERO 
VWIN*ZERO 

ALT-R-REI 

IF (NWIN»NE» 1 ) CALL WIND 

IF (NPUT «EQ*3> THEN ! INPUT WAS RELATIVE VELOCITY 

DO 8 1=1 >3 

8 XD < I ) *XD < I ) +XDA < I ) +XDW < I ) 

END IF 


IF < ABS < NPUT > * EG « 2 ) THEN 
SHDN=CHDN/CLAT 
CHDN*SQRT ( ONE-SHDN*SHDN ) 

IF ( NPUT *LT* ZERO) CHDN—CHDN 
END IF 

REI-RE*SQRT< 1 ♦-GF*SLAT*SLAT) 

IF<R»LT*REI) THEN 
IF(NMET»EQ*2) R»R*1000* 
R-REI+R 
END IF 

X ( 1 ) *R#CLON*CLAT 


DO 10 1*1 >3 

10 XDR< I)*XD< I )-XDA< I )-XDW< I ) 


A(1)*X(1) ! INITIAL INTEGRATION VARIABLES S 

A<4)=X<2) 

A<7)*X<3) 

A<2)=XD< 1 ) 

A(5)=XD<2) 

A<8)*XD<3) 


XS(1)=X<1) 

XS(2)=X(2) 


! SAVE INITIAL POSITION 

//? 


003560 

003570 

003580 

003590 

003600 

003610 

003620 

003630 

003640 

003650 

003660 

003670 

003680 

003690 

003700 

003710 

003720 

003730 

003740 

003750 

003760 

003770 

003780 

003790 

003800 

003810 

003820 

003830 

003840 

003850 

003860 

003870 

003880 

003890 

003900 

003910 

003920 

003930 

003940 

003950 

003960 


003970 

003980 

003990 

004000 

004010 

004020 

004030 

004040 

004050 

004060 

004070 

004080 

004090 

004100 

004110 

004120 

004130 

004140 

004150 

004160 

/\ A A H “7 S\ 


004180 

CALL VCROSS ( XN r XDR > X ) ! ENTRY PLANE NORMAL 004190 

CALL VUNIT(XNfXN) 004200 

004210 

TOTAL WEIGHT = BASIC WGT + ASCENT FUEL + ORBIT FUEL + RESERVE 004220 

WGT=WGTE+WGTF+WGTO+WGTR 004230 

FMASS=WGT/GO 004240 

FMD=ZERO 

FMDF*ZERO 

FMDM=ZERO 


SFC=ZERO 

THRST=ZERO 

XTHRM=ZERO 

RETURN 

END 


004250 

004260 
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TO COMPUTE CL AND CD AT VARIOUS 
ALTITUDES AFTER THE CRAFT HAS SEEN SLOWED TO MACH 0.7 II 


ORIGINAL PAGE IS 
OE POOR QUALITY 


CD 
CD2 = 
CD3 * 
CDO 
CL * 
CFU « 
cri: - 


•. >. or 

CDF £• 


i.1 

,*• *. « »r. 

.. W V 


•• «• l -v r 

LZ-fJ 
CGP.J . 

cdim ' 


COEFFICIENT OF DRAG 
SECOND METHOD FOR FINDING DRAG 

CD2 CORRECTED FOR ERRORS MENTIONED IN ASSUMPTION 13. 

ZERO LIFT DRAG COEFFICIENT 
COEFFICIENT OF LIFT 

FLAT PLATE SKIM FRICTION COEFFICIENT OF THE WING 
TUR SOLENT FLAT PLATE SKIN FRICTION COEFFICIENT OF FHE ROD* 
U - SKIN FRICTTOM drag COEFFICIENT 0T THE WING 
-PRESSURE IK -.C OF THE *uDr 

- F V3E DRAG ZOEFF’ZUN! HhCED ON THE NAXToOM 
FRONTAL . SOD- AREA 

- the ie; fj drag coltt.ii; is. sr of the- sod, 

- *VR ZERO LIFT 5F AG COEFFICIENT OF TJ»E WI:«C l IHC F". 
COfeTFItXEflT OF DRAG DOE IQ i I" 

- •;?SrT!*'ICi i ( ' -F DRAG OF THE tOCv DUE TO ANGLE OF ATT A" .*■ ' 

I.-GUCED DRAG COEFFICIENT CF LING DOG. CC.iOINA ' 1 0- ' 

- JN DRAG COEFFICIENT OF THE NACELLES 

- F RES SURE DRAG C r »EFFICINT DUE TO NACEl LZS 

- INTERFERENCE DRAG COEFFICIENT OF THE NACELLES 


cETA — ' 1 *• M SQUARED 1 
A - ASPECT PATIO 


It* 0.5 


m. 


ANGLE CF ATTACK (IN RADIANS) 


$ - WING SPAN 


CL - COEFICIENT OF LIFT 
3 - WING AREA ~ 

E = OSWALD EFFICIENCY FACTOR 

£2 - 2r,<i METHOD FOR FINDING OSW. EFFIC. FAC. I 

F * EQUIVALENT PARASITE AREA 

SREF * WING PLANFORM OR REFERENCE AREA 

SS * BODY WETTED AREA 

SB * BODY FRONTAL AREA 

SNUET = (JETTED AREA OF THE ENGINE NACELLES 
RHO - AIR DENSITY 

NOTE - RHO DEPENDS ON ALT WHERE M=0.7 IS BEGUN. IN ADDITION 
ITS CHANGE. IE- RHO(l) TO RH0<2>. DEPENDS UPON THE SHIP'S 
RATE OF DESCENT • ! 

V = CRAFT VELOCITY 
DR * DYNAMIC PRESSURE 
L * LIFT ( *WEI6HT ) 

MU * VISCOSITY 

LL » AIRFOIL THICKNESS LOCATION PARAMETER 
RE = REYNOLDS NUMBER BASED ON THE WING CORD 
REB * REYNOLDS NUMBER BASED ON THE BODY LENGTH 
REN = REYNOLDS NUMBER BASED ON THE ENGINE NACELLES 


ASSUMPTIONS 


j < 

4, 


e < 

/ < 

s, 


WEIGHT CHANGE AT THIS POINT IS NEGLIGIBLE. 

THE CRAFT IS FLYING THROUGH AN IDEAL STANDARD ATMOSPHERE 
EACH ITEPATIQN IS A 500 FOOT INCREMENT (30 SECOND INTERS 
SPEED REMAINS AT MACH =0.7 

WAVE DRAG COEFFICIENT IS NEGLIGIBLE AT M=0.7 
THICKNESS RATIO (t/c) IS 0.12 AND IT IS LOCATED 0 0.5C 
BODY LENGTH USED 192.7 FEET 
SHARP LEADING EDGE <i» - R1 LER = 0) 


i*n 


vty 


© 



4 • W4 f tt 4T »(• M* 41 tA 



TAI 


NEGLIGIBLE ON DRAG 
18 NEGLIGIBLE, 
FEET FROM THE BODY 
[NCLUDE MISCELLANEOUS CAU8E8 OF DR AO 
AS CANARDS* INTERFERENCE BETWEEN THE NACELLE8. 
AND SKIN R0U8HNEBS* 



DIMENSION RO ( 50 > * ALT (50) » CCD(50> »VEL(50) »T(50> *MU(50) * 
t RE (SO) *CFW(SO) * CDFW ( 50 > » CDOW ( 50 ) » REB ( 50 ) *CFB(50) * CDFB(50) » 
t CDFB(SO) »CDB2(50) *CD0B(50) »CD0UB(50> » CLAW (50) 

REAL LAM» L» H»MU»LLfLODMAX»LOD 
PI *3. 141592654 
LODMAX-O.O 

BO 1000 J*14*16.2 

AALT*35500 
PP=J 

P*P P-1'2.0 
AL*P«PI/ 180*0 
B AL -AL* l 80 • 0/P I 
WRITE <6*5>BAL 
FORMAT ( ' 1 ’ *T40* ALPHA = '*F5.1*/> . 

WRITE (6*10) 

FORMAT (T5* 'ALT CFT3 ' * T17* ' VEL CFT3 ‘ » T30 * ' DENS CSLUGS -'FT**33 ' * 
S 752* DYN PRESS CLBS/FT**23 ' *T77» 'CD' *T87* 
i CL'*T?6f 'ATM TEMP CR3 ' */) 


ORromr p AGE b 
0f P °0* CILAUTK 


DO 500 1=1*41 
M=0. 7 

BETA*<1 ,0-M**2>**0.5 

KAPPA* 1 .0 

ALT< I ) =AALT-500 • 0 

T( I ) =518. 69- < 0 • 00356) * ( ALT ( I ) ) ' 

MU(I)=1.458E-6*T( I ) **1 . 5/< T< I ) +110. 4) 

R0< I )*0. 0023769* < T ( D/518.69) *45.26669 
VEL(I)*M*(1.4*1617.0*T(I) )**0.5 
RE( I)*RO< I )*VEL< I )*47.3329/MU( I ) 
REN»R0(I)*VEL(I)*30.0/MU(I) 

REB< I ) =R0( I )*VEL ( I ) *192. 7/MU ( I ) 

CFU ( I )*-0.0014*(LOG10(RE( I ) )-9. 0237 ) 

CFB< I ) *-0.001 945* (L0G10( REB ( I)) -8. 419) 

CFN*-0 .00197*L0G10< REN )+0. 01613 

PI-3. 141592654 

A=3. 051 1 

B*154 . 0 

L*880U0 

8-7773. 0+625.0 

SB*PI*( 18* 0**2) 

SS*2*PI*18.0*192. 7 
SP*6937 • 2 
SNWET *6785 ♦ 85 
S*10545. 0+1525. 6 
LL-1.2 

DP-0 . 5*R0 ( I ) * ( VEL ( I ) **2 > *S 
CL-L/DP 

E=t. 0/(1. 0/0. 85+1. 0/1. 4+1. 0/0. 05) 

F*28. 0 

CLAU< I )*(2*PI*A ) /< 2. 0+SQRT < ( A*BETA/KAPPA)**2* 
i (1.0+1. 0/BETA**2) +4. 0) ) 

E2-1 ♦ 1*CLAW( I ) /< A*PI ) . 

CDL-GL**2/(PI#A*E2) 

SREF-7773. 0+625.0 
S-10545. 0+1525.0 ' 



) 



$ 


& 








COFW(I)-CFM( Z >*< 1 . 0+LL*&12 > *2 . 0*7773 . 0/10545 . 0 
CD0W(I)-CDFU(I> 

CDFB(I)-CFB(I)*S8/88 

CDPB ( I > »CFB(I > *< 60 . 0/< 192.7*36 . 0 > **34-0 . 0023* ( 192 . 7/36*0) ) *S8/S8 
CDB2 ( I ) -0 . 029/ < CDFB < I > **0 . 5 ) 

CDOB < I ) -CDFB (i ) +CDPB < I ) +CDB2 ( 1 ) 

CDOUB ( I ) «CDOU ( I ) +CDOB ( I > *SB/SREF 


CDAB»2*AL**2+0.62*1.2*SP/SB*AL**3 
CDIUB«CDL+CDAB*SB/SREF 
CDFN«CFN*SNUET/(PI*36. 0*6.0) 

CDPN=CFN* (60.0/(30. 0/6 . 0 > **3+0 . 0025* ( 30 . 0/6.0 ) ) 
% *SNWET/(PI*36. 0*6.0) 

CDIN^O . 085*6.0 

CH2-CD0WE •' I MCDTWB+CDIN+CDPN+CfcFN 

CD1’~ 1. > 0?*CD2 

CL2 -CL-i-'l- l • *AL 

1-CL2 'CBS . • 

CJ -LTC P2 f AL > / ( 0 ♦ 5*RC ( I ) *VE'L ( I ) **2*SREF > 

CL5/C03 

ALPHA-AL*130.0/PI 

CL4-S0RT ( (CD3-CD0WB( I) )*PI*A*E2> 

IF < AL .LE. 0.0 > CL4=-CL4 


LOD-CL 4/CD3 

IF (LOU .ST. LODMAX) ALF*ALPHA 
IF ( LOD .GT. LODMAX) LODMAX=LOD 
AALT -ALT ( I ) 


UR I TE( 6. 250) ALT ( I > . VEL ( I ) t RO ( I ) . DF> CD3» CL4 . T ( I ) 

250 . F0RMAT<T6.F7.1fT19»F6.2fT34.F10.8.T58fF6.2fT74»F6.3» 

* T85»F6.3fT100»F5.1) V 

300 CONTINUE 
1000 CONTINUE 

WRITE (6. 120) LODMAX »ALF 

120 FORMAT (///»T30» ' (L/D)»a;< » ' . F7.3.T38. ' AT ALPHA * '» 
S F3.2* ' CDEG3 ' »//) 

STOP 

END 
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- 1° 
©20 
30 
40 
50 
60 
70 
3 71 

72 

73 

74 

75 

76 
3 77 

78 

79 

80 
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110 

3120 
130 
135 
140 
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160 
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190 
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220 
3230 
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270 
3 280 
300 
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320 
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390 
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3410 
420 
430 
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REM Chuck Harnden and Mike O’ connor 

REM CALCULATION OF CLALPB AT DIFFERENT MACH NUMBERS 

REM * 3 

REM 

REM CLALPHW = LIFT CURVE SLOPE OF WING 
REM CLALPHC = LIFT CURVE SLOPE OF WING 

REM KNWB, KNBW, AND KNC ARE WING BODY INTERFERENCE FACTORS 

REM SEC s EXPOSED WING AREA OF CANARD 

REM SC s PROJECTED WING AREA OF CANARD 

REM SEW = EXPOSED WING AREA OF WING 

REM SW = PROJECTED WING AREA OF WING 

REM QRAT = AVERAGE DYNAMIC-PRESSURE RATIO ACTING ON THE AFT SURFACE 
REM BC = WING SPAN OF CANARD 
REM BW = WING SPAN OF WING 
REM DC = DIAMETER OF FUSELAGE AT CANARD 
REM DW = DIAMETER OF FUSELAGE AT WING 
REM AW = CROSS SECTIONAL AREA OF FUSELAGE 
CLALPHW = 3.364 

CLALPHC = 3.3315 * 

SEW = 7161 

SW = 10545 

PI = 3.1452 
SEC = 625 

SC = 1849 . 

KWBC = .96 

KWBW = . 99 

KNC = .0845 

QRAT = .9 
IVW = -3 

BC s 86 

BW s 190 

DC = 36 

DW s 36 

AW = 3.42 

INPUT "KBWC=";KBWC 
IF KBWC=0 THEN 430 
INPUT "KBWW=";KBWW 
INPUT **Mach Number =" ;M 
INPUT "BVW" ;BVW 

A= CLALPHC* ( SEC/SC ) *CLALPHW*QRAT*KWBC*IVW* ( BW/2-DW/2 ) 

B = ( 2*PI*AW* ( BVW/2-DC/2 ) ) 

CLALWCV = A/B/57 . 293 

C s ( CLALPHC/ 5 7 . 296 ) * ( KNC+KWBC+KBWC ) *SEC/SC 
D s ( CLALPHW/57 . 296 ) * ( KWBW+KBWW ) *QRAT*SW*SEW/( SC*SW) 

E = CLALWCV 
CLALPH = C+D+E 
PRINT "Mach number =";M 
LPRINT "Mach number = " ; M 
PRINT "CLalpha = CLALPH 
LPRINT "CLalpha ="; CLALPH 
GOTO 260 
END 


3 


M 


J 


rr-.* 


G 



_ Mach number = 1.5 
^ ® CLalpha = .2361887 

Mach number = 2 
CLalpha = .2126911 
Mach number = 2.5 
CLalpha = .206774 
„ Mach number = 3 
> - CLalpha = .202602 


* 3 


} 3 


$ © 


3 


© O 



• 10 REM CHOCK HARNDEN AND MIKE O’CONNOR 

> ® 20 REM ?(LOG<UWH 

30 REM CALCULATION OF DRAG COEF. AT AS A FUNCTON OF ALPHA AND MACH NU1 

40 REM 

50 REM CDO = DRAD COEFFICIENT AT ZERO LIFT 

60 REM CDOWB s ZERO LIFT DRAG COEFFICIENT OF WING BODY COMBINATION 

70 REM CDOP s ZERO LIFT DRAG COEFFICIENT OF THE EXPOSED TAIL PANELS 

) 3 80 REM CDOW = ZERO LIFT DRAG COEFFICIENT OF THE WING 

90 REM CDOB = ZERO LIFT DRAG COEFFICIENT OF THE BODY 
100 REM SB/SREF = THE RATIO OF BODY BASE AREA TO THE REFERENCE AREA 
110 SWETE= 14322 
115 SREF =10545 
120 LAMALE=. 7854 

> 2 130 TOC=.12 

140 SWE=7161 
150 SS=19079.52 
160 SB=1017. 87 
170 INPUT "MACH#'' ; M 
175 IF M=3. 5 THEN 300 
) 5180 INPUT “ CFW" ;CFW 

190 INPUT ‘*CDLE“ ; CDLE 
200 INPUT"CFB" ; CFB 
203 INPUT "Alpha"; ALPHA 
205 CDOB=CFB*SS/SB+.125 
210 CDWW=CDLE+16/3*(TOC)~2*SWE/SREF 
3 215 CDFW=CFW*SWETE/SREF 
220 CDOW=CDFW+CDWW 

230 CDLWB=. 07039+. 6791*ALPHA~2+. 00081*ALPHA A 4 

240 CDOWB=CDOW+CDOB*SB/SREF 

250 CD=CDOWB+CDLWB ' 

260 PRINT "Mach number =;M 
5 3 261 LPRINT "Mach number =;M 

270 PRINT "Alpha ="; ALPHA 

271 LPRINT "Alpha ="; ALPHA 

280 PRINT "CD = ";CD 

281 LPRINT "CD = ";CD 

290 IF ALPHA=1 . 0472 THEN 170 
5 3 293 GOTO 203 

300 END 


3 © 


m © 


l+t 


o © 



Mach number = ;M \><> 

® Alpha = 0 

CD = .3837323 

Mach number = ; M \<* 
Alpha = .17453 
CD = .404419 

Mach number =;M 
® Alpha = .34907 
CD = .4664926 

Mach number = ;M 
Alpha = .5236 
CD s . 5699733 
Mach number = ;M 
© Alpha = .69813 \ 

CD = .7149083 

Mach number =;M 7-<0 
Alpha = .87266 
CD s .9013608 
Mach number = ;M 
©Alpha = 1.0472 
CD a 1.129426 
Mach number =;M 'b-b 
Alpha = 0 
CD a .4728455 
Mach number =;M T&r%S 
© Alpha = .17453 
CD = .4935321 

Mach number =;M ' L,it 

Alpha = .34907 
CD s .5556057 f 

Mach number =;M L ' 

Q Alpha = .5236 
CD s .6590864 
Mach number = ; M 
Alpha = .69813 
CD = .8040214 

Mach number = ;M 
G Alpha s .87266 
CD = .9904739 

Mach number = ; M 
Alpha s 1.0472 
CD s 1.21854 
Mach number =;M 
G Alpha = 0 
CD a . 5218952 
Mach number =;M ? » 
Alpha = .34907 
CD = .6046555 

Mach number = ;M , 

C Alpha = . 5236 
CD = .7081361 

Mach number =;M /t ' 
Alpha = .69813 
CD = .8530711 

Mach number =;M 
© Alpha = 1.0472 
CD = 1.267589 

Mach number = ; M 3* 
Alpha = .87266 
CD = 1.039524 

Mach number =;M l ->” 
©Alpha = 1.0472 
CD = 1.267589 

Mach number = : M 


** 
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ci>‘ = .55123 

_ Mach number =;M I f 
P • Alpha = . 17453 

CD = .5719166 ^ 

Mach number =;M <*• 
Alpha = .34907 
CD s .6339903 
Mach number = ;M 
) 3 Alpha = . 5236 

CD = . 7374709 

Mach number = ;Mi 
Alpha s .69813 
CD = .8824059 

Mach number = ;M 2-* 

3 Alpha = .87266 
CD = 1.068858 

Mach number = ;M ? 
Alpha = 1.0472 
CD = 1 . 296924 
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PROGRAM FOR FLYBACK BOOSTER LIQUID ROCKET ANALYSIS 


PI=3. 14159 
XISP®320. 0 
TB00S®1500000.0 

TTOT2=6230 ♦ 0 ' 

GAMMA*! . 24 
XMBAR=22 » 0 
GC=32. 2 

GE=32 . 2 V 

RLITT=2. 3 
PA-1 *04 - 
RBAR=1545 ♦ 0 

DCTMB® ( GC*TBOOS ) / < X I SP*GE ) 

R=RBAR/XMBAR 

CP=(GAMMA*R>/(GAMMA-1.0) 

CV=CP/GAMMA 

UE=XISP*GE 

PE=PA 

PT0T2=PE/ ( ( 1 . 0- < ( UE**2 > / ( 2 . 0*GC*CP*TT0T2 )))**( GAMMA/ 

(GAMMA-1.0))) 

XME=SQRT ( (2.0/(GAMMA-l .0) )#( ( (PT0T2/PE) **( (GAMMA-1 .0)/GAMMA) )-l .0) ) 
TE=TT0T2-( (UE**2)/(2.*CP*GC) ) 

A=SQRT(GAMMA*GC*R*TE) 

TSTAR=TT0T2/( 1 ♦ + ( (GAMMA-1. )/2« ) ) 

ATHOA=SQRT ( GAMMA#GC#R*TST AR ) 

USTAR=ATHOA 

QR=(UE**2)/(2.*GC*(1.-(TE/TT0T2) ) ) 

QRBAR=QR*XMBAR 

GAMM1=SQRT(GAMMA*(2.0/(GAMMA+1.0) >**( ( GAMMA+1 ♦ 0 )/ (GAMMA-1 . 0) ) ) 

WRITE ( A f 79? ) TBOOS t XMBAR f RLITT 

FORMAT ( 'O' f5Xf ' T hrust* ' fF9. 1 fIXf ' lbf ' f5X» ' M= ' f F4 ♦ 1 f 17X f ' r =' * 

F3.1) 

WRITE ( A f 800 > R f CP f CV 

FORMAT ('O' f5Xf 'R®' fFA.3f 17Xf 'Cj*®' fF7.3f 13Xf 'Cv®' fF7.3) 

WRITE ( 6 f 801 ) XISP f UEf XME 

FORMAT ( ' 0 ' f 5X * ' Isf® * t F5 ♦ 1 f IX f * sec * 1 12X f * Ue = ' f F9 » 3f IX f ' ft/sec * f 
4XF / Me= / FF5.3) . 

WRI TE ( 6 f 802 ) DOTMB f PA f PE 

FORMAT ( '0' f5Xf 'hi®' fF9.3f1Xf ' lbm/sec ' f AX f ' P s® ' fF5.2f1Xf 'psi'f 
IIXf 'P e®' fF4.1f1Xf 'F si') 

WRITE(Af803)PT0T2f TT0T2 f TE 

FORMAT ( 'O'fSXf 'P t=' fF9.3f1Xf 'F> si' f9Xf 'T t®' fFA.IfIXf 'R'f 
12Xf 'T e® ' fF8.3f IXf 'R' ) 

WRITE ( A t 804 ) TSTAR f USTAR f QR 

FORMAT ( 'O' fSXf 'T ='fF8«3f1Xf'R'f12Xf'U ®' f F8. 3f IXf ' ft/sec ' f 
5Xf 'Q r=' fF12.3) 

DOTM=DOTMB 





ASTAR=<D0TM/GAMM1*SQRT< <R*TT0T2)/GC> >/<PT0T2#144;0> 

ASTAR=1.55 

PT0T2=<D0TM/GAMM1*SQRT< < R*TT0T2 > /GC > >/< ASTAR*144> 

DSTAR=SQRT ( ASTAR*4 . O/PI ) 

GAMM2= < GAMMA+1 ♦ 0 )/< 2 ♦ 0* ( GAMMA-1 ♦ 0 ) ) 

AEXIT=ASTAR/XME*( <2./<GAMMA+l . ) >*(l. + < (GAMMA-1, >/2, ) *XME**2) )**GAMM 
DEXIT=SGRT(AEXIT*4.0/PI> 

N=5 

WRITE < 6 f 809 ) N y 

809 FORMAT ( 'O' r5X> 'NumbeY of engines '»I2) 

WRITE ( 6 » 81 0 ) ASTAR » DSTAR 

810 FORMAT < 7X » 'A * ' > F7 . 3> 1X-* ' so ft '»5Xf'D *' »F6*3» 1X» 'ft ' ) 

WRITE(6r811 ) AEXIT fDEXIT ' 

811 FORMAT ( 7X » ' Ae= ' r F7 ♦ 3 > IX > ' so ft ' >5X> 'De=' >F6,3r 1X» 'ft ' > 

STOP 

END / 
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ABSTRACT 


Work was completed this quarter on the design of the 
Shuttle II orbiter. Included in this -final report is a brie-f 
explanation and continuation o-f the preliminary and 
intermediate design phases. During the -final phase o-f the 
project, the trajectory analysis was -finalized. Also, all 
A.S.T.S. groups agreed upon respective masses with regard to 
the center o-f gravity locations and vehicle mating. Overall 
c.g. locations -for assorted weight conditions were determined 
for Shuttle II. Nonintegral versus integral internal tank 
structures, one-piece wing planform, and differing internal 
structural makeup were investigated. An overall cost analysis 
was also performed to complete the Shuttle II project. 
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INTRODUCTION 


The information contained within this report deals with a 
design effort which was carried out on a new, improved Space 
Shuttle Orbiter, known as Shuttle II. This project was 
undertaken as part of a design study consisting of 
preliminary, intermediate, and final design stages. These 

stages are now reviewed. 

I. Review of Preliminary Design Stage. 

Design was begun fall quarter 1986. Initial design 
specifications set forth by NASA for the Shuttle II called 
f or s 

* Two-stage vehicle utilizing a fully reusable, winged, 
flyback booster. 

t Payload capable of 40,000 lb. transported to space 
station orbit (28.5 degree inclination at 270 NM altitude). 

* 40,000 lb. payload return capab ility. 

* Cargo bay having dimensions of 15 ft. diameter X 60 
ft. length. 

* Engine propellant (second stage) of liquid oxygen and 
liquid hydrogen. 

* Maximum velocity of 7000 fps. for staging. 

* Maximum acceleration of 3 G’s during ascent into orbit. 

* Projected first flight in 2005; therefore, technology 
and design freeze approximately 1997. 

Design work was coordinated with the Booster group. 

Using a gravity turn analysis, it was determined that the 
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Shuttle II should have an initial mass of one million pounds. 
A Shuttle II vehicle weight was theh calculated using a mass 
ratio o-f 4.5. 

It was subsequently -found that the calculated -final 
velocity was not great enough to achieve the specified 
in-orbit altitude. Because o-f this, thrust and velocity were 
varied with values -for velocity and mass ratio being tabulated 
in an unsuccessful attempt to achieve the correct final 
velocity. 

By taking into account the determined initial mass of one 
million pounds, the Shuttle II was tentatively dimensioned 
with increased wing and fuselage areas. A computer program 
was written which utilized the assumed mass ratio of 4.5 to 
calculate the mass of the Shuttle II, mass of the fuel, and 
the corresponding fuel volume. 


II. Review of Intermediate Design Phase. 

Work on the Shuttle II design continued into the 
intermediate phase winter quarter 1987. Trajectory was 
finalized for the ASTS by use of the program S2TRAJ.BAS. 
Through this program, it was found that Shuttle II achieved 
its desired trajectory and orbital requirements. Briefly 
investigated also was a secondary trajectory into orbit 
utilizing a technique similar to a Hohmann transfer. 
(Appendices A,E) 

Also written were computer programs which calculated 
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masses and volumes of the propellants as well as possible 
lengths and diameters of the propellant tanks. Corresponding 
empty and -full tank masses were also calculated. (Appendices 
B,C> 

A propulsion system -for the Shuttle II was selected based 
on such criteria as: number o-f engines required, engine c.g. 
locations, and required -firewall thicknesses. An Orbital 
Maneuvering System <OMS) and a Reaction Control System (RCS) 
were also selected. 

Finally, structural design was undertaken to determine 
st *ch considerations as: selection o-f construction materials 
with respect to vehicle structural loads, material cost 
comparisons, and overall Shuttle II con-figuration. 

III. Review o-f Final Design Phase 

Work on the -final design phase began spring quarter 1987. 
All ASTS groups agreed upon respective masses with regard to 
launch by use o-f the program S2TRAJ . BAS. Weight and center o-f 
gravity locations were calculated by use of the program 
ACCEPT. FOR. A cost analysis on production of the Shuttle II 
was also performed. Nonintegral versus integral internal fuel 
tank structures were investigated with a decision being made 
to use nonintegral structures. Alternatives were considered 
such as different placement of cargo bay and fuel tanks to 
facilitate a more advantageous c.g. location. 
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THEORY 


A. Trajectory 

Initial design specifications call -for the Shuttle II to 
achieve a circular orbit at a 270 NM altitude. To reach this 
orbit, a direct powered ascent was chosen. Important 
parameters of the trajectory such as time of flight, velocity, 
range, altitude, thrust, and G-forces were continuously 
calculated at time intervals from launch until attainment of 
desired orbit. These values were calculated using an 
iterative-type program, S2TRAJ . BAS. S2TRAJ executed these 
calculations with input of initial mass of the Shuttle II, 
engine thrust, and the altitude at which the gravity turn was 
initiated. (Fig. 1A, IB) 

A secondary method was considered which used a powered 
ascent up to some elliptical orbit with a perigee near 60 
nautical miles and an apogee of 270 NM. Upon reaching this 
apogee, an impulse is applied using the OMS to adjust the 
elliptical orbit into a circular orbit. The maximum velocity 
change required for this maneuver would be 376 fps. To derive 
at this value, necessary velocities for orbits with perigees 
from 55 NM to 110 NM and apogees from 220 NM to 270 NM were 
calculated. This value was compiled by the program ORB shown 
in Appendix E. As can be seen from the section in this report 
concerning propulsion, this velocity change is within the 
capability of the OMS. It is hoped that this method will 
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reduce the fuel required to attain orbit and in turn, the 
total mass required for the shuttle system. 

The program mentioned, which determined the actual 
trajectory, can also be used to determine the elliptical orbit 
for this secondary method. In order to determine if the 
burnout data supplied by this program actually attains orbit, 
the program ASTRO can be used (see Appendix F) . This program 
takes burnout data and calculates the orbital elements. The 
burnout data can either be in the 3-dimensional celestial 
system mentioned previously, or in the 2-dimensional system 
used by the trajectory program. Pertinent results of 
importance for this program are as follows: 

H P <Nf1? «<NM) e Vave(fos) Vdiff 


355. 

185. 

3710. 

.0229 

24723 

285 

350. 

190. 

3710. 

.0216 

24739 

268 

345. 

195. 

3710. 

.0202 

24755 

252 

340. 

200. 

3710. 

.0189 

24772 

235 


TABLE 1 


S 



Data -found in this table conforms to the necessary orbital 
requirements as set by NASA. These necessary requirements set 
by NASA -for the circular orbit phase are as follows: 



a 

3710 nautical miles 



e 

0 (circular orbit) 


Jl 

, Omega 

261 degrees (longitude of 

ascending 



node) 



i 

28.5 degrees (latitude of Cape 



Kennedy) 


• 

, Theta 

0 to 360 degrees (circular 

orbit) 

w 

, omega 

0 to 360 degrees (circular 

orbit) 


However, because of the uncertainty involved in this procedure 
and the deduced assumptions, it was decided to proceed with 
the direct powered ascent trajectory. It should also be 
stressed that a Hohmann transfer be considered and anaLyzed 
for the Shuttle II’ s final orbit achievement. 


6 



B. Propulsion 

From trajectory analysis, it was found that Shuttle II 
will be required to supply one million pounds of thrust at 
liftoff. The engine chosen for Shuttle II was the STME 481 
rocket engine. (Fig. 2) A minimum of three of these engines 
will be necessary to produce the required liftoff thrust. The 
Shuttle II will be equipped with four of these engines to 
provide an engine-out successful mission capability, and also 
to increase the lifespan and reliability of the engines since 
they will not be required to operate at their maximum 
capability. 

The examination of the efficiency of the Shuttle's 
engines becomes a pertinent matter in mission planning and 
cost analysis. Each engine is capable of producing 396,492 
pounds of thrust. The flight configuration with four engines, 
•nginfs, and two engines (an impossible configuration as 
shown below), have been tabulated as follows: 
lumber gf gnqjnqs Thrust E ach Produces V. of Wax Thrust 


4 

230,000.33 

63 

3 

333,333. 33 

84 

2 

500,000.00 

• 126 


TABLE 2. 



As can be seen by the -four and the three engine cases, the 
percent of the maximum thrust necessary to complete the 
mission are 63 percent and 84 percent, respectively. The 
percent difference between these two cases is 25 percent. The 
overall advantage yielded by the three engine versus four 
engine case is that more missions can be performed per engine 
in the four engine case. This can be seen in that the 25 
percent difference between the operating thrusts of the two 
cases yields that for every four three-engine missions, five 
four engine missions can be performed. Another important 
difference, as mentioned above, is the safety factor. The 
engine— out possibility for the three-engine mission would be 
fatal to the mission as can be seen from the two-engine data 

in Table 2. However, the engine— out possibility for the 

four— engine case would allow, with 16 percent reserve thrust, 
the success of the mission. Therefore, with the discussed 
advantages and the comparison of data in Table 2, it is 
intuitively obvious that the four — engine case be implemented. 

The engine configuration decided upon utilizes a 
rectangular engine placement. (Fig. 3) The dimensions of this 
rectangle are approximately 35 ft. in the vertical . and 32.5 
ft. in the horizontal direction. 

The maneuverability of the Shuttle II will be assisted by 
this rectangular engine configuration together with engine 
gimbaling. The engines will be gimbaled 10.5 degrees for 
pitch control and 8.5 degrees for yaw control. The limiting 
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area that the end of each nozzle can traverse is indicated by 
an ellipse with eccentricity o-f 0*3843 with reference to the 
nozzle as the center of the ellipse. 

Investigated also was an Orbital Maneuvering System (QMS) 
for the Shuttle II. The OMS for the current orbiter consists 
of two 6000 pounds force engines. This system carries enough 
fuel to produce a total velocity change of 1000 fps. A normal 
shuttle mission requires a velocity change of less than 400 
fps. f allowing 600 fps. as needed for any abnormal 
circumstances. On performing a mass analysis for Shuttle II 
it was decided to keep the OMS at two 6000 lbf. engines with 
an increased fuel capacity. It should be noted that with the 
increased mass of Shuttle II maneuver time will be increased. 
(Fig. 4) 

The Reaction Control System (RCS) is utilized for 
altitude correction of Shuttle II. After investigation, it 
recommended that the existing RCS in use on the present 
orbiter be utilized for Shuttle 1 1 as well. 
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C. Aerodynamics 

The aerodynamic considerations -for the Shuttle II are 
based upon its reentry requirements. These requirements 
dictate that Shuttle II be able to fly in the speed range of 
Mach 25 at reentry to Mach 0.28 at touchdown. This range 
demands hypersonic, supersonic, transonic, and subsonic flight 
configurations. 

The hypersonic flight regime, with a speed range of Mach 
5 and higher, is the least explored regime of the four. The 
accumulated data retrieved from the vehicles that have flown 
in this regime is limited. 

The supersonic regime (with Mach numbers from 5 to 1.2), 
the transonic regime (with Mach numbers from 1.2 to 0.8), and 
the subsonic regime (with Mach numbers below 0.8) are all 
speed regimes which have been thoroughly explored for many 
years by many flight vehicles. Because of this, data 
concerning flight configuration and parameters is easily 
obtained. However, it should be noted that the shape of 
Shuttle II will be heavily dependent on the characteristics of 
the reentry (hypersonic) regime. 

In consideration of the effects that the .hypersonic 
regime had on Shuttle II’s aerodynamics, it is important to 
note that the key point concentrated upon in this regime was 
overall shape. (Fig. 5) A breakdown of Shuttle II’s key 
components and why they were shaped and sized as they were is 
listed below: 
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* Nose Cap — Contoured -for hypersonic trim, performance, 
and heating 

* Fuselage - Sized by payload requirements 

* Double Delta Wing - Sized by 188 kt. landing design 

speed 

* Vertical Tail - Sized by subsonic stability 

* Flared Rudder — Speed Brake — Rudder sized by crosswind 
landing 

* OMS Pod - Sized by OMS Required fuel tankage 

* Aft Fuselage - Sized by Shuttle II main engines 

< Body Flap — Sized to protect main engines from reentry 
heating 

* Full Span Elevons/Ai lerons — Sized by hypersonic trim; 
pitchdown maneuver 

It is important to note that the shape of the Shuttle II 
wing was determined through reference to previous NASA studies 
on shuttle concepts. (Ref. 2) The planform shape decided 
upon was that of an initial sweepback angle of 79 degrees with 
an outboard sweepback of 63 degrees. (Fig. 6) 

It is also important to note that although the main 
component shapes were a function of the hypersonic reentry 
phase, that the subsonic and hypersonic regimes dictate the 
stability and control analysis of Shuttle II. 

An important component of the stability and control 
analysis is the correct determi nati on of overall center of 
gravity location. A breakdown of all major Shuttle II 
components by individual c.g. location is necessary in order 
that the overall c.g. location may be determined. 
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To analyze the major components, a total weight breakdown 
of all major structural sections of the existing orbiter (ie. 
wing group, tail group, etc.) was obtained from the Marshall 
Space Flight Center in Huntsville, Alabama. (Appendix G) 

By considering only the major components, an estimated 
weight analysis was determined by a ratio method. This method 
included comparison of the existing orbiter mass-to-area 
ratios to Shuttle II mass (unknown) -to-area (known) ratios. 
From this comparison, approximate component and section masses 
were found. (Fig. 7) 

A similar method was utilized in the calculation of the 
approximate c.g. locations with reference to an origin located 
at the Shuttle II nose. By utilizing geometric configuration 
comparisons, c.g. locations relative to the nose are 
represented with a rectangular coordinate system, XYZ. 

Two extreme conditions were considered for calculation of 
the c.g. locations. The first was that of launch with a 
40,000 pound payload, full fuel complement, and all provisions 
and support equipment. The other condition was that of 
landing with no cargo, all fuel expended, and minimum amount 
of provisions and support equipment to contribute to weight. 
(Figs. 8,9) 

Also considered was a third condition of landing with a 
40,000 pound payload, all fuel expended, and minimum 
provisions and support equipment. (Fig. 10) Figure 11 
illustrates the reference coordinate system and the resulting 
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c.g. shifts taking into account th» above three conditions. 

In conclusion, it is important to realize that although 
the Shuttle II design was based upon all four flight regimes, 
each regime in itself was a contributing factor to the success 
of the overall design conf i guration. The hypersonic regime 
dictated the shape and size, while the sub-hypersonic regimes 
dictated the location and size of components which in turn 
affect the stability and control of Shuttle II. 



D. Materials and Structures 

In selecting the materials -for the air-frame o-f the 
Shuttle II, the conventional tradeo-f-f studies considering the 
weight and cost o-f production had to accurately re-flect the 
compatibility requirements demanded o-f the structure by the 
external Thermal Protection System (TPS). The very low 
strength, brittle TPS tiles, when bonded to Shuttle II’s skin, 
demand that when exposed to 115 percent o-f maximum expected 
loads during ascent and 100 percent during atmospheric 
reentry, that the skin not buckle. This skin-buckling 
requirement make? aluminum and titanium equivalent in terms o-f 
weight criteria since both materials have a ratio o-f 
compression modulus to density approximately equal to 10<10E7) 
inches. Since -flight loads can accumulate high stresses. 
Shuttle II’ s structure favored the high-strength material such 
as titanium. Another consideration for selection of the 
structural material involved the thermal effects at 
atmospheric exit and entry at elevated temperatures. After 
considering such factors as TPS compatibility, heat 
capacitance, and strength as set forth by the buckling 
requirements, titanium proved to be the best rated material 
for the primary structure in terms of limiting total weight. 
However, the cost of titanium is approx i mat el y 300 percent 
greater than that of aluminum. Therefore, in terms of 
economic practicality, the selection of aluminum tuas .*•« as 
the primary structural material for Shuttle II resulting from 
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these cost end production risk comparisons. 

A major area of the structure where aluminum is not the 
primary structural material is the main engine thrust 
structure. The thrust structure supports one million pounds 
o-f main engine thrust and distributes the vertical stabilizer 
loads onto the -fuselage. Here, the possibility o-f truss 
compression/buckling was much too great for the aluminum 
structural supports. 

From Reference (15), NASA studies in minimum weight 
structural concepts were observed. These studies identified 
tradeoff allowances for the orbiter structure. Some of these 
studies also proved beneficial to Shuttle II’s structure. For 
Shuttle II, a space frame concept for the thrust structure as 
an alternative to the competing plate girder concept together 
with an aft carrythrough spar with the bulkhead instead of a 
floating carrythrough was found to alleviate much unneeded 
weight. Through these weight and cost studies, NASA engineers 
selected the composite material systems used in the payload 
bay doors, the OMS pod external shell, and thrust structure. 
However, the graphi te-epoxy composite used in the payload bay 
doors was found to fail due to moisture being absorbed by 
diffusion, accompanied by degradation at temperatures .i.ove 
250 degrees Fahrenheit. For Shuttle II’s structure, it was 
decided to us aluminum-lithium honeycomb panels in place of 
the graphite-epoxy composite. Figure 12 displays the primary 
structure with the corresponding materials used. 
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C. ' 3 


Also considered was the interrelation of the structure 
and aerodynamic characteristics. Emphasis was placed on such 
concerns as the size and mass of each fuel tank and its 
effects on the structural weight. Another concern mss the 
required increase in wing area and thus an increase in total 
material and weight of Shuttle -II. 

Also investigated were the advantages and disadvantages 
of nonintegral as compared to integral internal fuel tank 
structures. The results found are as follows: 

I. Nonintegral Internal Fuel Tanks 

Advantages 

1. Allows expansion space due to changes in temperature. 

2. Reduced tank masses due to thinner shells made 
possible by insulation of tanks with reliable, light-weight 
material. 

3. Allows space for conventional airplane-type 
spar/str inger structural conf iguration beneath tanks for 
support of wing/fuselage junction. 


Di sad vantages 

1. Propellant tanks develop exterior moisture. This 
necessitates adequate plumbing and drainage vents to remove 
moisture from tank section. 

2. Required truss support systems result in increased 
structural mass. 


3. Enlarged fuselage diameter also results in increased 
structural mass. 
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II. Integral Internal Fuel Tanks 


Advantages 

1. Less structural materials required -for truss support. 

2. Reduced -fuselage diameter results in less structural 
mass o-f -fuselage. 

3. No moisture removal system will be necessary because 
moisture will -form on exterior o-f -fuselage. 


Di sadvantages 

1. No available space for spar /stringer structural 
supports of wing/fuselage Junction. 

2. Increased mass of tanks due to corresponding increase 
in shell thickness of each tank necessitated by thermal 
protection requirements. 

3. Additional fuselage structural supports necessary to 
accomodate the increased tank masses. 

By comparing the advantages and disadvantages associated 
with these two concepts, it was determined that the 
nonintegral internal fuel tanks proved to be the most 
advantageous even though a moisture removal system is 
necessary. 

After aerodynamics had determined the best location 
within the fuselage to balance the center of gravity, the fuel 
tanks had to be positioned so as not to alter this location. 
The tanks will be made of 2219 Aluminum alloy with a skin 
thickness of 0.03166 ft. The four fuel tanks give a total 
empty mass of 32,301 pounds. The diameter and length of the 
three oxygen tanks are 15 ft. and 80 ft., respectively; while 
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those for the hydrogen tank are 30 ft. and 80 ft. (Figs. 


13, 14) 

Using these 

dimensions 

as design 

cri teria, 

it 

was 

found that Shuttle 

II 

will 

definitely 

increase 

over 

the 

present 

orb iter in 

not 

only 

material masses but 

also 

in 


structural masses. 

Accompanying this mass increase due to placement of fuel 
tanks within the fuselage, it was necessary to increase the 
wing area in order to supply the required lift. Again, this 
increase in structure demands additional material, and thus, 
more weight. 

Additionally considered was the wing/fuselage junction. 
It was decided to utilize a conventional airplane-type of 
design of spars and stringers continuously connected to form a 
single wing planform which the fuselage would then be set 
upon. (Figs. 15,16) 

It should be noted that at this point the Shuttle II 
structure consists of a significant amount of unused internal 
volume -for which a suitable purpose has not yet been found. 
With the advancement pf new technology and the discovery of 
new composites, the additional mass resulting from this unused 
volume may be diminished, thereby enabling more, favorable 
performance character i sti cs. 
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Di scussi on/Summary 


Final design -for the Shuttle II Mas primarily concerned 
with -finalization of trajectory, agreement between all ASTS 
groups on respective masses, weight and c.g. analysis, and 
weight and cost considerations with respect to materials and 
structures. 

The use of the program S2TRAJ enabled the finalization of 
the Shuttle II trajectory. Knowledge of this precise 
trajectory into orbit provided important values such ass 
initial mass <M) , mass at separation from booster. Shuttle II 
liftoff thrust, orbital inclination angles ($, 2P ) , and final 
velocity and altitude. Use of S2TRAJ by the other ASTS groups 
(Booster and Cargo Vehicle) enabled agreement on such issues 
as duration and magnitude of thrust during launch and correct 
vehicle mating locations with respect to weights and c.g. 
locations. The values obtained through the trajectory 
analysis dictated important design considerations such ass 
total amount of fuel required, overall configuration size, and 
required structural materials. 

Research into the desired propulsion unit for Shuttle II 
had to satisfy certain specifications, among them being the 
capability of successful completion of mission after one 
engine loss. Therefore, a total of four main engines for 
Shuttle II was decided upon. Slight redesign to the Orbital 
Maneuvering System in the form of increased fuel capacity was 
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decided upon -for maneuvering of the enlarged configuration. 
This decision for added fuel rather than an additional engine 
is also cost effective in terms of final weight. One 
disadvantage to this choice is the increased time required for 
each maneuver. 

The design phase concerning aerodynamics for Shuttle II 
was primarily concerned with the reentry portion of the 
mission profile. The main consideration was the shaping and 
sizing of the major Shuttle II components with respect to the 
particular speed regime through which the vehicle is 
traveling. Another major consideration was that of stability 
and control analysis. This analysis depends heavily upon 
correct determination of the overall vehicle c.g. location 
which was found using the program ACCEPT. FOR. 

The main guideline followed in the selection of materials 
for Shuttle II was a comparison to the existing orbiter and 
modern material technology. Taken into consideration were 
compatibility requirements, thermal effects, and . costs. 
Comparisons between aluminum and titanium were investigated, 
with aluminum being chosen as the primary structural material. 
Selection of composite materials was not made in favor of 
aluminum honeycomb panels for certain areas of Shuttle II such 
as: payload bay doors and the OMS external pod shell. Also 
investigated were nonintegral versus integral internal fuel 
tank structures with the decision being made to utilize 
nonintegral tanks. A cost analysis was performed for Shuttle 
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II which indicated that an approximate total cost of 
development and production would be 6.558 billion dollars. 
(Fig. 17) 

Finally, it should be noted that this design proposal 
(Fig. 18) is a major change over the existing orbiter. A 
change in not only dimensions, (Figs. 19,20), but also in 
concept to a degree. The Shuttle II will hopefully prove to 
be a much more adaptable, multi-role oriented component of the 
Advanced Space Transportation System than the present orbiter 
was ever envisioned to be. In conclusion, it is hoped to one 
day see the Shuttle II and Booster (Fig. 21) lifting off 
together bound for another successful mission in space. 


20 



LIST OF REFERENCES 


r?« 

Vk 

fei 



1.) 


2.) 

iP 

f - ■ i 

£ 




k 

3.) 

0 

4*) 



1 


r; 

5.) 

IS 

6.) 

I 

7.) 

m 


l.H 

Vji 

8*) 

P 

m 

9.) 

h 

10.) 

Ki 

11.) 

f. 0 

LL 


par* 

1-3 

f ■■■> 

12.) 

rfrr 


c— 

V'V; 

i 

13.) 

Uy 

r 

14.) 

Q: 


tr. 

15.) 


111 


jv* 


Bat* Roger. R.* *t. al. ? 'Fundamentals of Astrodvnamics' • 
Oov*r Publications* Inc.* N*w York* 1963* 

Compton* Harold R.* *t» al. ? 'Stability and Control Over 
th* Su** r sonic and Hypersonic Speed Rani» # « Shuttle Per- 
formance* Lessons Learned* NASA Conference Publication 
2283 r Pt* 1* 1983* pp 473-509. 

Gardner* James P.* et* al* ? 'Applied Astronautics* An 
Introduction to Space Flight'. Prentice Hall Inc** En- 
glewood Cliffs* N»J«* 1963* 

Glynn* Phillip C.*l Hoser* Thomas L.? 'Orbiter Structural 
Design and Verification' * Space Shuttle Technical Confer- 
•nett NASA Confvrtnct Publication 2342>Pt« If 1983* 

Hill* P*G*> Peterson* C.R*? 'Mechanics and Thermodynamics 
of Propulsion'. Addison-Weslew Publishing Company* 1945* 

Joels* K.M* I Kennedy* G.P.? 'The Space Shuttle Operator's 
Manual'* Ballantin* Books* New York* 1982* 

Kaplan* Marshall H.? 'Space Shuttle? America's Wings to th* 
Future'* Aero Publishers* Inc»* 1978* 

Merrell* Grawson* Capt* USN (ret.) * 'Handbook of Satellites 
and Space Vehicles'* D. von Nostrand Company* Inc.* New 
York* 1965* 

Row* A.E. ? 'Orbital Motion'. Halsted Press* 1982. 

Rupee* Harry 0*? 'Introduction to Astronautics'. Academic 
Press* New York and London* 1966. 

Stuhlinger* Ernst.? 'Space Science and Engineering'. McGraw- 
Hill Book Company* 1965. 

Sutton* G.P . $ Ross* D.M.? 'Rocket Propulsion Elements'. John 
Wiley ft Sons* New York and London* 1976. 

Tawlor* John U.R.* *d«? 'Jane's All th* World's Aircraft* 
1985-86'* Jane's Publishing Company* Limited* London* 

1986* pp 497-498. 

Thompson* W.T.? 'Introduction to Space Dynamics'. John Uilew 
ft Sons* New York and London* 1961. 

Underwood* Jimmy M. > Young* James C.? 'The Aerodynamic Chal- 
lenges of th* Design and Development of th* Space Shuttle 
Orbiter'. Space Shuttle Technical Conference* NASA Confer- 
ence Publication 2342* Pt. 1* 1983* pp 209-264. 


§ 


21 



1 . 

^ is*) 'Columbia Excaads Flight Goal* 7 . Aviation Waak * Sraca Tach 
_ nolodv. April 20» 1981 t pp 18-33. 

m 

17 .) 'Orbital Fliaht Handbook. Voluaa I 7 . Martin Coaranw. 1963. 


m 

y 


fl 


m 


i 


0 


T*? ft 

L? 



r:~ 

i 


f:3 

p-f 

*- ‘ ' 




iv< 

uw 


m 

i 



FIGURE X A» GRAPHICAL REPRESENTATION OF TRAJECTORY 
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FIGURE It: TABULATED DATA 


PROGRAM 1J 



INITIAL MASS 

7*400*000.00 

LBM. 

(INPUT) 

G-TURN ALTITUDE 

2*625.00 

FT. 


SHUTTLE II LIFTOFF THRUST 

1*000*000.00 

LBF ♦ 

(RESULTS) 

SHUTTLE MASS AT SEPARATION 

1 * 290*040*00 

LBM. 


FINAL ALTITUDE 

1*643*340.00 

FT. 



270.00 

NM 


FINAL VELOCITY 

23*668.00 

FT./S. 


PSI 

103.80 

DEGREES 


GAMMA 2 

103.80 

DEGREES 


SHUTTLE MASS IN ORBIT 

241*898.00 

LBM. 

PROGRAM 21 





SHUTTLE MASS IN ORBIT 

241*898.00 

LBM. 

(INPUT) 

MASS RATIO 

4.59 



MASS OF SHUTTLE FUEL BURNED 

1*008*142.00 

LBM. 


MASS OF OXYGEN 

864*121.90 

LBM. 

(RESULTS) 

MASS OF HYDROGEN 

144*020.11 

LBM. 


MASS OF PAYLOAD 

40*000.00 

LBM. 


VOLUME OF HYDROGEN TANK 

32*510.23 

FT. ”3 


VOLUME OF OXYGEN TANKS (3) 

12*099.16 

FT. "3 

PROGRAM 3: 





VOLUME OF HYDROGEN / ONE TANK 

32*510.23 

FT. ”3 

(INPUT) 

MAX* DIAMETER OF TANK 

30.00 

FT. 


MAX* LENGTH OF TANK 

80.00 

FT. 

(RESULTS) 

DIAMETER OF TANK 

25.07 

FT. 


LENGTH OF TANK 

74.42 

FT. 


VOLUME OF OXYGEN (EACH OF THREE 

TANKS) 4*033.05 

FT. "3 

(INPUT) 

MAX. DIAMETER OF TANK 

15.00 

FT. 


MAX* LENGTH OF TANK 

80.00 

FT. 

(RESULTS) 

DIAMETER OF TANK 

9.11 

FT. 


LENGTH OF TANK 

65.41 

FT. 

PROGRAM At 


* 



WALL THICKNESS 

0.03166 

FT. 


METAL DENSITY 

176.256 

LBM/FT~3 

(INPUT) 

VOLUME OF HYDROGEN 

32*510.230 

FT. "3 


INSIDE TANK DIAMETER 

25.035 

FT. 


MASS OF HYDROGEN 

144*020.000 

LBM. 

(RESULT) 

MASS OF EMPTY TANK 

16*338.000 

LBM. 


VOLUME OF OXYGEN PER TANK 

4*033.050 

FT. ”3 

(INPUT) 

INSIDE TANK DIAMETER 

9.075 

FT. 


MASS OF OXYGEN PER TANK- 

288*040.000 

LBM. 

(RESULT) 

MASS OF EMPTY TANKS 

5*3210.00 

LBM. 
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FIGURE 3 SQUARE POSITIONING OF ENGINES 


FIGURE 4 


OhS MASS ANALYSIS 


SHUTTLE I SHUTTLE II 


MASS SHUTTLE 

155*000 

LBM 

241*898 

LBM 

MASS PAYLOAD 

65*000 

LBM 

40*000 

LBM 

TOTAL MASS 

220*000 

LBM 

281*898 

LBM 

VELOCITY CHANGE 

1*000 

FT/S 

1*000 

FT/S 

MASS/VELOCITY 

220 

LBM S/FT 

282 

LBM S/FT 

SCALE FACTOR 

220/220 

= 1 

282/220 

= 1.28 

PROPELLANT 

N20H 

14*866 

LBM 

19*028 

LBM 

MMH 

9*010 

LBM 

11*353 

LBM 

TOTAL 

23*876 

LBM 

30*561 

LBM 


CHANGE IN PROPELLANT MASS 
6*685 LBM 



o 


I 





e: 

r 

r 


i 


i 

i 


i 

i 



[s 


i 




c. 


■ 

! - 


O 

ce 

u. 

UJ 

x 

CO 

CO 

h 

u 

i UJ CD 
I I- z 
i o mi 
' Ct h- 

u. <r 

? Ui 

C X 
r- 

>- 

Cl Cl £C 
C Ul 
U N 2 
Ji. Ml UJ 

V 

« * 

o 


>-4 UJ 

ocao 
uj o<t 
> 0 .X 

X z 

UI ^ < 
2 C 0 K 

< X 

x o> 

~ A 
-J 

< xo 

H UJ UJ 
whN 
A CO h 
K > CO 
O CO 

» 


CO 


r 


Ui 


u: 

CD 

2 

Mi 

Ni 



a 

t 

mi 

m 

X 

UJ X 

c 

2 mi 

2 

mi CL 

V 

X 1 - 


U 

X 

•> CD 

i£ 

H- U 2 


UiH N 


O 2 »- 


co o <r 


1 CO UJ 


o^x 


Mi Ui 

m 

hCifi 


C >- 2 

Oi 

dCX 4 

im 



vice » 

T 

CO O UJ 

**■» 

UJ u. u 


2 2 


uj a «r 


2 UJ X 


Ml « cc 


li.DC 


o u. 


un- ct 


h 2 UJ 


<CO Cl 


ii:a 


UJ 


a * 


o 


X 




27 


FIGURE 7 


$ 

WEIGHT ANALYSIS OF MAJOR SHUTTLE COMPONENTS 
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FIGURE 8 

X CASE ♦ !> 

MAXIMUM LOAD 



fjENTER OF GRAVITY LOCATION 

r~ 

CHART WITH DESCRIPTION OF 

COMPONENT 

AND ITS MASS 

u 

INSCRIPTION OF COMPONENT 

X-COORD 

Y-COORD 

Z-COORD 

MASS 

( CARGO 


61*00 

0.00 

10.00 
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pXYGEN TANK # 1 

W/FUEL 
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-9.50 

2.00 
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OXYGEN TANK * 2 
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FIGURE 9 
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18.00 

7.00 

2.00 

0.00 

(PERSONNEL 

i * 

18.00 

0.00 

9.00 

2873.00 

RESIDUAL FLUIDS 

rr 

18.00 

-6.00 

7.00 

15227.00 

Ucs PROPELLANT 

165.00 

8.00 

5.00 

0.00 

rOMS PROPELLANT 

| 

168.00 

-6.00 

5.00 

0.00 

T5 

i 

V^E COORDINATES WITH RESPECT TO 
OF THE SHUTTLE OVERALL C*G* ARE 

i( C*G* » 89.25 

V Y C.G. - -0*45 

(Z C.G* - 2.57 

THE SHUTTLE 

AS follows: 

NOSE. 




•TOTAL MASS » 249946*00 

t ; : 'ii 




ORIGINAL PAGE IS 
of poor QUALrna 



FIGURE 10 J CASE *3f MINIMUM LOAD 




fcj 

i CENTER OF GRAVITY LOCATION 

CHART WITH 

DESCRIPTION OF 

COMPONENT 

AND ITS MASS 

II 

DESCRIPTION OF COMPONENT 

X-COORD 

Y-COORD 

Z-COORD 

MASS 

J .: 1 

''•CARGO 

61*00 

0.00 

10.00 

0.00 

rpXYGEN TANK # 1 WO/FUEL 

59*00 

-9.50 

2.00 

5210.00 

i:r 

OXYGEN TANK # 2 WO/FUEL 

59.00 

0.00 

3.00 

5210.00 

(OXYGEN TANK * 3 WO/FUEL 

t.'.;* . 

59.00 

9.50 

2.00 

5210.00 

HYDROGEN TANK WO/FUEL 

130*00 

0.00 

5.00 

16338.00 

t 3TME ENGINE * 1 

172.00 

-8.00 

- 2.00 

7140.00 

pTME ENGINE * 2 

172.00 

8.00 

- 2.00 

7140.00 

5TME ENGINE * 3 

174.00 

8.00 

4.00 

7140.00 

f.STME ENGINE * 4 

174.00 

-8.00 

4.00 

7140.00 

OMS ENGINE * 1 

f..' --T1ir 

173.00 

13.00 

4.00 

1050.00 

{•OMS ENGINE * 2 

173.00 

-13.00 

4.00 

1050.00 

{ .REACTION CONTROL SYSTEM 

K- „ , r 

165.00 

0.00 

9.00 

1404.00 

WING GROUP . 

131.00 

0.00 

-12.00 

15950.00 

[BASIC STRUCTURE 

94.00 

0.00 

-5.00 

37455.00 

TAIL GROUP 

169.00 

0.00 

34.00 

7828.00 

[• THRUST STRUCTURE 

165.50 

0.00 

6.00 

3874.00 

f BODY FLAP 

183.00 

0.00 

-10.00 

939.00 

v induced ENVIRONMENT 

18.00 

0.00 

5.00 

30230.00 

[FRONT LANDING GEAR 

21.00 

0.00 

-5.00 

2245.00 

RIGHT LANDING GEAR 

131.00 

24.00 

-10.00 

2600.00 

auil .... ||B am n -n-T- ui -- - - . 

jfwEFT LANDING GEAR 

131.00 

-24.00 

-10.00 

2600.00 

,-fUEL CELLS 

22.00 

8.00 

-5.00 

276.00 

' REACT ANT DEWARS 

21.00 

-8.00 

-5.00 

498.00 

vir, 
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V-V: 


|<E ACT ANTS 

23*00 

1 

1 

01 

♦ 

o 

o 

-6.00 

655.00 

f PATTERIES 

24*00 

4*00 

-5.00 

675.00 

ililLEC. CONV. AND DIST. 

125*00 

-1.00 

-8.00 

8476.00 

SURFACE CONTROLS 

F; 

172.00 

3*00 

-10.00 

3932.00 

t 

AVIONICS 

15*00 

0*00 

5.00 

4956.00 

Environmental controls 

18*00 

-5.00 

10.00 

4625.00 

PERSONNEL PROVISIONS 

/ • ?: 

18*00 

7.00 

2.00 

0.00 

Personnel 

18.00 

0.00 

9.00 

2873.00 

r 'RESIDUAL FLUIDS 

18*00 

-6.00 

7.00 

15227.00 

1 RCS PROPELLANT 

165*00 

8.00 

5.00 

0.00 

jlMS propellant 

168.00 

-6.00 

5.00 

0.00 

s 

THE COORDINATES WITH RESPECT TO 
PC’ THE SHUTTLE OVERALL C.G. ARE 

C»G« » 94*63 

Y C*G. * -0.53 

£Z C.G. » 1*15 

‘TOTAL MASS ■ 209946.00 

THE SHUTTLE 
AS FOLLOUS: 

NOSE* 




[5 

I 

p" . 

I 

B 

i., 

I 
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FIGURE 11 t CENTER OF GRAVITIES UNDER SPECIFIC 

LOADING CONDITIONS IN THE 
XYZ- COORDINATE SYSTEM 
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CONVENTIONAL ALUMINUM 2219 STRUCTURE 
MAXIMUM TEMPERATURE 350F 
ALUMINUM-LITHIUM ALLOY SECONDARY STRUCTURE 
PROTECTED BY REUSABLE SURFACE INSULATION 
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Figure 16 SIDE CUTAWAY VIEW OF STRUCTURES 



FIGURE 17 J 

LAUNCH VEHICLE COST NOOEL 
STAGE (HE 
COST ESTIMATE 
OATEt 

INFLATION t 
UT RESEM: 


23-HV-87 ! 

1.229 FHH 82$ TO 84$ j 


1.23 



CER'S 


COMPLEXITY 

COST 



1N0EP. MAR 

• 

FACTOR 

(84 $ MILLIONS) 

COST ELEMENT 

DOTI£ 

TFU 

DOTirE 

TFU 

DDT&E 

TFU 



■ ■— 

— 

•••• 


" M T ^ 

STRUCTURES/TPS 

73297.10 

73297.80 

1.00 

1.00 

450.84 

194.43 

THERMAL CONTROL 







EWIRGHENTAL CONTROL 

4423.01 

4423.00 

1.08 

1.00 

28.33 

2.18 

BASE HEAT SHIELD 

939.88 

939.00 

1.00 

1.00 

9.93 

1.30 

WING/TAIL/LEAD IN6 EDGE 

23743.00 

23743.00 

1.00 

1.00 

348.14 

44.93 

LANDING GEAR 

7443.88 

7443.00 

1.80 

1.00 

54.32 

14.11 

AVIONICS 

4934.88 

4934.88 

1.80 

1.00 

914.74 

84.84 

ELECTRICAL POUER 

9131.80 

9151.08 

1.00 

1.00 

408.14 

99.47 

PROPULSION (LESS EN61NES) 

274.00 

274.08 

1.00 

1.00 

37.48 

1.88 

SEPARATION PROVISIONS 

1133.00 

1133.00 

1.00 

1.00 

29.97 

9.84 

SURFACE FLIGHT CONTROLS 

3932.88 

3932.08 

1.00 

1.00 

231.84 

78.82 

AUXILARY POUER UNIT 

1.88 

8.88 

1.00 

1.00 

0.88 

0.00 

HYDRAULICS 

1484.88 

1404.88 

1.00 

1.00 

34.47 

18.72 

SUBTOTAL 

124134.00 

124134.08 



2428.33 

445.39 

STRUC. TOOLING 

194.43 . 


1.00 


424.84 


SYS. TST. HROU. 4 ASSS18LT 

443.39 


1.00 


483.48 


SYS. TST. OPS. 

483.48 


1.00 


222.37 


SUBTOTAL 





3488.42 

443.39 

6SE 

3488.42 


1.00 


343.03 


SUBTOTAL 





4243.44 

443.39 

SEfcl 

4243.44 

443.39 

1.00 

1.00 

410.14 

37.31 

SUBTOTAL 





.4433.81 

302.70 

PR06RAH MANAGEMENT 

4433.81 

302.78 

1.00 

1.00 

137.73 

17.84 

SUBTOTAL 





4813.34 

320.33 

EN61NES (OU*Tv*ISP*Pc) 

391234493.30 

1343018782.00 

1.00 

1.00 

42.94 

0.04 


SUBTOTAL 
FEE (ICO 

PROG. SUPPORT (37. DB/., 27. PROD.) 
SUBTOTAL 


4834.30 320.39 


479.91 

144.09 


833.37 

4357.87 


72.88 

11.87 


5702.30 403.34 


90.80 

494.14 


COST CONTINGENCY (IE) 
TOTAL 
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Figure 18 3-D VIEW OF SHUTTLE II 
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APPENDIX A 

PROGRAM 1 WITH FLOWCHART 
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PROGRAM 1 


■ - 5 REM 

•10 REM 

15 ALT-0 

,„20 DT-O* 1 

| > 30 RANG-0 

■ 40 ISP1-320 

50 ISP2-380 

f,60 INPUT "ENTER INITIAL MASS"*M1 

[65 M-Ml 

70 INPUT "ENTER ALTITUDE FOR G-TURN" > AGT1 

,75 AGT-AGT1 

K 80 TH-1.35*M 

1 90 INPUT "ENTER LIFTOFF THRUST (SHUTTLE II ENGINES) " rTH21 

95 TH2-TH21 

[•100 TH1-TH-TH2 
Pi 10 MD0T1-TH1/ISP1 
120 MD0T2-TH2/ISP2 
p 130 ISP— TH/(MD0T1+MD0T2) 

[ 140 MDOT-TH/ISP 

150 AMC1-0 
.... 160 MC2-0 

[170 Y-0 

1 ^ 180 X-0 

190 GO-32.2 

[ 200 RO— 2 • 092E+07 

p 210 VO-O 

220 V-VO 
P 230 T-0 

K240 PRINT 

v 250 PRINT"TIME V RANG ALT G'S MASS 

260 PRINT 

k 270 F*-"##*.* *****.* ****** ****** *♦** ******* 

k 280 M2»M-MD0T1*DT-MD0T2*DT 
290 G-GO* < RO/ < R0+ ALT ) ) "2 

f 300 DV«ISP*GO*LOG < M/M2 ) -G*DT 

P 310 AG1-DV/G0/DT 

320 IF AG1>3 THEN 400 
f 330 CHY-0 . 5* ( VO+VO+DV ) *DT 

p 340 IF Y+CHY>AGT THEN 580 

1 350 Y-Y+CHY 

360 ALT-Y 

IP 370 V-VO+DV 

U 380 AG-AG1 

390 GOTO 460 

f 400 DV-3*G0 

U 410 M2*M/EXP< <DV+G*DT) /ISP/GO) 

420 A-DU/GO 

430 V-VO+DV 

p, 440 Y-Y+O ♦ 5* ( VO+V ) *DT 

m 450 ALT-Y 
460 T-T+DT 

p 470 DM-M-M2 

m 480 MDT-DM/DT 
490 TH-MDT*ISP 

p 500 MC1-MC1+TH1/ISP*DT 
p 510 MC2»MC2+(TH-TH1 )/ISP2*DT 
l ' ; 520 MC-MC1+MC2 

I 


THRUST* 

******** 
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'530 PRINT USING F*»T» V*RANG» ALTr AG>M»TH 
540 REN 
550 M-M2 
560 VO-V 
570 GOTO 280 
! 580 PRINT 

:590 PRINT’TIME V RANG ALT G'S MASS THRUST* 

600 PRINT* " 

r- 620 M=M+MD0T1*DT+MD0T2*DT 
! 630 DPSI«(PI/180>*0.1 

-640 PSI0-PI/120 

650 PSI-PSIO 
(660 GAM-1 

l. 670 Z0«SIN<PSI0>/<1+C0S(PSI0>> 

680 PRINT*TIME V RANG ALT PSI GAM THRUST MASS G'S* 

f 700 Ft-*##* ***** ####### ####### ##♦# #*♦# ####### ####### *♦##* 

| 710 PRINT 
720 N-TH/M 

. 730 PRINT USING F*!T» V?RANG» ALT »PSI*180/PI rGAMfTHf Mr AG 

f 740 C-V0/<Z0~<N-1))/(1+Z0~2) 

l - 750 PSI-PSI+DPSI 

760 Z-SIN<PSI)/(1+C0S<PSI> ) 

(£ 770 V1«C*<Z~<N-1> )*< 1+Z~2> 

b 780 IF Vl>7000 THEN 1390 
785 V-Vl 

« 790 DT- C/G*Z~< N-l ) * < 1/ < N-l ) +Z~2/ < N+l > > -C/G*ZO~ < N-l > * ( 1/ < N-l ) +Z 0~2/ ( N+l ) ) 

800 DX-0*5*<V0*SIN<PSI0)+V*SIN<PSI)>*DT 
810 DY-O* 5* ( VO*COS < PSIO ) +V*CQS < PSI ) ) *DT 

820 AG«<V-VO>/DT/GO 
pi 830 IF AG>3 THEN 1060 
1 * 840 X-X+DX 

845 Y-Y+DY 

{ 850 THETA-ATN < X/ < R0+ Y ) ) 
l 860 GAMMA-PSI -THETA 

870 GAM-GAMMA* 180/P I 
{ 880 ALT- <Y+RO> /COS < THETA )-R0 
900 RANG-R0*THET A 

930 T-T+DT 
, 940 PSIO-PSI 

r 950 ZO-Z 
t 960 VO-V 

970 IF ALT< 100000 THEN 1020 
I 980 ISP1-340 

990 ISP2-462 

1000 MD0T1-TH1/ISP1 
r 1010 MD0T2-TH2/ISP2 

1020 M»M-MD0T1*DT-MD0T2*DT 
L - 1030 IF M<0 THEN 1390 
1 040 G-GO* < RO/ ( RO+ALT ) ) "2 
1050 GOTO 720 

{'■■ 1060 N»3*G0/G+< < l-Z0"2)/< 1+Z0~2> ) 

1070 C*VO/ZO~ < N-l ) / < 1 +Z0~2 ) 
p: 1080 V-C*Z"<N-1)*<1+Z~2> 

« 1090 D1-C/G*Z~< N-l ) * < 1/ < N-l ) +Z~2/ < N+l ) ) -C/G*ZO~ ( N-l ) * < 1/ ( N-l ) +Z0~2/ < N+l ) ) 

1100 DX-0.5*(V0*SIN<PSI0)+V*SIN(PSI>>*DT 
1110 DY»0.5*(V0*C0S(PSI0)+V*C0S(PSI))*DT 

1120 X-X+DX 
L 1130 RANG“RO*THET A 

fS 

I... 
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•iv“ 

' 1140 
1150 
1160 
1170 
1180 
i 1190 
'1200 
1210 
'1220 
i 1230 
1240 
r 1250 

I 1260 

'1270 
1280 
(1290 
ti 1300 
1310 
| 1320 
[1330 
1350 
r 1360 
| 1370 
> 1380 
1390 
1:1391 
l, ,1392 
1400 
,1420 
L 1421 
"1422 
r 1423 
F (1424 
L 1430 
1431 
[1432 
( : 1435 
1440 
1450 
1460 
1470 
1480 
j 1490 
t" 1500 
1510 
: 1520 
1530 
1540 
1550 
1560 
‘ 1570 
1580 
(V- 1590 
[>1600 
1610 
r 1620 
1630 
u 1640 

I 


Y-Y+DY 

ALT* < Y+RO ) /COS < THETA ) -RO 

T-T+DT 

PSIO-PSI 

ZO-Z 

VO»V 

TH-N*M 

TH1-TH-TH2 

HD0T1-TH1/ISP1 

MD0T2*TH2/ISP2 

M»M-MD0T1*DT-MD0T2*DT 

MC1-MC1+MDQT1*DT 

MC2«MC2+MD0T2*DT 

MC-MC1+MC2 

IF M<0 THEN 1390 

6*60* < RO/ < R0+ ALT ) ) ~2 

PSI*PSI+DPSI 

THETA-ATN < X/ < RO+Y ) ) 

6AMMA-PSI-THETA 

6AM*6AMMA*180/PI 

PRINT USIN6 F$*T>V»RAN6*ALT»PSI»6AM?THrM»A6 
IF V>7000 THEN 1390 
Z-SIN<PSI)/<1+C0S<PSI> ) 

60T0 1060 
PRINT 

PSI-PSI-DPSI 


PRINT VI 

PRINT'TIME V RAN6 ALT PSI 6AM 

PRINT* — 

PRINT 

PRINT* SEPARATION 

PRINT 

PRINT* 

PRINT*ti»« v ran* *lt Mi *«»2 

PRINT 

F$«**** ***** ******* ******* ***** ***** 

M-<1-.15*M1/M>*M 


THRUST MASS 6'S" 


thrust »*ss si's 

******* ******* *.** 


■ 


ISP-ISP2 


TH-TH2 

MDOT-TH/ISP 

VO-V 

PSIO-PSI 

6AM2-PI/2 + THETA 
6AMM«6AM2*180*/PI 
ZO*SIN(PSIO)/< 1+C0S<PSI0) ) 

N-TH/M 

PRINT USIN6 F*JTfVfRAN6f ALTfPSI*180./PIf6AMM»TH»M»A6 
C* VO/ < Z0~ ( N- 1 ) ) / < 1 +Z0"2 ) 

PSI-PSI+DPSI 
THETA-ATN < X/ ( RO+Y ) ) 


6AM2-PI/2+THET A 
6AMM-6AM2* 1 80 ♦ /P I 
IF PSI>6AM2 THEN 2080 
Z*SIN<PSI )/( 1+C0S<PSI ) ) 

V«C* ( Z~ < N-l ) ) * < 1+Z~2 ) 

IF V>25000 THEN 2080 

DT-C/6*Z~< N-l ) * < 1/ < N-l ) +Z'“2/ < N+l > > -C/6*Z0~ < N-l > * < 1/ < N-l > +Z0~2/ ( N+l ) ) 
DX-.5*<VO*SIN<PSIO)+V*SIN<PSI)>*DT 
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1650 
1660 
1680 
1700 
,1701 
1702 
1710 
1720 
' 1730 
1740 
1750 
r 1760 
1770 
1780 
1790 
i 1800 
! 1810 
1820 
! 1830 
i 1840 
1850 

r I860 
[1880 
1 1881 
1882 
1890 
1900 
1910 
(1920 
1930 
‘1940 
.1950 
1960 
t 1970 
1980 
1990 
2000 
2010 
,2020 
2025 
1 2040 
2045 
2050 
: 2060 
2070 

f 2080 
[ 2081 
2090 
r 2100 
■ 2110 
ho 21 11 
2112 
$2113 
tS2114 
2115 

p 2116 

E2117 


m 

m 


DY- ♦ 5* ( VO*COS < PSI 0 ) + V*CQS < PSI ) > *DT 

AG«<V-VO)/DT/GO 

X-X+DX 

Y-Y+DY 

THETA-ATN < X/ < RO+Y ) ) 

RANG-RO*THETA 

ALT- < Y+RO ) /COS < THETA ) -RO 

T-T+DT 

PSIO-PSI 

ZO-Z 

VO-V 


M*M-MDOT*DT 
IF H<0 THEN 2080 
G-GO* < RO/ < R0+ ALT ) ) ~2 


GOTO 1520 

N-3*G0/G+< ( 1-Z0*'2) /< 1+Z0''2) ) 

C-VO/ZO" < N-l ) / < 1 +Z0~2 ) 

V*CtZ^( N— 1 )$( 1+3^2 ) 

DT-C/G*Z~ ( N-l ) * < 1/ < N-l ) +Z~2/ < N+l ) ) -C/G*ZO" < N-l ) * < 1/ < N-l ) +Z0~2/ < N+l > 


DX«.5*<V0*SIN<PSI0)+V*SIN<PSI))*DT 
DY- . 5* < VO*COS < PSIO ) +V*COS < PSI ) ) *DT 
X-X+DX 
Y-Y+DY 

THETA-ATN < X/ < RO+Y ) ) 


> 


RANG-RO*THETA 

ALT- < Y+RO ) /COS < THETA ) -RO 


T-T+DT 

PSIO-PSI 

ZO-Z 

UO-V 

TH-N*M 

MDOT-TH/ISP 

M«M-MDOT*DT 

IF M<0 THEN 2080 

G-GO* ( RO/ ( R0+ ALT ) ) "2 

AG-3 

PSI-PSI+DPSI 
THETA-ATN < X/ < RO+Y ) ) 

GAM2-PI/2+THETA 

GAMM-GAM2*180*/PI 

PRINT USING F$»T i> V » RANG » ALT * PS I *180* /PI *GAMM»TH»M> AG 
IF PSI>GAM2 THEN 2080 
IF V>24000 THEN 2080 
Z»SIN(PSI >/( 1+C0S<PSI ) ) 

GOTO 1800 

REM 

PRINT 

PRINT* INITIAL MASS J*>M1 

PRINT* ALT FOR G-TURN J*?AGT1 

PRINT ‘THRUST < SHUTTLE II ) J*JTH21 
PRINT "FINAL VELOCITY **»V 
PRINT 'FINAL MASS t*JM 
PRINT*FINAL PSI t * *PSI#180/PI 
PRINT'FINAL ALTITUDE J**ALT/6080 
PRINT 'FINAL GAMM ♦ * rGAMM 
INPUT "TRY AGAIN Y-l N-2*?TRY 
IF TRY<2 THEN 15 
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FLOWCHART FOR TRAJECTORY ANALYSIS 


PROGRAM •S2TRAJ* 


Q START ) 


ENTER INITIAL MASSf SHUTTLE II ENGINE; 
THRUST f AND GRAVITY TURN INIATION 
ALTITUDE 


!- 


<D 


CALCULATE TIME OF FLIGHT f 
VELOCITY f RANGE f ALTITUDE r 
G-FORCEf MASS* AND THRUST. 


T 


PRINT TIME OF FLIGHT f VELOCITY 
RANGE f. ALTITUDE f G-FORCE f MASS 
AND THRUST. 


< 



IS THE G-FORCE GREATER 
THAN THE ALLOWED SAFETY FACTOR 
OF THREE? 



YES 



■*© 


NO 


DOES THE 

ALTITUDE EQUAL THE 
SET G-TURN INITIATION 
ALTITUDE? 



m. 


YES 


CALCULATE TIME OF FLIGHT f VELOCITY f 
FLIGHT PATH ANGLE FROM VERTICLEf THRUST f 
MASSf AND G-FORCE USING FORMULAS FOR 
GRAVITY ASSISTED TURN MANEUVER 


© 


© 
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YES 



IS THE 
HASS LESS THAN 
ZERO ? 



YES 


3 

NO 

l 

SUBTRACT OFF MASS OF BOOSTER 
STAGING OCCURS 




CALCULATE TIME OF FLIGHT * VELOCITY * 
RANGE* ALTITUDE* FLIGHT PATH ANGLE* 
THRUST* MASS* AND G-FORCE. 





© 
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PRINT FINAL TIME» VELOCITY* RA 
FLIGHT PATH ANGLE t THRUST t 
G-FORCE 
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APPENDIX B 

PROGRAM 2 WITH FLOWCHART 


56 



nnnn non o o ooooooo 


PROGRAM 2 


THIS PROGRAM CALCULATES PERTINENT MASSES AND VOLUMES 
OF CERTAIN COMPONENTS OF THE SPACE SHUTTLE II PROGRAM* 


REAL MR f MI f MSH f MC f M02H2 * MH2 f M02 
MC-40000* 

200 PRINT *f 'INPUT INITIAL MASS AND MASS RATIO* ' 
READ *fMIfMR 

MSH-<MI/MR)-MC 
M02H2- < MC+MSH ) * < MR-1 . > 

THE MASS RATIO OF 0/H IS 6.0 

MH2-M02H2/7.0 

M02-6.0*MH2 

THE DENSITY OF 0 IS 71.42 LBM/FT 
THE DENSITY OF H IS 4*43 LBM/FT 


VH2-MH2/4.43 

V02-M02/71.42 

VT-VH2+V02 


100 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
READ *fA 
IF<A«EQ 
PRINT * 
PRINT * 
PRINT * 
GOTO 200 
STOP 
END 


'MASS INITIALLY 
'MASS OF SHUTTLE 
'MASS OF CARGO 
'MASS OF 02 AND H2 
'MASS RATIO 
'MASS OF H2 

'MASS OF 02 

/ / 

'VOLUME OF 02 AND H2 
'VOLUME OF H2 
'VOLUME OF 02 


MI 

MSH 

MC 

M02H2 

MR 

MH2 

M02 

V02H2 

VH2 

V02 


If 'POUNDS MASS 
;f 'POUNDS MASS' 


' f MI f 'POUNDS MASS' 

' fMSHi 
'fMCf'I 
'fM02H2f 'POUNDS MASS' 
' F MR F 'POUNDS MASS' 
'fMH2f 'POUNDS MASS' 
f'fM02f 'POUNDS MASS' 

'fVTf 'CUBIC FEET' 
'fVH2f 'CUBIC FEET' 
'fV02f 'CUBIC FEET' 


'WANT TO CALCULATE NEW VALUES?' 
' 1-YES 2-NO' 

2) GOTO 100 
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FLOW CHART FOR PROGRAM 2r MR. FOR 


i- ; ! 


w 


I 

b 


(begin J 

1 

INPUT 

INITIAL MASS AND MASS RATIO. 


CALCULATE 
MASS OF SHUTTLE 
TOTAL MASS OF 02 AND H2 
MASS OF 02 
MASS OF H2 
VOLUME OF H2 
VOLUME OF 02 
TOTAL VOLUME 


OUTPUT 

INITIAL MASS 
MASS OF SHUTTLE 
MASS OF CARGO 
TOTAL MASS OF 02 AND H2 
MASS RATIO 
MASS OF H2 
MASS OF 02 

TOTAL VOLUME OF H2 AND 02 
VOLUME OF H2 
VOLUME OF 02 


INPUT 
DECISION TO CONTINU 
A 



-N 




rr 

■ 


I..-. 
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APPENDIX C 

PROGRAM 3 WITH FLOWCHART 


K:' 

r. 

& 





r.- : 

i- - 


iv 
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PROGRAM 3 


THIS PROGRAM TAKES A GIVEN VOLUME OF FLUID t ALONG WITH VARIABLES 
FOR MAXIMUM DIAMETER AND LENGTH » AND FINDS DIMENSIONS FOR ALL 
POSSIBLE TANK SIZES THAT COULD ACCOMIDATE THE VOLUME OF FUEL* 


REAL L 

OPEN(UNIT»l>FILE='TUOLOX.DAT' rSTATUS«'OLD' ) 

PI-3.141592654 

101 WRITE(6f 102) 

102 FORMAT </////) 

PRINT *» 'INPUT VOLUME ( FT • "3 )» MAX .DIAMETER (FT. ) r MAX .LENGTH (FT.)' 

READ *fV»DMAXfLMAX 

URITE(6f103> 

103 FORMAT ( ' VOLUME RADIUS DIAMETER LENGTH' ) 

WRITE(6» 104) 

104 FORMAT ( ' > 

R-DMAX/2. 

99 X-<V-<4./3.)*PI*R**3)/(PI*R**2> 

- L«X+2*R 
D«2*R 

URITE(6r 105) VrRrDrL 

105 FORMAT (1X»F9.2»F11*4»F11*4»F12.4) 

URITE(1f*)L»D 

p a ^ 99$R 

IF ( L » GE * LM AX ) GOTO 100 
GOTO 99 
100 PRINT *»' 

PRINT Mf / / 

PRINT *»' TRY AGAIN??' 

PRINT *f' YES-1 NO-2 ' 

READ *fN 

IF(N»EQ*1)G0T0 101 
CLOSE UNIT-1 
STOP 
END 
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FLOW CHART FOR VOL. FOR 



L ^AD N / 



YES 
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APPENDIX D 

PROGRAM 4 WITH FLOWCHART 


£ 

t-: 

W 

( - 

I 

6 

I 

t 

I 

r* 

i'i: 

I 

1 

I 
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PROGRAM 4 


THIS PROGRAM CALCULATES THE WEIGHT OF THE FUEL TANKS GIVEN THE W^LL 
THICKNESS* THE DENSITY OF THE METAL* THE VOLUME OF FLUID*THE INSIDE 
DIAMETER* AND THE FUEL MASS* 


REAL LO * MASSTANK * MTANKFUEL * MLIQUID 

10 PRINT ft* 'ENTER WALL THICKNESS (FT*)* AND METAL DENSITY(LBM/FT~3>?' 
READ «*WTH*RHO 

PRINT ** 'ENTER VOLUME OF FLUID (FT. ~3 ) * AND INSIDE DIAMETER (FT ♦ )?' 
READ «*VI*DI 

PRINT #* 'ENTER FUEL MASS ?' 

READ **MLIQUID 
PI-3. 141592654 
D0-DI+2.*WTH 
RI-DI/2* 

X«(VI-(4*/3*)*PI#(RI)**3)/(PI*RI*#2> 

L0«X+2.*WTH 

RO-DQ/2* 

V0-(4./3. >*PI:M*0**3+X*PI*R0**2 
VMETAL-VO-VI 
MASSTANK-RHO*VMETAL 
MTANKFUEL-MASST ANK+MLIQUID 


WRITE(6*100) 

100 FORMAT (' TANK AND FUEL MASS TANK MASS FUEL MASS') 
WRITE(6*99) 

99 FORMAT ( ' '> 

WRITE ( 6 * 98 ) MTANKFUEL * MASSTANK * MLI QU ID 
98 F0RMAT(F13.2*12X*F8.2*4X*F9.2) 

WRITE(6*97) 

97 FORMAT (//////) 

PRINT ** ' ' 

PRINT **' ' 

PRINT **' TRY AGAIN ? YES-1' 

READ **NN 

IF ( NN ♦ EQ ♦ 1 ) GOTO 10 

STOP 

END 
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FLOW CHART FOR METAL. FOR 







APPENDIX E 

PROGRAM 5 WITH FLOWCHART 
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ORBITAL VELOCITY PROGRAM 


THIS PROGRAM IS DESIGNED TO CALCULATE THE VELOCITIES 
REQUIRED FOR VARIOUS ELLIPTICAL ORBITS 

REAL KM 
R=3440*0 
HA=220 . 

HP=55* 

NNN=1 

f KM=62746* 8 
WRITE<5r 001 ) 

FORMAT(//f 9X> 'Ha' >3 X» 'Hp' r5x? 'a' t 7xt 's' t5xr 'Vp' t5x? 'Vs' r5x» 
S' Vcs' ' Vdiff ' ) 

DO 105 K*1 r 12 

RP=R+HP 

RA=R+HA 

f • DO 100 I»l»ll 
A=<RP+RA)/2* 

E=( A~RP)/A 

r VP=SQRT <2*#KM#( 1 */RP-l »/<2*#A> > )#6080. 

| ; VA=SQRT ( 2 ♦ *KM* < 1 ♦ /RA-1 . / ( 2 ♦ *A ) ) ) #6080 ♦ 

VCA=SQRT ( KM/RA ) #6080 ♦ 

VDIFF*VCA-VA 
fu NNN=NNN+1 

m WRITE ( 5 f 220 > HP » HA » A » E f VP » VA » VCA > VDIFF 

220 FORMAT </f5Xr2F6*l»F8»lfF7«4» 3F8 ♦ 0 » F7 ♦ 0 ) 
rr HA-HA+5* 
iA RAaRA+5. 

IF ( NNN ♦ GT ♦ 28 ) THEN 
URITE(5f 221 > 

FORMAT (//////) 

NNN*0 

END IF 

JO CONTINUE 
HP=HP+5* 

HA=220 . 

105 CONTINUE 
l END 


fc; >; 
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FLOW CHART FOR ORB ♦ FOR 
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APPENDIX F 

PROGRAM 6 WITH FLOWCHART 



I 


I 

- 

l>* 
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ASTRO* FOR 


THIS PROGRAM IS DESIGNED TO CALCULATE THE 
ORBITAL ELEMENTS GIVEN THE BURNOUT DATA 
OR THE BURNOUT VECTORS GIVEN THE ORBITAL 
ELEMENTS 


MAIN PROGRAM 

CONTROLS THE TYPE OF DATA TO BE OBTAINED 

REAL IfKfM 

COMMON/DATA/X r Y t Z r XDOT t YDOT f ZDOT f A t ECCEN f THET A t OMEGA 1 I 1 U 
50 WRITE<6f 100) 

iOO FORMAT <5Xf 'IF YOU WISH TO FIND BURNOUT DATA TYPE 1 'f/f5Xf 
B' IF YOU WISH TO FIND THE ORBITAL ELEMENTS TYPE 2'f/f5Xf 
O' INPUT DATA MUST BE IN UNITS OF RADf NM AND NM/SEC') 
READ<5f*)J 
TF<J.EQ.1)THEN 
CALL ORBIT 
ELSE 

CALL BURN 
END IF 
WRITE(6f 110) 

LiO FORMAT <//f5Xf 'IF YOU UISH TO CONTINUE TYPE 1') 

READ<Sf*)K 

P IF(K*EQ*1)G0 TO 050 
U END 

SUBROUTINE ORBIT 

!J SUBROUTINE ORBIT 

CONTROL INPUT OF ORBITAL ELEMENTS AND 
OUTPUT OF BURNOUT DATA 

B 

REAL I 

p COMMON/DATA/X r Y t Z f XDOT f YDOT t ZDOT t A f ECCEN f THET A f OMEGA f I f U 

WRITE <Af 300) 

$00 F0RMAT</f5Xf' INPUT OF ORBITAL ELEMENTS' f/f5Xf 'A*' > 

. READ < 5 f-¥) A 

t- URITE(6 f305) 

Lj5 F0RMAT</f5Xf 'ECCENTRICITY-') 

READ<5f«) ECCEN 
f WRITE<6 f310) 

I 10 F0RMAT</f5Xf' THETA-') 

READ<5f«) THETA 
i : WRITE<6f315> 

f . 5 FORMAT < / f 5X f ' OMEGA- ' ) 

READ<5f«) OMEGA 
_ WRITE<6f320) 
pO FORMAT </f5Xf'I-') 

K READ(5f*)I 

WRITE <6? 325) 

FORMAT < / f 5X f ' W— ' ) 

& READ(5f*)W 


!'-*• 
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CALL BURNOUT 

WRITE<6f 330) Xf XDOT * Y f YDOT f Zf ZDOT 
'10 FORMAT <//f9Xf 'RADIAL COMPONENTS' f5Xf 'VELOCITY COMPONENTS ' r 

0//»2X? 'XDIR' f6XfF10»2f 15Xf F7*3f // f 2Xf ' YDIR' f 6Xf F10*2f 15Xf F7*3f 
8//f2X»'ZDIR'>6X*F10.2*15X*F7.3) 

RETURN 

END 

SUBROUTINE BURN 

SUBROUTINE BURN 

CONTROLS INPUT OF BURNOUT DATA AND 
OUTPUT OF ORBITAL ELEMENTS 


REAL I 

COMMON/DATA/X »YrZr XDOT f YDOT » ZDOT » A r ECCEN f THET A f OMEGA 1 1 1 U 
WRITE<6f 200) _ 

; >0 FORMAT</f 5Xf ' IF THE COMPONENTS OF THE BURNOUT DATA ARE' 

8 'NOT KNOWN TYPE If ELSE TYPE 2*') 

READ<5f*)J 

IF<J.EQ.1)THEN 

PRINT *f 'INPUT THE ALTITUDE' 

READ ( 5 f HO R 

( ; PRINT *f' INPUT THE VELOCITY MAGNITUDE' 

E READ<5f*)V 

PRINT *f' INPUT THE ANGLE PSI' 
n- READ<5f*)PSI 

k X- < 3440 ♦ +R ) tCOS < PS I ) 

Y» ( 3440 ♦ +R ) *S I N < PS I ) *COS ( - . 4974 ) 
Z-<3440f+R)*SIN<PSI)#SIN<-.4974> 

F XDOT— l.*V*SIN< PSI) 

1 YDOT«V#COS < PSI )*COS<-. 4974) 

ZDOT— 1 . *V*COS < PSI ) *SIN < - . 4974 ) 

y> : ELSE 

F WRITE<6f 400) _ ,, 

400 F0RMAT(/f5Xf 'INPUT OF THE BURNOUT DATA IN COMPONENT FORM'»//» 

? 05X f ' X- ' ) 

F READ<5f*)X 

WRITE<4f 405) 

405 FORMAT < / f 5X f ' Y- ' ) 

| READ(5f*)Y 

URITE(6f 410) 

410 FORMAT (/fSXf'Z-') 
f ; READ<5f*)Z 

1 WRITE< 6»415> 

415 FORMAT < / f 5X f ' XDOT * ' > 
r READ < 5 r * > XDOT 

WRITE<6f 420) 

420 FORMAT < / 1 5X r ' YDOT- ' ) 

READ < 5ft) YDOT 
j;* WRITE<6f 425) 

425 FORMAT < / f 5X f ' ZDOT- ' ) 

READ (5ft) ZDOT 
r END IF 

( CALL ORB-ELEM 

WRITE < 6 f 440 ) A f ECCEN f THET A f OMEGA f I f U 
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0 FORMAT <//f5Xr'A*'fF12*3r//r5Xf 'E-' rF6.4r//r5Xt 'THETA *' t 
8FA*4»//»5Xf 'OMEGA*' fF6«4f//f5Xf'I*'»F£»4f//f5Xf'U*'fF7*4> 
><?9 RETURN 
END 

SUBROUTINE ORB-ELEM 

SUBROUTINE ORB-ELEM 

CALCULATES THE ORBITAL ELEMENTS USING 
THE BURNOUT DATA 


REAL I» KM 

COMMON/DATA/X t Y f Z t XDOT t YDOT r ZDOT r A r ECCEN t THET A r OMEGA 1 1 r W 
KM-62744.8 

V-SQRT < XD0T#*2+YD0T**2+ZD0T#*2 ) 

R-SQRT < X##2+Y**2+Z#*2 > 

C0NST-R#V**2/KM 
A»R/( 2-CONST) 

H-SQRT ( < Y*ZDOT-YDOT*Z > #*2+ < XDOT*Z-X*ZDOT > ##2+ ( X*YDOT-XDOT*Y ) *#2 ) 
ALPHA*ASIN< <X*XDOT+Y*YDOT+Z#ZDOT)/<R#V) ) 

ECCEN*SQRT < ( CONST-1 ♦ ) ##2*C0S < ALPHA > ##2+SIN < ALPHA > *#2 > 
boO THETA*ACOS< <H**2/KM/R-1 ♦ ) /ECCEN) 

OMEGA-AT AN < < Y*ZDOT-YDOT*Z ) / < X*ZDOT-XDOT#Z ) > 

' I*ACOS<<X#YDOT-Y#XDQT)/H) 

L. U»AT AN < < -1 . #X*SIN < OMEGA ) #COS < I ) + Y*COS < OMEGA ) *COS <I)+Z*SIN<I)>/ 

8 < X#COS < OMEGA ) +Y*S I N < OMEGA ) ) > -THET A 
RETURN 
[ v END 

SUBROUTINE BURNOUT 
f2 SUBROUTINE BURNOUT 

H CALCULATES THE BURNOUT DATA FROM THE 

ORBITAL ELEMENTS 

i 

REAL I r N»L1 r Ml *N1 » L2»M2 »N2» KM 

I . COMMON/DATA/X f Y r Z r XDOT f YDOT t ZDOT t A f ECCEN t THET A t OMEGA t AI r AU 

KM-62746.8 
EL*A#<1.-ECCEN**2> 

R-EL/ < 1 . +ECCEN*COS < THETA ) ) 

B* ASSORT ( 1 . -ECCEN##2 ) 

N*SQRT<KM/A##3) 

V-SQRT < KM# ( 2 . /R-l ♦ / A ) ) 

I E* ACOS < 1 . /ECCEN-R/ < A#ECCEN ) ) 

L 1-COS < OMEGA ) *COS < AU ) -SIN < OMEGA ) #SIN < AW ) «COS ( AI ) 

MI-SIN < OMEGA ) *COS < AU ) +COS < OMEGA ) «SIN ( AU ) «COS ( AI ) 

N1*SIN< AW) #SIN( AI ) 

L2— 1 . #COS ( OMEGA ) #SIN-< AU ) -SIN < OMEGA ) #COS < AU > #COS < AI > 

1 M2— 1 . #SIN ( OMEGA )#SIN< AU) +COS < OMEGA ) #COS ( AU > #COS ( AI ) 

N2-C0S ( AU ) #SIN ( AI ) 

X*A#L1 *COS < E ) +B#L2#SIN ( E ) -A#ECCEN*L1 
Y*A#Ml#COS ( E ) +B*M2#SIN < E ) -A#ECCEN#M1 
Z-A#Nl#COS < E > +B*N2#SIN < E ) -A*ECCEN#N1 
t CONST»N*A/R 

p XDOT -CONST# ( B#L2#C0S (E)-A#L1#SIN(E) ) 

YDOT-CONST# < B*M2*C0S < E ) -A#M1#SIN < E ) ) 

ZDOT -CONST# ( B#N2#C0S ( E ) -A*N1 #S I N ( E ) ) 
li RETURN 71 

^ END 


FLOW CHART FOR ASTRO. FOR 
MAIN PROGRAM 









ORIGINAL PAGE IS 
OF POOR QUALTCB 


FLOW CHART FOR ASTRO, FOR 
SUBROUTINE BURN 



FLOW CHART FOR ASTRO, FOR 
SUBROUTINE ORB-ELEM 
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FLOW CHART FOR ASTRO* FOR 
SUBROUTINE ORBIT 



(return ) 


FLOW CHART FOR ASTRO. FOR 
SUBROUTINE BURNOUT 


(begin j 


CALCULATE 

XrYf ZrXBOT »YBOT rZDOT 


(return]) 
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APPENDIX G - REFERENCE MATERIAL 
RECEIVED FROM MARSHALL SPACE FLIGHT CENTER 

1 . MASS REPORT 
2* DESIGN DATA 
3. GROUP WEIGHT STATEMENT 


75 




ORIGINAL’ PAGES 13 
OF POOR QUA LTPH 


>******************* 
•HASS REPORT* 
« •* ••• • •••«•>• •*>••• • * *•• 


t.v 


I? 


AEROSHEU. CRAPHITE-POLYIHIDE SSTO SOT LOX-RP/TOTAL 
INITIAL GUESS VERSION S. 1 


. 1.0 VINO GROUP 
. 2.0 TAIL GROUP 
0.0 BOOT GROUP 

BASIC STRUCTURE 
THRUST STRUCTURE 
RP-1 TAMC 
LOX TANK 
LH2 TANK 
BOOT FLAP 

4.0 IWUCED ENVIR0T#€NT 
S.B L AWING GEAR 

6.8 PROPULStCN 

7.0 PROPULSION. RCS 
B.O PROPULSION. OHS 
0.0 PRU£ POWER 

o.i fuel cells 
0.2 reactant oevars 
0.3 REACTANTS 
0.4 BATTERIES 
IB.B ELEC CQNV AW OISTR 
It.B CONTROLS 

11. 1 HYOAAULICS 

11.2 SURFACE CONTROLS 

13.0 AVIONICS 

14.0 EMVIRO ACNTAL CONTROL 

15.0 PERSOMCL PROVISIONS 

10.0 NARGIN 
C RT WEIG HT 

17.0 PERStMCL 

10.0 resiooal fluios 

LANDED WEIGHT V/0 CARGO 

20.0 CARGO (RETUmCDI 
LAWED WEIGHT 
ENTRY WEIGHT 

23.8 ACPS PROPELLANT 

RCS 

24.8 CARGO DELIVERED 

26.8 ASCENT RESERVES 

20.8 INFLIGHT LOSSES 

27.8 ASCENT PROPELLANT 

RP-I 1S02U3. kg 

LIFT OFF WEIGHT 


IS566. Vg 
4340. Wg 
1806. kg 
8S00. kg 
7680. kg 
414. kg 


ITS’ Kg 
22*. K« 
207. kg 
306. kg 

B. kg 
248S. kg 


1326. 

I28C8. 


kg 


4334. 

1234. 

38405. 


!3712. 

3377. 

18200. 

1265. 

637. 

054. 


456B. 

2405. 


2240. 

2000. 


75B7. 

m44S - .Ha. 

6007. 

I006SS. k a 
15420. 
125075. kg 
I 250 75. kg 
14134. 


LH2 68074. kg 


kg 


v® 

s 

S 

kg 


kg 

kg 

Jf 

§ 

kg 

J* 

kg 


kg 


127§. kg 
1031. kg 
1062542. kg 
LOX 8442US. kg 

1204056. kg 


<0.806 ( f 


<0.000 1) 


<8.009 

<0.000 


t ) 
1 ) 


<0.080 ) ) 


LV/S 

V-..V 


1 /*.; 
r :-i. 
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ORIGINAL PAGE IS 
OR POOR QUALITY 


tsioN e. t 



I cum 
I CUM 
CUM 

CUM 


I N 
i N 

i 


NODE 2 


76401 ss.ee N 

oeer678.ee n 
e.eeiS27ae 


cun 

kg/cu* 

cun 

kg/cuM 

cun 

kg/cuM 


ORIGINAL PAGE IS 
OF POOR QUALITY 




z w 


o 

ac 

a 


v (ft o 

UXOte 
WOW . 
mm (ft 
(ft z 


«• (ft 


u •- 
v* 



N o 


z 

(• 


xxs«ooac« 
ac 9 mm u<<3C 

***»<QO*® 


UZ222«X»a>KZ233 “ 

§4444 i2203**»ee 2 

^ 8 »>»>•■ ft Z 22 Z ft J 


I 0 W 8 >* 

•tv 


*- 2 IKX 


U (3 


^ y ^ ^ _ 

SjjjJ 0 4*1444 5 *“ 

9 4 4 4 ft U • •• 

c ^ £ a ft k & 
ftftft^u J 



O XX 
ftUKW 
O C ft 

IlftU 


o 

y 


fluom 4 AOft ®u^y 


yOiM < y ft uOm» 


SDT2-SH-0120-33 
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Mi|«HT X«C6 V*CC I»C Q 

FUNCTION FUNCTION OfSC'l'TlON (|U) (IN) (IN) (IN) 


ORIGINAL PAGE IS 
OE POOR QUALITY 


) 





TM UMHRftUl) BODY 041$, 0 . 004, • 0,0 HO, 4 

FORWARO (HOOVI 1700*0 174,0 0,0 JJ0,| 

NIOOOOV 1771,0 011,0 0,0 200,1 

HIDBOOV-PaVIOAO BOV DOORS *70,0 007,1 0,0 O4S,0 
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ABSTRACT 


This report describes progress made to Improve Initial design 
considerations for the Intermediate-Orbit Cargo Vehicle. I.t provides 
a comprehensive overview of mission profile, structural design and 
cost analyses, and relates these to the overall feasibility and 
usefulness of the proposed system. 
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INTRODUCTION 


The decision to Initiate design of a new cargo vehicle was based on 
predictions of Imminent near-space equipment needs to support the space 
station. As the US presses on towards the 21st century, many new aerospace 

projects will require payload capabilities beyond that of present launch 

\ 

vehicles. By devising a vehicle capable of transporting cargo to an 
Intermediate orbit and In essence "parking* It In that orbit, the US will 
achieve the capability of having "orbital warehouses" where material and 
supplies for orbital operations can be stored. 

The mission requirements established for this vehicle are shown In Figure 
1. They call for an unmanned vehicle that will place a payload of 150,000 lbs 
In a 150 nautical mile orbit. The payload bay size must be 25 ft by 90 ft or 
33 ft by 100 ft. The staging velocity must be a maximum of 7000 ft/s. The 
fly-back-booster which will launch the cargo vehicle Into orbit must see a 
maximum of 3 g's. The proposed configuration Is shown In Figure 2. 

After separation from the fly-back-booster, the restriction that 
acceleration be held to within 3 g’s for human tolerance Is removed since the 
cargo vehicle Is unmanned. The vehicle will sustain a maximum of 5 g's. All 
preliminary specifications were based on data obtained from the Saturn and 
Apollo projects. Subsequently, as a result of continuing research and follow- 
on design considerations, the accuracy of design specifications has Increased 
drastically. Research on the cargo vehicle Includes precise trajectory 
analysis, estimates of structural sizes and weights, fuel tank sizing, and 
estimated costs of production. 

As stated earlier, the cargo vcehlcle's mission requires that It be 
placed In orbit at an altitude of 150 nautical miles. After It has achieved 
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the required orbit, the nose cone will detach from the cargo section and will 
burn up as It re-enters the atmosphere. The cargo will then be removed 
through the open end of the payload bay. An orbital tranfer vehicle (OTV) 
will attach to a docking mechanism (Figure 3) similar to the one used on the 
Saturn Command Module. Docking will be achieved by maneuvering the OTV close 
enough to the multiple docking adaptor (MDA) 'located on the first pqyload 
frame's forward spider beam so that the probe on the OTV engages with a drogue 
on the multiple docking adaptor. When the probe comes into contact with the 
drogue It will be guided Into a socket at the bottom of the drogue. Three 
capture latches In the probe head will secure the OTV to the first payload 
frame. The OTV will withdraw the first payload frame from the payload bay 
(Figure 4), towing the remaining frames forward until the second frame Is In 
place In the space or Igunally occupied by the first frame. The OTV will then 
separate the first payload frame from the second payload frame and proceed to 
the space station with the first payload frame in tow. After transporting the 
first payload frame to the space station, the OTV will return to the cargo 
vehicle and engage the MOA on the second patload frame's forward spider beam. 
The OTV will remove the second payload frame, tow it to the space station, and 
then return to the cargo vehicle. It will repeat this off-loading procedure 
until all the cargo Is secured atthe space station. After cargo off-loading 
Is completed, explosive bolts will separate the upper and lower sections, each 
of which will be equlped with an MOA, for transportation to the space station. 
This will save the fuel tanks, flight Instruments, engines, and empty payload 
bay for future use. 

Although It would be more economical to tow the entire cargo vehicle, 
cargo Included, to the space station in one trip, considerations for the capa- 
bilities of the OTV dictate that the components be transferred individually. 
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TRAJECTORY 


A direct ascent trajectory will be used to place the vehicle In the 
parking orbit. In this approach, the engines burn continuously fro* lift-off 
to orbital Insertion. When burnout occurs, the vehicle Is at circular 

satellite speed and the angle between the orbital direction and the velocity 

\ 

vector of the vehicle Is zero. When using this approach, angle and velocity 
are critical. If the angle or velocity Is incorrect, the vehicle will not 
attain the proper orbit. Therefore, extra fuel space has been made available 
so that the engines can be used for orbital corrections. Reaction controls 
located In the nose cone will be used to pltsh the vehicle so that errors In 
the ascent trajectory can be corrected by re-firing the engines for short 
periods of time. 

Using a TRAJECTORY program, the trajectory data were calculated. For 
brevity, a summary Is given here and Important events are Illustrated in 
Figure 5. The program and the flow chart are shown In Appendix A. 

Vehicle Mass - 383,247.3 
Payload Mass - 150,000 lbs 
Propellant Mass - 1,108,133 lbs 
Structural Mass • 233,247.3 lbs 
Initial Mass 

at Separation - 1,491,380.3 lbs 
6-Turn Altitude - 1800 ft 
Thrust - Cargo - 1,531,100 lbs 
• Total at 

Lift-off - 10,221,000 lbs 
Separation Velocity • 7080 ft/s 
Separation Altitude - 46.92 nml 


- 3 - 



Separation Range - 38.87 nml downrange 
Burnout Velocity - 24,352 ft/s 
Burnout Altitude • 159.4 rani 
Burnout Range - 720 nml downrange 

Referlng to Figure 5, point A denotes the commencement of the gravity 

\ 

assisted turn. Point B denotes the point of booster separation and point C 
denotes the point of orbital Insertion. 

The next step was to Investigate the advantages of the conventional bum- 
coast-burn trajectory. In this approach, the vehicle's engines bum at a very 
high rate up to a certain point. This maneuver places the vehicle Into an 
ascent ellipse with apogee at the desired orbital altitude. When the vehicle 

reaches the apogee of the ascent ellipse, the engines are again fired to 

» 

accelerate the vehicle to circular orbital speed and to align the vehicle with 
the direction of the orbit. It was Initially presumed that lift-off mass 
might be reduced significantly through use of this approach. However, 
subsequent analysis Indicated .that use of the burn-coast-burn trajectory would 
not result In a significant reduction of lift-off mass. 



STRUCTURES 


Continuing development of the cargo vehicle led to extensive research In 
the area of vehicle structures. Originally, structural weights were assumed 
and structural materials were not even considered. In addition, two cargo 
bays were Initially Included in the design. However, because of design re- 
specification, only one cargo bay Is used In the final design. Specific 
calculated structural sizes and weights of the lower and upper sections are 
listed below and shown In Figure 6. 

LOWER SECTION 
Aluminum 2014 Tb Alloy 

Overall Size - 89.0 ft (height) 

- 34.0 ft (diameter) 

Fuel Tanks - 67.0 ft (height) 

- 33.0 ft (diameter) 

Engine Compartment 

and Engine Cones - 19.0 ft (height) 

- 34.0 ft (diameter) 

Instrument Unit and 

Separation Wall - 3.0 ft (height) 

- 34.0 ft (diameter) 

Overall Weight - 109,913 lbs (Est. 113,368.78 lbs) 

Fuel Tanks (Dry) - 36,500 lbs 

Skin - 29,845 lbs at .25 In thickness) 

Instrument Unit - 6,000 lbs 
Thrust Structure - 7,000 lbs 
Engines (4) - 28,568 lbs 

Other - 2,000 lbs (Pipes, Insulation, Wiring, Bulkheads) 
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UPPER SECTION 


Aluminum 2014 Tb Alloy 

Overall Size - 101.0 ft (height) 

- 34.0 ft (diameter) 

Skin - 101.0 ft (height) 

• 34.0 ft (diameter).^ 

Payload Structure • 100.0 ft (height) 

- 33.0 ft (diameter) 

Spider Beams (2) - 1.0 ft (height) 

- 33.0 ft (diameter) * 

Overall Weight - 123,334.0 lbs 

Skin - 47,065.0 lbs at .375 In (thickness) 

Payload Structure > 37,715.0 lbs 

Spider Beams (2) - 4,000.0 lbs 

Payload Containers - 9,127.0 lbs 

Cone - 12,000.0 lbs 

Interstage - 6,000.0 lbs 

Restraints - 7,424.0 lbs 

TOTAL STRUCTURAL SIZE AND WEIGHT 
Total Structural Mass - 233,247.0 lbs 
Total Structural Height - 223.0 ft 
Lower Stage - 82.0 ft 
Payload Stage - 101.0 ft 
Interstage - 3.0 ft 
Cone - 30.0 ft 
Engine Cones - 7.0 ft 
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The structure of the cargo vehicle can be divided Into three main 
categories: the lower section, the Interstage, and the upper or payload 
section. 

LOWER SECTION 

The lower section ,as shown In Figure 7, Is composed of an Instrument 
unit, a tall unit, a propellant tank unit, and an engine unit. 

The Instrument unit (Figure 8) Is an unpressurized, cylindrical, load- 
supporting structure. It will be 34.0 ft In diameter and 3.0 ft high. It 
will weigh approximately 6,000 lbs. 

The Instrument unit will be constructed In three segments to facilitate 
shipping and packaging. The structure consists of Inner and outer stage 
skins, a honeycomb core, and forward and aft Innerstage rings. Splice plates 
connect the three segments. After assembly of the Instrument unit, a door 
designed to act as a load-carrying member during flight provides access to the 
electronic equipment Inside the structure. 

The Instrument mounting panels, which are used for mounting various 
components of systems housed In the Instrument unit, attach to welded brackets 
around the Inner periphery of the Instrument unit. 

The Instrument unit houses electrical and mechanical equipment. This 
equipment Is used for guidance and control and monitors the vehicle's flight 
performance from lift-off to orbit. The Instrument unit houses several 
systems such as guidance and control, measuring and telemetry. Instrument uni* 
command, tracking, and electrical systems. 

The tall unit (Figure 9) Is composed mainly of the thrust structure. 

This structure provides a mounting point for the four engines and the 
Instrument unit. Engine thrust Is transmitted through the thrust structure to 
the Instrument unit and then up through the propellant tank system. The 
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forward end of the tall unit Is enclosed by fire walls and the aft end Is 
enclosed by heat shield panels. 

The thrust structure consists of the barrel assembly and thrust supports. 

The instrument unit fire wall panels are also housed In the thrust structure. 

The bottom of the barrel assembly will contain the heat shields for the four 

\ 

engines. 

The heat shield panels seen In Figure 9 enclose the aft end of the tall 

i 

unit assembly to form an engine compartment that protects the engines from the 
heat of engine exhaust. Each panel Is a layer of corrosion-resistant steel 
honeycomb sandwiched between two layers of corrosion-resistant steel sheet. A 
plaster-type thermal Insulation Is trowelled Into the cellular openings of the 
honeycomb. 

The flame shield closes the area between the engines. The flame shield 
Is made of corrosion-resistant steel. All parts of the flame shield that come 
In contact with the engines are Insulated. 

Flame curtains will be Installed between the heat shield panels and the 

# 

four engines' thrust chambers. This will prevent engine exhaust flame from 
entering the engine compartment. The curtain will be made of silicon rubber 
Insul at Ing jnaterl al sandwiched between a double thickness of fiberglass cloth 
tape on the side closest to the engines. 

The propellant tank unit (Figure 7) will house the propellant tanks. 

These propellant tanks will Incorporate the principle of the common bulkhead, 
which comprises the top half of the liquid oxygen tank and the bottom of the 
liquid hydrogen tank. 

The common bulkhead causes problems because the H* is so much colder than 
the 0 t ; therefore, the bulkhead will be well Insulated between the two 
ellipsoid sphere halves forming It. The Insulation will be a honeycombed 

i 

I 
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phenolic material that Mill be fitted between the base of the IH« tank and the 
top of the LOX tank. It Is estimated that the common bulkhead will save about 
12.0 ft In length and 12,250.0 lbs In weight. For Information on tank sizing 
see Appendix B. 

The propellant pipes for the H 2 will run -through the center of the 0 X 
tank and Instrument unit to the engines. They will be insulated. The pipes 
for the 0 a tank will also run through the center of the Instrument unit. 

Fins will not be used on the cargo vehicle. It was determined that 
during the Initial boost phase the combined booster and cargo vehicle could be 
stabilized by the use of the combined thrust vectoring capabilities and the 
lifting surfaces of the booster. With 10.5 degrees total pitch and 8.5 
degrees horizontal yaw, the resultant force can be made to counteract the 
moment produced by the adverse difference in the center of pressure and the 
center of gravity. 

After separation from the booster, the atmospheric conditions at 47 
nautical miles altitude are such that the aerodynamic effects are very minute. 
Therefore, stability must be maintained using non-aerodynamlc effects (l.e. 
thrust vectoring). 

The engine unit (Figure 7) consists of the propulsion system. The lower 
stage will be powered by a cluster of four engines. Each engine develops 
395,000 lbs thrust for a total thrust of 1,580,000 lbs for the whole rocket 
(see Figure 10). The four engines will be glmballed for thrust vectoring. In 
addition, each engine will have an extendable nozzle for added thrust after 
separation. 
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INTERSTAGE 


The Interstage forms a structural Interface between the lower section and 
the upper or payload section. It Is a cylindrical, skin-stringer-type 
aluminum structure. It Is 3.0 ft high and 34.0 ft In diameter. It will weigh 
approximately 6,000.0 lbs. 


UPPER OR PAYLOAD SECTION 

. The payload section (Figure 11) will be of aluminum skin-stringer-type 
construction. It wlll-be 101.0 ft tall and have an outside diameter of 34.0 
ft. It will be constructed In two parts, an outside shell and an Internal 
frame which will hold the payload. 

The outside shell will have an aluminum skin, 8 horseshoe beam longitudi- 
nal members, and nine rings. The 8 horseshoe beam members and skin will 
provide most of the structural Integrity of the system. They will also 
provide support for the payload frame tracks. The nine rings will provide 
support for the system. The 8 horseshoe beams will be supported by aluminum 
rings on both ends. The aft ring will be attached to the Interstage. The 
outside shell will be 101.0 ft long and have an Inside diameter of 33.5 ft 
(see Figure 10). 

The payload frames (Figure 12) will each have 8 I-beam type longitudinal 
members and will feature circumferential aluminum bands placed at regular 
Intervals to provide additional structural Integrity. The. fore and aft ends 
of each payloid frame will be capped with spider beam supports. The aft spider 
beam of the rearmost payload frame will help transfer weight and thrust 
between the Interstage and the payload bay. Each payload frame will have an 
Inside diameter of 33.0 ft and an outside diameter of 33.25 ft. 

The 8 horseshoe beams of the outside shell will act as tracks for 
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movement of the payload frame I-beam members during cargo off-loading. The 
tracks and I-beams will be coated with Teflon and sprayed with silicon to 
provide smooth, virtually frlctlonless movement. 

The nose-cone will be made of aluminum skin-stringer-type construction. 
It will serve as an aerodynamic fairing and will also house the vehicle's 
reaction controls. The nose-cone will be 34.0 ft In diameter at the base and 
30.0 ft high. It will have the shape of a tangential ogive. 
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FUEL SECTION 


Knowing that the total mass of the propellant Is 1,108,133 lbs and that 
the mixture ratio Is 1 part oxygen to 6 parts hydrogen, the masses of the fuel 
and oxidizer were found. The mass of the oxygen was determined to be 
approximately 949,828 lbs and the mass of thd hydrogen was calculated to be 
approximately 158,304 lbs. The density of the oxygen will be 71.2 lbm/ft 5 and 
the density of the hydrogen will be 4.43 lbm/ft*. By dividing the mass of 
each propellant component by Its particular density, the volume necessary to 

r 

accommodate the required amount of propellant can be determined. 


949,828 lbs 

VEo, ■ -- - 13,340.28 ft 

71.2 lbm/ft* 


V 


158,304 lbs , 

» 35,734.698 ftr 

4.43 lbm/ft* 


The oxygen will be held In an ellipsoid shaped tank with radius equal to 16.5 
ft and height equal to 23.4 ft. For details on the tank sizing see Appendix B. 

The hydrogen will be held In a cylindrical tank joined In tandem with the 
oxygen tank (Figure 7). The upper end of the cylinder will be closed with a 
curved cap equal In volume to half that of the oxygen tank. The bottom end of 
the cylinder will be joined to the oxygen tank, which will Intrude half Its 
volume upon the cylinder, reducing the cyll idler's volume by an amount equal to 
half the volume of the oxygen tank. The ne effect of the upper cylinder cap 
In tandem connection with the oxygen tank will be to provide space for a 
volume of hydrogen equal to that available In the standard, unenhanced 
cylinder. To determine the total weight of the fuel tanks, the heights of the 
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ellipsoid and the enhanced cylinder were combined. As stated previously, the 
ellipsoid height was determined to be 23.4 ft. The height of the cylinder was 
calculated from the cylindrical volume equation (Appendix B) and was found to 
be 41.78 ft. Thus the total fuel section height will be -65.18 ft, rounded up 
to 67.0 ft to accommodate orbital maneuvering purposes. 
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COST ANALYSIS 


The costs for the vehicle were determined using the Launch Vehicle Cost 
Model program at the Marshall Space Flight Center. A compilation of all data 
used In the program and the results are presented In Figure 13. 

The structural weight used In the program was 187,288 lbs and Includes 
the vehicle skin, fuel tanks, Internal supports, pipes, Insulation, bulkhead, 
wiring, frame, forward and aft skirts, tunnels, Interstage, and nose fairing. 
The base heat shield protects the thrust structure, propellant tanks, and 
other stage elements from excessive heat and gases given off by the engines 
during their burn. A weight of 1,000 lbs was used for this component In the 
program. The avionics, whose weight was estimated to be 250.0 lbs, will be 
composed of the following components: (1) Guidance, Navigation , and Control, 
and (2) Communications and Data Handling. The electrical power equipment 
weight used In the program was 500.0 lbs. This weight will be comprised of 
batteries of fuel cells, a power distribution system, relays, Inverters, power 
transformers, cables, and wiring. The propulsion components other than the 
engines were estimated to weigh a total of 1500.0 lbs and will Include a 
propellant feed system, recirculation system, propellant management system, 
pressurization system, plumbing, valves, lines, and glmballlng bearings. The 
total cost for the vehicle elements other than engines Is $1.12473 billion 
during the design, development, test, and evaluation (DDT&E) stage and $41/. 21 
million for the total first unit (TFU). 

The structural tooling, systems test hardware and assembly, and systems 
test operations will cost an estimated $2.60864 billion for DDT&E. Ground 
support equipment (GSE) during DDT&E will have an estimated cost of $453.44 
million, while systems engineering and Integration Is estimated at $322.55 
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million during DOT&E and $34.12 million for the TFU. The estimated cost for 
the program management Is $129.65 million during DDT&E and $16.47 million for 
the TFU. 

The engines and Incidental costs are the final considerations. The 
engine parameter used In the program was calculated by dividing the product of 
weight, vacuum thrust, specific Impulse, andx chamber pressure by 10,000,000 
(this parameter Is known as the "composite variable"). Use of this parameter 
Jn the program yielded estimated costs of $44.1 million during DDT&E and 
$40,000 for the TFU. The Incidental costs Include a contractor's fee of 
$498.17 million during DDT&E and $65.5 million for the TFUf program support 
costs of $121.7 million during DDT&E and $10.67 million for the TFU; and cost 
contingencies of $626.74 million during DDT&E and $81.6 million for the TFU. 

A simplicity factor of 1 was used In the cost analysis program since our 
vehicle was designed based on previous Saturn project data. Also, surface 
flight controls and hydraulics were neglected due to the elimination of 
aerodynamic fins from the final design. 
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CONCLUSIONS 


The goal of the ongoing research described In this report Is to perfect 
the design and production of an unmanned, Intermediate-orbit cargo vehicle for 
use with fly back booster. Recent work has led to significant enhancements In 
the projected capabilities of the vehicle, Including an Improved mission 
profile, a cost-efficient development and production plan, and an Innovative 
approach to cargo storage and transferral. 

The overall design of this cargo vehicle meets and exceeds all required 
standards set by NASA. It can be seen that aerospace projects of this nature 
will lead to the advancement of the US space program In the near future. 
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GROUND RULES FOR CARSO VEHICLE 

•UNMANNED 

-PLACE 150,000 1b PAYLOAO IN 150 NAUTICAL MILE ORBIT 
-PAYLOAO BAY: 25*0 X 90'L OR 33'D X 100'L 
-PROPELLANT: LIQUID OXYGEN/LIQUID HYDROGEN 
-FIRST FLIGHT IN 1998} DESIGN FREEZE IN 1990 


Figure 1. Mission Requirements 
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FIGURE 5. Graphical Representation of Trajectory 
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FIGURE g INSTRUMENT UNIT 
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Figure 10. Space Transportation System Main Enqlne 



Fig. II. Cutaway of: Payload Section 
With Fbyldad Frame 
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Figure 12. TV© payload frames in docked configuration. 
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This is a general purpose trajectory analysis program utilizing 
the gravity of the Earth to achieve the\ proper attitude for insertion 
into a low Earth orbit* This maneuver is called a gravity assisted 
turn maneuver* Since this program was used to model the trajectory 
for an unmanned cargo vehicle; the forces due to high acceleration 
were of secondary importance; and not limited* The cargo vehicle 
was used in conjunction with a reuseable flyback booster which 
returns to a landing strip after staging from the cargo vehicle* 

The program is set up so that staging does not occur until the 
total configuration is accelerated to 7000 feet per second* As 
another limit the program ceases execution if all of the mass is 
burned as fuel for obvious reasons* 


ALT =0 
DT=0 ♦ 1 
RANG=0 
ISP1=320 
ISP2=380 

INPUT-ENTER INITIAL MASS' ?M1 
M=M1 

INPUT -ENTER ALTITUDE FOR G-TURN' ; AGT1 

AGT=AGT1 

TH=1.35*M 

INPUT-ENTER LIFTOFF THRUST (CARGO ENGINES )•; TH21 

TH2=TH21 

TH1=TH-TH2 

MD0T1=TH1/ISP1 

MD0T2=TH2/ISP2 

ISP=TH/(MD0T1+MD0T2) 

MD0T=TH/ISP 

MC1=0 

MC2=0 

Y =0 

X=0 

60=32.2 
R0=2 ♦ 092E+07 
00=0 
0=00 
T=0 
PRINT 

PRINT-TIME 0 RANG ALT G'S MASS THRUST 

PRINT 

F$=-**#.* ******* ****** ****** *♦** ******* ****** 

M2=M-MD0T1*DT-MD0T2#DT 
G=G0# ( R0/ ( R0+ALT ) >'*2 
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300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 
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490 

500 
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520 
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550 

560 

570 

580 

590 
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620 

630 

640 

650 

660 

/ — * -\ 

<Li / v/ 

68C 

700 

710 

720 

730 

740 

750 

760 

770 

730 

735 

790 

300 

810 

820 

S30 


DV=ISP#G0*L0G < M/M2 > -G*DT 

AG1=DV/G0/DT 

REM IF AG1>3 THEN 400 

CHY=0 ♦ 5* < V0+V0+DV ) *DT 

IF Y+CHY>AGT THEN 580 

Y=Y+CHY 

ALT =Y 

V=YO+DV 

AG=AG1 

GOTO 460 \ 

DV=3*G0 

M2=M/EXP< <DV+G*DT>/ISP/GO> 

A=DV/GO 

V=VO+DV 

Y=Y+0.5*CV0+V>*DT 

ALT=Y 

T=T+DT 

DM=M-M2 

MDT=DM/DT 

TH=MDT*ISP 

MC1=MC1+TH1/ISP*DT 

MC2=MC2+(TH-THl)/ISP2*nT 

MC=MC1+MC2 

PRINT USING F$?T » V»RANG» ALT » AG»M»TH 

REM 

M=M2 

VO=V 

GOTO 280 
PRINT 

PRINT 'TIME V RANG ALT G'S MASS THRUST* 

PRINT*—. * 

M=MfMBQTl#BT +MDQT2#DT 

BPSI-< PI/180 )*0» 1 

PSI0-PI/120 

PGI=PSIO 

GAM=1 

2as=siN(Psio)/(i+cos(Psio) ) 

PRINT'TIME V RANG ALT PSI* GAM THRUST MASS G'S 

F$=**#« ***** ******* ******* **»* **♦* ******* ******* #.4 

PRINT 

N==TH/M 

PRINT USING F$ ? T * V » RANG 7 ALT - P7I*180/PI r GAM r TH ? M 7 AG 
C=VC/< Z0"<N-1) >/<l+Z0~2; 

?s:=psi+dpsi 

Z=SIN(PSI)/( 1+C0S(PSI) ) 

, n=C*(Z''<N-l) )*<1+Z~2) 

IF Vi:: 7000 THEN 1390 
y=vi 

BT= C/G*Z~ < N-l > * < 1/ < N-l > +Z~2/ ( N+ 1 ) ) -C/G*ZO~ < N-l ) * < 1 / ( N-l ) +Z0~2/ < N+ 1 > > 
BX=0 ♦ 5* < VOfcSIN < PSIO ) +V*SIN < PSI ) ) *DT 
DY=0 . 5* < VOfcCQS < PSIO ) +V*COS < PSI > ) * DT 
AG = ( V-VO ) /DT /GO 
IF AG>3 THEN 1060 
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840 X*X+DX 

84S Y=Y+DY 

850 THETA*ATN(X/<RO+Y) ) 

860 GAMMA*PSI-THETA 

870 GAM=GAMMA*180/PI 

880 ALT-<Y+R0) /COS < THETA ) -R0 

900 RANG=RO*THETA 

930 T*T+DT 

940 PSIO=PSI 

950 ZO=Z 

960 VO=V 

970 IF ALT<100000 THEN 1020 

980 ISP1=340 

990 ISP2=462 

1000 MD0T1=TH1/ISP1 

1010 MD0T2=TH2/ISP2 

1020 M=M-MDOTl*DT-MDOT2*nT 

1030 IF M<0 THEN 1390 

1040 G=GO#<RO/<RO+ALT) >~2 

1050 GOTO 720 

1060 N=3*G0/G+< <1-Z0~2)/<1+Z0~2> ) 

1070 C=V0/Z0~<N-1>/<1+Z0~2> 

1080 U=C*Z''(N-1>*< 1+Z~2) 

1090 D1=C/G*Z~ ( N-l ) * < 1/ < N-l ) +Z~2/ < N+l > > -C/G*Z0~ ( N-l > * < 1/ ( N-l ) +Z0~2/ < N+l ) > 

1100 DX=0.5*(V0*SIN(PSI0>+V*SIN<PSI> >*DT 

1110 DY=0.5#(V0*C0S<PSI0>+V*C0S<PSI> )*DT 

1120 X=X+DX 

1130 RANG==R0#7HETA 

1140 Y*Y f DY 

1150 ALT= < YtRO ) /COS < THETA ) -RO 

1160 T=T+DT 

1170 PSIO^PSI 

1180 ZO=Z 

1190 VO=V 

1200 TH--=N*M 

1210 TH1=TH-TH2 

1220 MDQT1 -TH1/ISP1 

1230 MD0T2=TH2/ISP2 

1240 M=M-MD0T1#DT-ND0T2$DT 

1250 MC1=MC1+MDQT1*DT 

1260 MC2=MC2+MDCT2*D7 

1270 MC~MC1+MC2 

1230 IF M<0 THEN 1390 

1290 G=G0*CR0/;R0+ALT) >"2 

1300 PSI=PSI+DPSI 

1310 THETA=ATN<X/(RO+Y) ) 

1320 GAMMA*PSI-THETA 

1330 GAM=GAMMA*180/PI 

1350 PRINT USING F$ J T * 'Jr RANG r ALT r PSI t GAN > TH ? M » AG 

1360 IF V>7000 THEM 1390 

1370 Z=SIN< PSI )/<l+COS(PSI) ) 

1380 GOTO 1060 

1390 PRINT 
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1391 

1392 
1400 
1 4*?0 

ORIGINAL PAGE IS 
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PSI^PSI-DPSI 
PRINT VI 

PRINT'TIME V RANG ALT PSI GAM THRUST 

PPTM7* - 

MASS G'S* 

T1IT I 

1 Ti. V 

1421 

PRINT 


1422 

PRINT* SEPARATION 

■ 

1423 

PRINT 


1 424 

PPTWT* 

-r- -T- * 

4. tin 

1430 

r r\ j. r« » 

PRINT’time v rsnS sit psi . Sam2 thrust 

mass S's* 

1431 

PRINT \ 


1432 

***** ******* ******* ***** ***** ******* 

******* ***** 

1 435 

M=M-1020000 


1440 

ISP=ISP2 


1450 

' TH=TH2 


1460 

MD0T= TH/ISP 


1470 

V0=V 


1480 

PSI0=PSI 


1490 

GAM2=PI/2 + THETA 


1300 

GAMM=GAM2*130./PI 


1510 

zo-sin ■: psi o > / < l+cos ( psi o > > 


1520 

N=TH/M 


1530 

PRINT USING F$ J T ?V» RANG* ALT PSI#180*/PI » GAMM* TH * M * AG 


1540 

C=V0/ (Z0'*(N-1) )/(l+Z0"2) 


1550 

PSI=PSI+DPSI 


1560 

THET A=ATN (X/ ( R0+Y ) ) 


1570 

GAM2=PI/2+T'rlETA 


1580 

GAMM=GAM2*180*/PI 


1590 

IF PSI1GAM2 THEN 2030 


1600 

Z=SIN ( PSI ) / < 1+CCS < PSI ) ) 


1610 

V=C3*c ( Z~ C N-l > ) * ( 1+Z''2 ) 


1620 

IF V>250C0 THEN 2080 


1 630 

DT=C/G*Z~ (N-l)*(l/< N- l ) ~!-Z"*2/ (N+l) > -C/G*Z0~ I K -1 ) * ( 1/ < f 

2 -1 ' fZ0~2/(Nrl > > 

1640 

DX~ . 5#CV0#SIN(PSI0) +VKSIN ( PSI ) ) *DT 


1650 

DY= ♦ 5* ( V0*CQS < PSIO ) +V#C0S ( PSI ) ) #DT 


1660 

AG= < V-VO ) /DT/GO 


1680 

X-X+DX 


1 

• / V \/ 

Y=Y+DY 


1701 

THETA=ATN ( X/ ( RO+Y ) ) 

• 

1702 

RfiNG=RO*THETA 


*? "J « 

.i. / j. .* 

ALT = ( Y+RO > /COS ( THETA > -RC 

1 

^ -«n M 
-*' Xm \J 

T-TfDT 


1730 

PSIO- PSI 


1740 

“T ... -» 

Xm W *■’ Lm 


1730 

V0---V 


1760 

M*M-MD0T*BT 


1770 

IF M <0- THEN 2080 


1780 

G=G0* ( RO/ ; R0+ ALT ; ) "2 


1790 

GOTO 1520 


2080 

REM 


2081 

PRINT 

1 

2090 

PRINT * INITIAL MASS J*)M1 


2100 

PRINT* ALT FOR G-TURN S'JAGTl 


2110 

PRINT* THRUST (CARGO VEHICLE) J*iTH21 


. • 
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r j r j r j r j r j r j m 


111 PRINT'FINAL VELOCITY t'JV 

112 PRINT'FINAL MASS t * 5M 

113 PRINT'FINAL PSI J ' *PSI*180/PI 

114 PRINT'FINAL ALTITUDE ♦ * j ALT/60G0 

115 PRINT'FINAL GAMM J'JGAMM 

116 INPUT'TRY AGAIN Y=1 N*2'»TRY 

117 IF TRY<2 THEN 15 
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FLOWCHART FOR TRAJECTORY ANALYSIS 
PROGRAM • CTRA J* 

r 
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B 


IS 

THE ALTITUDE 
'GREATER THAN 100>000* 
FEET ? 


DECREASE FUEL MASS FLOW TO 
MAINTAIN CONSTANT THRUST, 


IS THE MASS 
LESS THAN ZERO 
y 


IS THE 

VELOCITY GREATER 
THAN 7000 FEET/ 
SECOND ? 


IS THE 

MASS LESS THAN 
ZERO ? 
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PRINT FINAL TIME» VELOCITY t RANGE r ALTITUDE* 
FLIGHT PATH ANGLE* THRUST * MASS* AND 
G-FORCE 
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APPENDIX B 


FUEL TANK DATA CALCUIATICNS 


1 ,108,133 Ibf 
Mpayload * 150,000 lbs 
Nitncturt * 233*247 lb*' 


H02 J% 2 " “fuel 

I ^02 /l %2 “ 6 

Mh£( 1,108.133)/(7) -138.304.71 lba-> - *4*. 828.29 lbs 


/•» - 71.2 lbs/ft* Ah. * 4-43 lbs/ft 1 
^ * *%/ O2 * 33.724.898 ft* 

Vqj - 13. 340.28 3 ft* 

OOCR3H TANK: E l lip s oid (similar to tuo saucers placed convex to one another) 
V - ytabe - jtfa 3 b because a - c 
b * 11.70 ft— > h « 2b • 23.40 ft 


HYDROGEN TANK: cylinder with ellipsoid convex upper cap and ellipsoid concave 


V « IT R*h - va J h 
h - 41.78 ft 


ORIGINAL' PAGE IS 

TOTAL HEIGHT CP TWO TANKS IN TANDEM: ° F POOR QUALITY 

hfcotal * a 5*40ft >41,78 ft - 83. 18 ft 


APPENDIX C 
C.6. 1030106 
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C.G. Location (Dry) 


lomoonent 

weight (lb) 

Engines 

28568 

Engine Compartment 

15250 

instrument Unit 

6000 

Fuel Section 

60095 

Interstage 

6000 

Upper Stage 

105334 

Payload 

150000 

Nose Cone 

12000 

Total C.6. Location (Dry) « 10471 feet 


C.G. Location (W»H 

Engines 

28568 

Engine Compartment 

15250 

Instrument Unit 

6000 

Liquid Oxygen 

949828 

Fuel Section 

60095 

Liquid Hydrogen 

... 158304 

interstage 

6000 

Upper Stage 

1 05334 

Payload 

150000 

Nose Cone 

12000 


Total C.G. Location (Wet) « 50.21 feet 


Location (ft) 
5 

, 6 
’ 13.5 

37.59 

81.68 

133.68 

133.68 
19418 


5 

6 

13.5 
26.7 
37.59 
59.29 

81.68 

1 33.68 

1 33.68 
19418 
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